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I.  Introduction 


Shielded  cables  have  been  used  extensively  on  aircraft,  ground  and  space- 
missile  systems  to  reduce  the  electromagnetic  coupling  (crosstalk)  between 
electrical  equipments  which  are  interconnected  by  wires.  The  wires  which 
interconnect  these  electrical  and  electronic  devices  are  generally  routed  in 
densely-packed,  cable  bundles.  The  unintentional  electromagnetic  coupling  or 
crosstalk  between  these  wires  may  be  of  sufficient  magnitude  to  degrade  the 
performance  of  the  equipments  which  the  wires  interconnect.  In  order  to  re¬ 
duce  this  level  of  crosstalk,  shielded  cables  and  twisted  pairs  of  wires  have 
been  employed. 

In  an  earlier  report  [1],  the  prediction  of  crosstalk  to  braided-shield 
cables  was  investigated.  The  braided-shield  cable  consists  of  a  circular, 
cylindrical  shield  which  is  composed  of  belts  of  wires,  interwoven  to  provide 
flexibility  and  an  interior  wire  (the  shielded  wire)  located  on  the  axis  of 
the  shield  as  shown  in  Fig.  1-1.  The  interior  wire  or  shielded  wire  is  of 
radius  rw  and  the  shield  has  radius  rs>  The  shield  thickness,  tg,  is  approx¬ 
imately  equal  to  the  diameter  of  the  wires  making  up  the  braid  which  have 
radius  r^,  i.e.,  tg  =  2r^.  The  shielded  cable  usually  has  an  overall, 
insulating  jacket  of  thickness  t_. ,  and  the  shield  is  woven  in  B  belts  of  wires 
with  a  weave  angle  of  0^  as  shown  in  Fig.  1-1.  Each  belt  contains  W  wires. 

If  the  length  of  the  shield  is  denoted  by  then  each  braid  wire  is  of 

length  *^s/cos  0  .  The  medium  internal  to  the  shield  and  surrounding  the 
w 

interior,  shielded  wire  is  a  dielectric  with  permittivity  e  =  £rev  and  perme¬ 
ability  u  3  Uv  where  the  permittivity  and  permeability  of  free  space  are 
denoted  by  z and  p  ,  respectively. 

In  terminating  a  braided-shieru  cable,  for  example,  at  a  connector,  the 


Braid  Wires  of 
Radius  rt 


braid  is  often  stripped  back  exposing  the  interior,  shielded  wire  as  shown 

in  Fig.  1-1.  The  braid  is  terminated  (usually  to  a  ground  plane  such  as  an 

aircraft  structural  skin)  via  another  wire  which  will  be  referred  to  as  the 

pigtail  wire.  These  exposed  terminal  sections  will  be  referred  to  as  pigtails 

throughout  this  report  although  this  use  of  the  term  is  not  standard.  The 

term  pigtail  is  sometimes  used  elsewhere  to  denote  the  pigtail  wire.  The 

length  of  the  left  pigtail  section  will  be  denoted  by  and  that  of  the 

pL 


right  section  will  be  denoted  by 


it,. 


„t>*  The  total  length  of 


the  cable  is  denoted  as  t  =  it  T  +  it  +  it 

^  pL  s  pR 

When  braided  shield  cables  are  terminated  in  connectors  in  this  fashion, 
the  shielded  wire  is  directly  exposed  over  the  pigtail  sections  to  other  wires 
in  the  cable  bundle  which  are  also  terminated  at  the  connector  as  shown  in 
Fig.  1-2.  If  it  is  desired  to  carry  a  shield  connection  through  the  connector, 
a  separate  wire,  the  pigtail  wire,  connects  the  shield  braid  to  an  additional 
connector  pin.  To  avoid  bunching  of  these  shields  at  the  connector,  the  pig¬ 
tail  sections  at  a  connector  are  generally  of  different  lengths.  The  various 
design  handbooks  [4,5,6]  seem  to  recommend  against  using  pigtail  sections.  It 
appears  to  be  common  nevertheless  [7,8],  In  fact  the  author  has  observed  pig¬ 
tail  sections  in  excess  of  6  inches. 

It  was  shown  previously  [1,2]  that  the  coupling  over  the  pigtail  sections 
may  constitute  the  dominant  coupling  mechanism  to  the  shielded  wire.  Even 
though  the  total  length  of  the  pigtail  sections,  #£  +  may  constitute 

only  a  minor  fraction  of  the  total  cable  length,  it  was  shown  in  [1,2]  that, 
depending  on  the  values  of  the  cable  termination  impedances,  the  contribution 
to  the  received  voltages  at  either  end  of  the  cable  (the  voltages  of  the 


shielded  wires)  due  to  the  pigtail  sections  can  be  larger  than  the  contribution 


3 


over  the  (much  longer)  shielded  section.  In  this  situation,  the  shield  simply 

serves  to  reduce  the  exposed  length  of  the  interior  wire  )  from 

p  L  pR 

what  it  would  be  if  no  shield  were  present  ).  Thus  the  shield  still  pro¬ 
vides  a  reduction  in  crosstalk.  However,  it  was  shown  in  [1,2]  that  if  the 
pigtail  sections  were  eliminated,  one  could  realize  an  additional  reduction  in 
crosstalk  of  as  much  as  30dB  (over  certain  frequency  ranges).  Thus  the 
effectiveness  of  the  shield  in  reducing  crosstalk  can  be  substantially  reduced 
by  the  pigtail  sections. 

An  additional  result  shown  in  [1,3]  was  that  the  crosstalk  to  or  from  a 
braided-shield  cable  including  the  effects  of  pigtail  sections  could  be  pre¬ 
dicted  quite  accurately  with  the  multiconductor  transmission  line  model  (MTL) 
whose  basic  features  are  described  in  [9].  Certain  other  features  of  the 
coupling  which  are  distributed  parameter  in  nature  could  not  be  predicted  with 
lumped  circuit  models.  For  example,  it  was  shown  in  [1,2]  that  the  coupling 
to  a  shielded  cable  in  which  only  one  end  is  grounded  depended,  in  some  cases, 
very  strongly  on  which  end  of  the  shield  was  grounded.  This  is  clearly  a  dis¬ 
tributed  parameter  phenomenon,  and  the  MTL  model  predicted  this  result  quite 
accurately.  In  addition,  it  was  shown  in  [2]  that  the  coupling  to  a  shielded 
wire  may  depend  quite  strongly  on  the  orientation  of  the  pigtail  wires.  Simple 
lumped  models  failed  to  predict  these  effects,  whereas  the  MTL  model  predicted 
these  results  quite  accurately. 

In  order  to  describe  the  coupling  through  cable  shields  we  first  consider 
a  cable  having  a  solid  shield.  The  coupling  from  an  exterior  field  to  the 
interior,  shielded  wire  for  solid-shield  cables  is  described  by  a  surface 
transfer  impedance,  z^,  per-unit-of  cable  length.  For  an  infinitely  long, 
solid,  conducting  cylinder,  suppose  the  cylinder  carries  a  total  current  I 


5 


.  The  return  path  for 


directed  in  the  axial  (x)  direction  where  1=1.  +  I 

in  out 

I.  is  within  the  cylindrical  surface  while  the  return  path  for  I  is 
in  r  out 

outside  the  surface.  The  conducting  cylinder  has  finite  conductivity  a  and 
therefore  the  current  I  flowing  along  the  cylinder  will  induce  electric  fields 
on  the  inner  and  outer  surfaces  of  the  cylinder,  and  respectively, 

which  are  directed  in  the  axial  direction  of  the  cylinder.  The  currents  and 
induced  electric  fields  may  be  related  as  [10-13,  19,  22] 


E. 

in 

z.  . 

IX 

•1 

z . 

10 

r~  “i 

i. 

in 

E  . 
out 

z  . 

OI 

z 

oo 

*out 

V/m 


(1-D 


The  terms  z ^  and  zqq  may  be  thought  of  as  per-unit-length  self  impedances 


of  the  surfaces: 


z .  . 
ii 


ft/m 


z 

oo 


E  . 
out 

^out 


ft/m 

I.  =0 
m 


(l-2a) 


(l-2b) 


whereas  the  terms  z.  and  z  .  are  called  surface  transfer  impedances 

io  oi 

(per-unit-length) : 


z 


io 


z  . 
oi 


out 


E  „ 
out 

I. 

in 


ft/m 

I.  =  0 

in 


ft/m 


I  =0 
out 


(l~3a) 


(l-3b) 


These  surface  transfer  impedances  relate  the  current  on  one  surface  of  the 
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cylinder  to  the  induced  field  on  the  other  surface. 

The  two  surface  transfer  impedances  z.  and  z  .  will  be  assumed  to  be 

10  01 

equal  and  will  be  designated  as  z ,  i.e.. 


z 


T 


z , 

10 


z  . 
01 


(l-4a) 


The  self- impedances  of  the  shield, 
same  and  designated  as  z  ,  i.e.. 


z..  and  z  will  also  be  taken  to  be  the 
11  oo 


z 


S 


z 

oo 


(l-4b) 


Certainly  for  cylinders  with  wall  thicknesses  which  are  sufficiently  small, 
these  t  '.1  be  reasonable  approximations.  Thus  (1-1)  becomes 


(1-5) 


The  surface  transfer  impedance  accounts  for  skin  effect  -  the  tendency 
for  a  current  to  concentrate  on  a  conductor  surface  nearest  its  source.  At 
D-C,  any  current  flowing  in  the  shield  will  tend  to  be  uniformly  distributed 

over  the  shield  cross  section.  As  the  frequency  is  increased  the  current  tends 

to  concentrate  towards  the  shield  surfaces.  Thus  for  increasing  frequencies 
less  current  diffuses  through  the  wall  to  induce  an  electric  field  on  the 

opposite  surface.  In  the  limit  as  the  frequency  is  increased  without  bound 

there  is  perfect  isolation  between  the  inner  and  outer  walls  of  the  shield. 

If  the  frequency  is  sufficiently  small  so  that  the  shield  thickness  is  less 
than  a  few  skin  depths,  the  surface  transfer  impedance  is  approximately  equal 
to  the  impedance  of  the  shield.  In  fact,  the  two  should  converge  as  the  fre- 


*; 
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quency  is  reduced  to  zero.  We  will  often  refer  to  this  surface  transfer 


impedance,  z^,  for  solid  shields  as  the  diffusion  impedance  and  denote  it  as 

z,  =  for  the  above  reasons, 
d  T 

The  surface  transfer  impedance  represents  the  coupling  of  some  field 
external  to  the  cylinder  to  the  field  internal  to  the  cylinder  and  occurs  via 
diffusion  through  the  finitely  conducting  shield.  This  surface  impedance  is 
employed  to  predict  the  coupling  from  some  field  incident  on  a  coaxial  cable 
to  the  loads  connected  to  the  end  points  of  the  coaxial  cable  in  the  follow¬ 
ing  manner  [19].  As  an  example,  consider  the  coaxial  cable  above  the  ground 
plane  shown  in  Fig.  l-3£a)  The  coaxial  cable  consists  of  a  conducting  cylinder 
and  a  concentrically-located,  conducting  wire.  An  incident  field,  such  as  a 
uniform  plane  wave,  illuminates  the  cable  and  induces  a  current,  Iq  ,  flowing 
along  the  shield  and  returning  through  the  ground  plane.  For  the  purposes  of 
computing  this  induced  current,  Iout>  it  is  universally  assumed  that  the 
interior  of  the  shield  and,  in  particular,  1^  have  no  effect  on  the  external 
circuit  [19].  In  other  words,  Iout  is  traditionally  calculated  as  simply  the 
current  induced  on  an  isolated  cylinder  above  ground.  Thus  we  assume  a  uni¬ 
lateral  effect  -  outside  to  inside  -  to  simplify  the  calculations.  Then  the 
effect  of  this  induced  current,  Iq  ,  and  consequently  the  incident  field,  on 
the  internal  structure  is  manifested  as  a  per-unit-length  voltage  source  in 
the  cable  interior  equal  to  z  iQut  (V/m) .  The  equivalent  circuit  represent¬ 
ing  a  small,  Ax  section  of  the  shield  and  interior,  shielded  wire  is  shown  in 
Fig.  l-3(b). 

The  per-unit-length  impedance  of  the  interior  wire  is  represented  by  z^ 
and  the  per-unit-length  inductance  and  capacitance  of  this  Ax  section  of  the 
line  are  represented  by  l  and  c,  respectively.  The  current  I  produces  a 
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Figure  1-3.  Field-to-wire  coupling  involving  shielded  cables. 


voltage  drop  z  Ax  I .  along  the  interior  surface  of  the  shield  as  indicated 
by  (1-1)  and  (l-2a),  and  the  effect  of  the  external  shield  current,  I  , 
induced  by  the  incident  field  is  represented  by  a  voltage  source  z ^  Iq  Ax  as 
indicated  by  (1-1)  and  (l-3a) . 

From  the  equivalent  circuit  in  Fig-  1-3 (b),  we  may  derive  the  transmission 
line  equations  as  [9] 

:r(x)  =  -  (z  +  z  +  juiJl)  I.  (x)  +  z  I  (x)  (l-6a) 

dx  w  s  J  in  T  out 

dI  (x) 

"~f~  =  -jwc  V  (x)  (l-6b) 

Note  that  the  effect  of  the  incident  field  appears  in  these  transmission  line 
equations  as  a  driving  term,  zTIQut* 

This  is  an  example  of  the  use  of  the  surface  transfer  impedance  for 
solid,  cylindrical  shields.  For  braided  shields,  additional  coupling  occurs 
via  the  small  holes  in  the  shield  introduced  by  the  overlapping,  wire  braid. 
This  construction  of  the  shield  braid  introduces  small,  diamond-shaped  holes 
between  the  belts  of  wires.  These  holes  allow  other  coupling  mechanisms  to 
occur  between  the  outside  environment  and  the  interior,  shielded  wire  which 
were  not  present  for  a  solid  shield.  For  the  solid  shield,  the  coupling  from 
the  exterior  to  the  interior  occurred  only  by  diffusion  through  the  metal. 

For  the  braided  shield,  additional  inductive  and  capacitive  coupling  occur  via 
the  holes  in  the  shield  as  discussed  by  Latham  [15,16]  and  Vance  [17]. 

Latham  provided  the  following  development.  Suppose  the  braided  shield 
and  interior  wire  are  perfectly  conducting.  Consider  a  small  Ax  section  of 


the  cable  as  shown  in  Fig.  1-4.  Select  a  planar  surface  of  length  Ax  between 
the  interior  wire  and  the  shield  as  shown  in  Fig.  1-4 (a).  Faraday's  law 
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written  around  this  contour  is 


c  s 


=  -jwi|> 


Write  the  flux  through  che  surface  as 


i|i  =  i  Ax  1^  -  Ax  Ig  (1-8) 

where  I  is  the  total  current  carried  by  the  interior  wire  and  I  is  the  total 
c  s 

current  carried  by  the  shield.  Defining  the  voltage  between  the  interior  wire 
and  the  shield  as 


wire 


shield 


(1-9) 


we  obtain 


V (x  +  Ax)  -  V (x)  ■  -  jw(£  Ax  I  -  £„Ax  I  )  (1-10) 

cels 

Dividing  (1-10)  by  Ax  and  taking  the  limit  as  Ax  ->  0  yields 


=  -  ju>£  Ic(x)  +  jW£T  Is(x) 


(1-11) 


£  =  -JJ-  £n  (— )  (1-12) 

4tt  r 

w 

is  the  per-unit-length  self  inductance  of  the  cable  where  y  is  the  perme¬ 
ability  of  the  interior  medium.  The  quantity  £T  is  the  per-unit-length 
transfer  inductance  which  is  zero  for  solid,  circular,  cylindrical  shields 
which  have  the  interior  wire  on  the  axis  of  the  shield.  Holes  in  the  shield 
as  well  as  shields  having  a  conducting  direction  not  along  the  x  axis  (such 
as  tape-wound)  cause  £^  to  be  nonzero. 
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Similarly,  consider  Fig.  1-4 (b).  A  closed  surface  is  constructed  around 
the  interior  wire.  Writing  Ampere's  law  for  this  surface 


^  H  •  d!  =  joje  f  1  • 


(1-13) 


we  obtain 


Ic(x  +  Ax)  -  Ic(x)  =  jui  Qc(x) 


(1-14) 


where  Qc  is  the  total  charge  contained  on  the  interior  wire  over  Ax.  The 
charge  on  the  interior  wire  and  on  the  shield,  Qg,  can  be  related  to  the 


voltage  as 


V  Ax  =  s  Q  -  s„  Q 
c  T  xs 


Solving  (1-15)  for  and  substituting  into  (1-14)  we  obtain  as  Ax  ->  0, 


(1-15) 


dlc(x) 


-  ju^V(x)  -  jw  ~  Qs(x) 


(1-16) 


For  a  circular,  cylindrical  shield 


3  "  2^  </) 
w 


(1-17) 


The  transfer  elastance  sT  is  due  to  the  penetration  of  the  electric  field 
lines  through  the  shield.  The  terms  £  and  s^  are  determined  for  braided 
shield  cables  in  terms  of  the  magnetic  and  electric  polarizabilities  of  the 
holes  in  the  braid  [16]. 

Vance  [16]  has  chosen  to  include  the  diffusion  term  in  the  voltage  change 
expression  for  braided  shields  as  was  present  for  solid  shields.  For  a 


braided  shield,  Vance  assumed  that  all  strands  of  the  braid  were  connected, 


electrically,  in  parallel.  For  example,  if  there  are  B  belts  of  wires,  W  wires 


per  belt  and  the  braid  is  woven  with  a  weave  angle  Q  ,  then  the  per-unit-length 
D-C  resistance  of  the  braid  is 


rDC"  ir'r'  z'  o'  B  W  cos©  n,/m 
b  w 


(1-18) 


where  r^  is  the  radius  of  each  strand  and  a  is  the  strand  conductivity.  Vance 
presumes  that  this  braid  impedance  is  modified  with  increasing  frequency  in 
the  same  manner  as  the  solid  cylinder.  Thus  the  per-unit-length  diffusion 
impedance,  z^,  is  taken  by  Vance  to  be 


yd 

Zd  rDC  sinhyd 


(1-19) 


where  d  is  the  diameter  of  the  braid  wires,  d  =  2r^,  and  y  =  (1  +  j)/6  where 
-1/2 

6  =  (Trfpa)  is  the  skin  depth.  The  self  impedance  of  each  strand  of  the 

braid,  z  ,  is  chosen  here  to  be  computed  as  the  impedance  of  an  isolated  round 
s 

wire.  The  total  self  impedance  of  the  braid  is  then  chosen  as  the  impedance 

of  all  strands  in  parallel. 

Comparing  Fig.  1-3  and  Fig.  1-4 (a)  we  note  that 

I  -I.  (l-20a) 

c  in 

I  =  I  -I. 
s  out  in 


or 


I 


I  =  I  (l-20b) 

in  c 

I  =  I  +  I 
out  s  c 

Substituting  (l-20b)  into  (1-6)  and  adding  the  transfer  inductance  and  transfer 
elastance  we  obtain 
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?Tl 


dV  (x) 
dx 


(z.w  +  Zg  +  jco«.  -  zd>  Ic(x)  +  (zd  +  joj«-T)  Ig(x)  (l-21a) 


dlc(x) 

dx 


-  jwc  V (x)  -  jiuc  s  Qs(x) 


(l-21b) 


where  the  diffusion  impedance  is  zd  and  the  total  transfer  impedance  for 
the  braided  shield  now  becomes 

ZT  =  Zd  +  U//m  (1-22) 


The  purpose  of  this  report  is  to  extend  the  above  concepts  to  multi¬ 
conductor  cables  in  which  the  individual  wires  may  be  unshielded  or  shielded 
with  solid  or  braided  shields.  The  results  will  be  incorporated  into  a 
digital  computer  program,  SHIELD,  which  may  be  used  to  predict  crosstalk  in 
cables  which  have  shielded  wires.  A  description  of  the  program  is  given  in 
Section  III.  A  User's  Manual  is  given  in  Section  IV.  The  shields  may  have 
pigtail  sections  as  described  previously.  The  pigtail  section  lengths  need 
not  be  the  same.  Two  choices  of  reference  conductor  for  the  line  voltages 
are  available.  The  cables  may  be  over  a  ground  plane  or  within  an  overall, 
circular  cylindrical  shield. 

The  distributed  impedances  of  all  conductors  due  to  their  finite  con¬ 
ductivities  are  included  in  the  model.  This  allows  the  modeling  of  common 
impedance  coupling  as  is  illustrated  in  Fig.  1-5.  Consider  two  wires  and 

a  reference  conductor  shown  in  Fig.  l-5(a).  At  a  sufficiently  small  frequency, 
the  current  in  the  generator  wire  can  be  computed  as 


IG  "  Z 


og  +  zj£g 

The  majority  of  this  current  passes  through  the  reference  conductor  pro¬ 


d-23) 


ducing  a  voltage  drop  across  the  reference  conductor  total  impedance  Z^,  of 
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©  Generator  Wire  IG 


©  Receptor  Wire 


- u _ 

^Zor  Nor 

V«4R  | 

-  VC|  + 
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Reference  Conductor 


Net 

Received 

Voltage 


Common  Impedance 
Coupling  Contribution 


>1 _ 

X7  \  Vc  ommon 

/  Electromagnetic  Coupling 
Coupling  Contribution 
Assuming  Perfect  Conductors 

Frequency 

(b) 


Figure  1-5.  Common  impedance  coupling. 


(1-24) 


vc 


Vg 


This  common  impedance  voltage  is  then  coupled  directly  to  the  receptor 
circuit  producing  voltages 


J0R 


0R  Z0R  + 


*R 


'iR  '  "  ZOR  +  Z*R  VC 


(1-2 5a) 


U-25b) 


These  low-frequency  voltages  produce  a  "floor"  in  the  low-frequency  cross¬ 
talk  as  shown  in  Fig.  1-5 (b). 

Each  cable  shield  may  be  grounded  to  the  reference  conductor  (via 
the  pigtail  wires)  at  one  end,  both  ends  or  neither  end.  The  lengths  of 
each  pigtail  section  on  each  shielded  cable  may  be  chosen  by  the  use*.  For 
example,  a  shield  may  have  no  pigtail  at  one  end  and  a  pigtail  at  the  other. 

In  section  II,  the  derivation  of  the  line  equations  is  presented. 
Comparisons  of  the  code  predictions  with  experimental  data  are  provided  in 
Section  V.  A  listing  of  the  code  is  given  in  the  Appendix. 
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II.  Derivation  of  the  Line  Equations 

The  basic  technique  is  illustrated  in  Fig.  2-1.  A  line  is  shown  for 
illustration  in  Fig.  2-1  (a)  as  having  2  unshielded  wires  and  2  shielded  cables 
having  pigtail  sections  of  various  lengths.  The  line  is  divided  into  sections 
of  uniform  cross  section.  The  chain  parameter  matrices  (which  in  this 

case  are  12  x  12)  are  obtained  for  each  section  and  the  overall  chain  para¬ 
meter  matrix  of  the  line,  4>,  is  determined  as  the  product  (in  the  appropriate 
order)  of  the  chain  parameter  matrices  of  the  individual  sections  as  described 
in  [1].  Once  the  overall  chain  parameter  matrix  for  the  entire  line  is  de¬ 
termined,  the  two  terminal  constraints  are  incorporated  and  the  terminal 
voltages  (and/or  currents)  of  the  line  are  determined.  This  is  the  essential 
technique  involved  in  the  program. 

2 . 1  Shielded  Section 

First  consider  the  case  of  u  unshielded  wires  and  s  shielded  cables  (with 

no  pigtail  sections)  as  shown  in  Fig.  2-2.  The  line  axis  is  the  x-coordinate 

axis  and  the  phasor  conductor  currents  for  sinusoidal  steady-state  excitation 

of  the  line  are  defined  in  the  +x  direction.  The  current  in  the  reference 

conductor  is  the  negative  sum  of  these  U  +  2S  conductor  currents.  U  conductor 

currents,  I  are  defined  for  the  unshielded  wires  am.  2S  currents,  I- .  and 
wi  wi 

I  are  defined  for  the  shielded  cables  as  shown  in  Fig.  2-2.  The  phasor 
si 

voltages  of  the  U  +  2S  conductors,  V  V.,  and  V  are  defined  with  respect 
°  wi  Oi  si 

to  the  reference  conductor. 

To  facilitate  deriving  the  line  equations,  it  is  convenient  to  redefine 
the  voltages  and  currents  of  the  shielded  cables  as  shown  in  Fig.  2-3.  The 
shielded  wire  voltage  is  defined  as 


(b)  : 

Figure  2-1.  Illustration  of  the  basic  computational  technique  for 

incorporation  of  pigtail  sections.  / 

> 
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//////// 

Reference  Conductor 
(Ground  Plane  or  Overall  Shield) 

U  S 

X  Iwj  +  Y  (X*i  +  ) 

i  =  1  i  =  I 


Figure  2-2.  Definition  of  conductor  currents  and  voltages 


2 


E 


u  s 

2  Iwi  +  2  ( I s i  +  ^wi  ) 


i=l  i *  I 


Figure  2-3,  Derivation  of  the  voltage  change  equation  for  a  shielded  cable. 
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and  the  voltage  of  the  shield  remains  unchanged.  The  shielded  wire  currents 


are  redefined  such  that  a  current  in  the  shield  has  I  .  and  I».  with  a  return 

SI  0)1 

exterior  to  the  shield  and  I*  with  a  return  interior  to  the  shield. 

wi 

Defining  a  contour  as  shown  in  Fig.  2-3  between  the  shield  and  the 
reference  conductor  and  between  x  and  x+Ax  and  writing  Faraday's  law  around 
this  contour  we  obtain,  in  the  limit  as  Ax  0  [1], 


dV  . 

— =  -  ju  f  .  -  z  .  (I  .  +  I'.)  +  z  , .  (I~.) 

dx  si  si  si  m  di  oil 


(2-2) 


U 


+  jwJl_.  I-.  -  zn  (Z  I  .  +  S  (I  .  +  I-.)} 

J  Ti  oil  0  .  .  o)i  .  '  si  o)i 

1=1  1=1 


where  tp  ■  is  the  per-unit-length  magnetic  flux  penetrating  the  surface 
bounded  by  the  contour  due  to  all  currents  in  the  system.  This  is  determined 
by  assuming  the  conductors  are  perfect  conductors  and  the  shields  are  solid 
with  no  penetrations.  The  other  per-unit-length  quantities  on  the  right-hand- 
side  of  (2-2)  are  the  self  impedance  of  the  shield,  z^,  the  diffusion 
impedance  of  the  shield  z^,  the  self  impedance  of  the  reference  conductor  z^ 
and  the  transfer  inductance  of  the  braided  shield  (due  to  holes  in  the  braid) 

Jl  ..  If  the  shield  is  solid,  is  removed  and  z  .  and  z,.  are  changed. 

Ti  ’  Ti  si  di  e 

The  voltage  change  expression  for  the  shielded  wire  voltage  (with  respect 
to  the  shield),  V-  ,  is  obtained  in  a  similar  fashion: 

U)1 


dV-  . 

— =  -jw  1|)  .  -  Z-.  I/'.  -  z  .  I-.  +  z  I~.  +  z  I  . 
dx  u)i  o)i  wi  si  o)i  di  o)i  di  si 


j^Ti  ^i  +  j^Ti  lsi 


(2-3) 


Here  ij)~  is  the  per-unit-length  magnetic  flux  penetrating  a  similar  contour 


conductor)  as 


(2-9b) 


(2-9c) 


Then  equations  (2-2),  (2-4)  and  (2-5)  become 


v  =  -  jw  i  -  zcl  -  z0  y  I 


(2-10) 


where  the  dot  denotes  ordinary  derivative  with  respect  to  x  and  the  (l!  +  2S)  x 


fU  +  2S)  matrix  Z  is 
„c 


Zc=  ° 


(Zs  -  Zd  -  jmLT) 


0  (Z  -  Z ,  -  juL  )  [2 (Z  -  z  -  jmT)  +  Z~] 
\s  _d  J  _T  _s  _d  „T  _w 


(2-11) 


The  submatrices  in  (2-11)  are  diagonal  with  the  main  diagonal  entries  given  by 


[  Z  ].. 

_W  11 


i  =  1, - ,  U 


(2-12a) 


[  Z  ].. 

_s  11 


i  =  1,  — ,  S 


(2-12b) 


1  Zd  hi 


i  =  1, - ,  S 


(2-12c) 


[  Z-  ].. 

w  11 


i  =  1, - ,  S 


(2-12d) 


l  h  ]ii 


i  =  1, - ,  S 


(2-12e) 


and  all  other  entries  are  zero. 
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We  denote  the  entry  in  the  i-th  row  and  j-th  column  of  a  matrix  M  by 
[  M  ]  .  The  matrix  U  is  the  (U  +  2s)  x  (U  +  2S)  unity  matrix  with  ones  in 
every  position: 


111  -  1 


111 

,  i  ' 

i  i  t 

•  1  • 

i  ‘  i 


1 

i 

i 

i 

i 


(2-13) 


111 


Once  again,  if  the  i-th  shield  is  solid,  is  removed  from  L_  and  z,.  and 

1 1  ~  r  di 

are  computed  for  a  solid  shield. 

The  (U  +  2S)  vector  ^  relates  the  per-unit-length  magnetic  fluxes  passing 
between  each  conductor  and  the  reference  conductor  to  the  line  currents  via  the 
(U  +  2S)  x  (U  +  2S)  per-unit-length  inductance  matrix  L: 


±  -  L  I  (2-14) 

where 

[  L  ]  =  £  (2-15) 

for  i,j  =  1,  - ,  U+2S. 

Computation  of  the  entries  in  L  are  described  in  [1]  and  [24].  The  self 
inductance  of  a  shield  is  computed  as  though  the  shield  were  a  solid  wire. 

The  mutual  inductances  between  a  shield  and  an  unshielded  wire  are  similarly 
computed.  The  only  variation  from  this  obvious  scheme  is  the  computation  of 
the  mutual  inductance  between  a  shield  and  its  shielded  wire.  This  is  equiv¬ 
alent  to  the  self  inductance  of  the  shield  as  described  in  [1].  The  compu¬ 
tation  of  the  per-unit-length  transfer  inductances  is  described  in  [22,23], 

The  primary  definition  of  the  per-unit-length  inductance  matrix  is 
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(2-16) 


where  is  a  vector  of  per-unit-length  magnetic  flux  passing  between  each 
wire  and  the  reference  conductor,  is  a  vector  of  per-unit-length  magnetic 
flux  passing  between  each  shield  and  the  reference  conductor,  and  is  a 
vector  of  magnetic  flux  passing  between  each  shielded  wire  and  the  reference 
conductor.  From  the  definition  in  (2-16)  one  can  easily  obtain  the  entries 
in  L  [1].  With  these  definitions,  the  voltage  change  expression  in  (2-10) 
becomes 

V(x)  =  -  (Zc  +  zq  U  +  ju  L)  J.(x)  (2-17) 

It  now  remains  to  determine  the  second  transmission-line  equation  -  the 
current  change  expression. 

The  current  change  expressions  are  similarly  computed.  Essentially  we 

need  to  relate  the  voltage  of  each  conductor  to  the  per-unit-length  charge 

carried  by  that  conductor.  The  per-unit-length  charges  and  conductor  voltages 

are  defined  in  Fig.  2-4.  These  definitions  are  directly  analogous  to  those 

for  the  line  voltages  and  currents  in  Fig.  2-3.  The  shield  carries  per-unit- 

length  charges  q  .  and  q~.  on  its  outer  surface  and  -q-.  on  its  inner  surface. 

6  Msi  Mwi  Wi 

Voltages  V-.,  V« .  and  V  .  are  defined  as  in  Fig.  2-3. 
wi  wi  si 

The  basic  relationship  is: 
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The  current  change  expression  is  obtained  by  inverting  (2-20)  as 


5  =  s"1  V 


(2-21) 


=  C  V 


where  we  denote  the  per-unit-length  capacitance  matrix  as 


C  =  S 


(2-22) 


Multiplying  both  sides  of  (2-21)  by  ju  we  obtain 


jojg_(x)  =  jwC  V(x) 


(2-23) 


Defining 


I(x)  =  -  jw^(x) 


(2-24) 


we  obtain  the  current  change  expression 


I.(x)  =  -jioC  V(x) 


(2-25) 


The  holes  in  the  braid  of  a  braided  shield  cable  introduce  additional 
terms  in  S.  We  may  treat  the  unshielded  wires  and  shields  of  the  cable  as  a 
set  of  U  +  S  wires  in  a  homogeneous  medium  (free  space  here)  so  that  [1] 


S  S 
_ ww  ,ws 


S  S 
_sw  _ss 


L  L 
„ww  ~ws 


L  L 
„sw  „ss 


(2-26) 


which  is  symmetric,  i.e.,  L  =Lt,L  =Lt,L  =Lt  and  L  =  L  t  where 

_ww  „ww  _ss  _ss  _ws  ,ws  „sw  _ws 

we  denote  the  transpose  of  a  matrix  M  by  MC.  That  is,  the  upper  left  block 
of  L  in  (2-16)  are  computed  as  though  the  shields  comprise  a  set  of  isolated 


wires  and  q_~  =  0  in  (2-20) 
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In  order  to  determine  the  remaining  entries  in  S,  consider  the  definition 
of  S  in  (2-20)  and  the  voltage  and  charge  definitions  in  Fig.  2-4.  In  (2-20) 
setting  ^  =  0  we  obtain  from  the  first  equation 


V  =  S  ~  q^ 

— W  „WW  <'*W 


%  =  =  2 


From  Fig.  2-4,  we  clearly  have 


(2-27) 


S  -  =  S 
_ww  ,ws 


(2-28) 


Similarly 


V  =  S  -  q~ 
— s  _sw  -V 


q  =  q 


so  that  from  Fig.  2-4  we  have 


(2-29) 


s  -  =  s  +  s„ 

.sw  ,.ss  „T 


(2-30) 


The  matrix  ST  is  diagonal  and  has  the  per-unit-length  transfer  elastances 


sT^  as  determined  by  Latham  [15,16]  on  its  main  diagonal: 
^T^ii  =  STi 


(2-31a) 


LVij  =  ° 


(2- 31b) 


For  the  last  block  row  of  S  we  determine  the  entries  S  S'  ,  with 

_ww’  _ws  _ww 

the  definitions  in  Fig.  2-4.  For  example,  in  terms  of  V^, 


V'  =  S'  q  +  S*  q  +  q* 
— w  _ww  -Hw  ..  ws  -^s  _ww  -^w 


(2-32) 


V-  =  V  ~  +  V 
— v  — w  — s 


(2-33) 
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So  it  remains  to  determine  S'  ,  S'  ,  S''.  Setting  q  =  q.  =  0  in  (2-32)  we 

_ww  .ws  .  ww  °  s  -n#  — 


V'  =  S'  q 
— w  ~  ww  -Hj 


(2-34) 


=  ^  =  0 


S'  =  0 

VWW 


(2-35) 


(Actually  S'  is  not  zero  due  to  penetrations  for  braided  shields.  But  when 

WW 

added  to  Sgw  according  to  (2-33),  it  is  felt  that  this  will  be  negligible 
compared  to  the  entries  in  Sgw  and  is  accordingly  neglected.)  Now 


V-  =  S'  q 
--W  „ws  -‘■s 


(2-36) 


\  %  =  %  =  ^ 

Here  penetrations  through  braid  holes  probably  should  be  included  so  that 
[15,16] 


sr  =  ST 

.ws  _T 


(2-37) 


Similarly, 


v'  =  s'' 

— w  .ww 


I  %  -  %  =  4 

But  this  is  simply  the  interior,  self-elastances  of  each  shielded  cable  [1]: 


fcn(r  ,/r'.) 

rS''  1  = 

~ww  ii  2i re. 

l 


(2-39a) 


[§"]••  =  o 

.ww  1J 

l*j 


(2-39b) 


where  r  ,  is  the  interior  shield  radius,  r'.  is  the  radius  of  the  shielded 
si  wi 
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wire  and  is  Che  interior  dielectric  permittivity.  Define 


c  =  C  „ 

.INT  .ww 


With  the  above  results  and  (2-33)  we  obtain 


S~  =  S  +  S~ 

_ ww  .sw  .ww 


=  S 
.sw 


S-  =  S  +  S- 
„ws  .ss  „ws 


=  +  ST 

_  SS  „  I 


S  ~  ~  =  S  *  4-  S  ~  a 
,ww  „sw  .ww 


!sw  +  ^INT 


~  s  +  s_  +  s 

.as  „T  _ INT 


(2-40) 


(2-41a) 


(2-41b) 


(2-41c) 


S  =  S 


<?ss  +  5t> 


(2-42) 


?ws  <?ss  +  ?T>  (?SS  +  ?T  +  ~INT)” 


In  order  to  show  that  this  result  is  reasonable,  suppose  all  shields 
are  solid  so  that  S_  =  0. 


Then 


In  order  for  the  results  to  be  reasonable  for  solid  shields,  the  structure  of 
the  per-unit-length  capacitance  matrix  C  must  be  [1] 


where 


(2-45) 


is  diagonal.  The  matrix  in  (2-44)  can  be  inverted  in  block  form  using  the 
results  in  [25].  Matching  blocks  in  (2-43)  we  obtain 
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2 . 2  Pigtail  Sections 

If  a  pigtail  occurs  over  a  section  of  line  as  illustrated  in  Fig.  2-1, 
certain  entries  of  L  and  Zcin  (2-17)  and  S  in  (2-20)  associated  with  the 
shield  attached  to  the  pigtail  and  its  shielded  wire  are  computed  differently. 
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With  regard  to  the  entries  in  L,  the  pigtail  wire  is  simply  treated  as  another 

wire  instead  of  a  shield.  Thus  those  entries  associated  with  the  shielded 

wire  and  the  shield  in  L  ,  L  - ,  L--,  are  simply  computed  as  th'.  ugh  the 

_ss  _SW  _WW  r  J  r  a 

shielded  wire  and  the  pigtail  wire  are  ordinary  wires  and  one  does  not  enclose 
the  other.  The  self  impedances  of  the  shield  in  Zg  are  now  computed  for  the 
pigtail  wire,  and  the  corresponding  entries  in  and  L^,  are  zero. 

A  similar  process  applies  to  the  entries  in  the  per-unit-length  elastance 
matrix  S  in  (2-20).  Those  entries  in  S  - ,  S  - ,  associated  with  the  pig- 

_  ^WW  ^  SW  ^  WW 

tail  wire  and  its  shielded  wire  are  computed  as  though  these  were  simply  two 
wires  with  neither  one  enclosing  the  other. 

2 . 3  The  Chain  Parameter  Matrices  of  the  Sections 

Consider  a  typical  uniform  section  of  the  line  of  length^.  The  trans¬ 
mission-line  equations  are 


V (x)  =  -  (Z  +  z  U  +  jo  L)I(x)  (2-48a) 

—  „c  o  _  — 

I(x)  =  -  joo  C  V (x)  (2-48b) 


we  may  compute  the  chain  parameter  representation  of  this  uniform  section  of 
line  [1]: 


— (XR} 


hi  h2 

"v(xL)  ‘ 

hi  -22, 

_I(V- 

(2-49) 


where  x  is  the  coordinate  of  the  right  end  of  the  section,  x  is  the 
R  L* 

coordinate  of  the  left  end  of  the  section  and  the  section  length  is 


(2-50) 
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We  first  compute  the  eigenvalues  and  eigenvectors  of  the  product  Y  Z 


where 


Y  =  juiC 


Z  =  Z  +  z  U  +  1u)L 
~c  o_ 


(2-51a) 

(2-51b) 


The  (U  +  2S)  eigenvectors  are  arranged  in  a  (U  +  2S)  x  (U  +  2S)  matrix  T  such 
that 


-1 


T  Y  Z  T  =  y‘ 


(2-52) 


where  y2  is  a  (U  +  2S)  x  (u  +  2s)  diagonal  matrix  containing  the  eigenvalues 
of  Y  Z  on  its  diagonal . 

The  entries  in  $  are  then  determined  from  [9]  as 


?11 


\  Y  1  T  [e?*  +  e  T  1  Y 


1  „-l 


:i2 


=  -  2  l  Ty  [ 


^  1  T  1 


1  t 


*21 


-  \  T  ^  -  e'l^j  y  1  T'1  Y 


122 


2  T  [  e 


+  e“T  t"1 


(2-53a) 

(2-53b) 

(2-53c) 

(2-53d) 


The  (U  +  2S)  x  (U  +  2S)  diagonal  matrix  y  is  the  square  root  of  y2  and  the 


(U  +  2S)  x  (U  +  2S)  diagonal  matrices  and  e  have  entries  e  1  and 

-Yi* 


on  the  main  diagonals,  respectively. 

2.4  The  Overall  Chain  Parameter  Matrix  of  the  Line 


Once  the  chain  parameter  matrices  of  the  sections  are  determined,  the 
overall  chain  parameter  matrix  of  the  line  is  the  product  (in  the  appropriate 
order)  of  the  individual  chain  parameter  matrices  of  the  sections  as  is 
illustrated  in  Fig.  2-1  [1], 


v '  -1 
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The  final  step  in  this  process  is  the  incorporation  of  the  terminal 
conditions  at  the  two  ends  of  the  line  into  the  overall  chain  parameter  matrix 
of  the  line  of  total  length  in  order  to  solve  for  the  terminal  voltages 
and/or  currents  of  the  line.  The  general  terminal  configuration  is  shown  in 
Fig.  2-5.  The  terminal  configuration  at  the  right  end,  x  =X ,  is  similar. 

Each  unshielded  wire  and  each  shielded  wire  is  connected  to  the  reference  con¬ 
ductor  and  each  other  by  an  admittance.  A  current  source  is  also  attached  to 
each  of  these  wires  (and  directed  from  the  reference  conductor  to  the  wire) . 
Each  shield  may  be  connected  to  the  reference  conductor  via  a  short  circuit 
or  an  open  circuit. 

In  order  to  incorporate  the  terminal  conditions,  we  first  rearrange  the 
overall  chain  parameter  matrix  of  the  line.  At  this  stage,  the  overall  chain 
parameter  matrix  of  the  line  (which  is  the  product  of  the  individual  chain 


Rearrange  the  rows  and  columns  of  this  matrix  such  that 


In  each  vector,  the  first  U  +  S  entries  contain  the  unshielded  wire  and 

shielded  wire  voltages.  The  next  U  +  S  entries  contain  the  unshielded  wire 

and  shielded  wire  currents.  The  next  S  entries  contain  the  shield  voltages, 

V  ,  for  those  shield  terminations  which  are  open  and  the  shield  currents,  I  , 
— s 

for  the  remaining  shield  terminations  which  are  shorts.  The  next  S  zero 
vectors  are  constraints  imposed  by  the  zero  shield  voltages  for  those  shields 
which  are  shor t-circuited  and  the  remaining  shield  currents  for  those  shields 
which  are  open  -  circuited . 


generalized  Norton  equivalent  representation  relating  the  unshielded  wire  and 
shielded  wire  currents  and  voltages: 


i  (o) 

w 


-4(0) 


"  "  ?o 


V  (0) 


V-(0) 


+  iso 


(2-57) 


The  main-diagonal  entries  of  Yq  are  the  sum  of  all  the  admittances  connected 
to  the  corresponding  wire  while  each  off-diagonal  entry  is  the  negative 
of  the  admittance  connected  between  the  appropriate  wires.  This  is  the  usual 
nodal  admittance  matrix  of  lumped-circuit  theory.  The  entries  in  1^  are  the 
values  of  the  appropriate  current  sources  connected  to  each  wire.  The  con¬ 


straint  at 


x  =  is  similar: 


V~(*  ) 
— w 


(2-58) 


Substituting  (2-57)  and  (2-58)  into  (2-55)  we  obtain 


<*i2!o'*n>  1  i  I  '*13  I  0 

I  '  i 


^22^0  *21> 


^42^0  -41^ 


^32*0  *31* 


V0) 

V.(°) 


vw(Jf) 

Va(*) 


Is(°) 


1JD 

Is<*>. 


*12  -SO 


is/  +*22i-S0 


ill  T 

~  42  -SO 


\b  I 
~  32  -so 


(2-59) 


Equations  (2-59)  may  be  reduced  by  noting  that  the  last  set  of  equations 


may  be  separated  from  the  first  three. 
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Thus 


<-^2lYo~^21') 


^42!o~!41) 


*  I 


0  1 


vs<*) 

Is^). 


^32^0  !^31) 


v°> 

^(°) 

13 

~~23 

v*> 

1 

CO 

O' 

T 1 

V0) 

_  Is<0)_ 

i 

r ' 

T 

V0) 

V.(°) 


+  <p 
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vs(0) 

is(0). 


(?12  -SO 


-Sj£  +^22  -SO 


^42  —SO 


+  ^32  -SO 


(2-60a) 


(2-60b) 


These  may  again  be  simplified.  The  first  set  of  equations  in  Equation 
(2-60a)  may  be  solved  for 


V  (*) 
w 

l 

✓—s 

o 

_ 1 

1 - 

o 

'w' 

»p 

i _ 

V~  (t  ) 

(l,12lo''?ll) 

V.(0) 

+  ?13 

Iwl 

■w 

—w 

L  ~s  -» 

+  ^12  —SO 


(2-61) 


Substituting  into  the  last  two  sets  of  equations  in  (2-60a)  yields  the  only 
simultaneous  set  of  equations  which  need  be  solved: 


U+S 


I 

( 


(^22lo"^21  YX  ~12lo  +  I*  hl)  |  ^  ^13  !^23) 

_ , _ I _ 

I 

(*42?o-?4i>  !  -♦ 


:43 


V  (0) 

— W 

V-(0) 


vs(°) 

1.(0) 


-S  Jt  +  (?22 

-?**i2>  —so 


?42  -SO 


(2-62) 


Once  the  U  +  2S  equations  in  (2-62)  are  solved  for  the  terminal  voltages 

at  x  =  0,  V  (0)  and  V-(0),  the  terminal  voltages  at  x  =  X  ,  V  (; £)  and 
*  — w  — W  w 


V^(;£),  are  obtained  directly  from  (2-61), 


42 


III.  Description  of  The  Program 


The  contents  and  operation  of  the  code  will  be  described  in  this  chapter. 
The  cards  (lines)  in  the  program  deck  are  sequentially  numbered  in  columns 
73-80  with  the  word  SHLD  in  73-76  and  the  card  number  in  77-80.  The  program 
is  written  in  Fortran  IV  language  and  is  double  precision.  A  listing  of  the 
program  is  contained  in  Appendix  A.  Steps  required  to  convert  the  program  to 
single  precision  are  given  in  Appendix  B  and  a  flow  chart  is  given  in  Appendix 
C.  In  this  flow  chart,  the  numbers  on  the  left  and  right  of  the  individual 
boxes  denote  the  beginning  and  ending  card  numbers  of  the  corresponding  portion 
of  the  code  listing,  respectively.  The  program  was  implemented  on  an  IBM 
370/165  digital  computer  at  the  University  of  Kentucky  using  the  WATFIV 
compiler,  and  the  single  precision  version  should  be  implementable  on 
other  computers.  The  program  requires  certain  function  subprograms  and  sub¬ 
routines  which  will  be  described  in  this  chapter  and  are  supplied  with  the 
main  program.  In  addition,  the  program  requires  two  subroutines  from  the 
IMSL  (International  Mathematical  and  Statistical  Library)  package  [26],  The 
first  of  these  subroutines,  LEQT1C,  solves  a  set  of  simultaneous,  complex 
equations  and  the  other,  EIGCC,  determines  the  eigenvalues  and  eigenvectors 
of  a  general,  complex  matrix.  Any  other  general  purpose  subroutines  of  this 
type  may  be  used.  See  [27]  for  a  description  of  LEQT1C  and  EIGCC  and  their 
argument  lists. 

3 . 1  Main  Program  Description 

Cards  1-37  contain  comments  concerning  the  general  applicability  of  the 
program.  Cards  38-113  contain  urray  dimensions  and  type  declarations.  The 
total  number  of  unshielded  wires  is  NU  and  the  total  number  of  shielded  cables 
is  NS.  Certain  combinations  of  these  variables  which  are  used  in  examples  of 
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dimensioning  are  given  in  cards  48-53.  Either  NU  or  NS  may  be  chosen  by  the 
user  to  be  zero.  However,  if  either  is  chosen  to  be  zero,  certain  dimensions 
must  be  set  to  1  and  not  zero.  See  cards  85  and  97. 

Cards  114-116  define  two  common  blocks  used  to  transfer  certain  constants. 
Cards  117-135  define  certain  constants.  Most  of  these  are  obvious.  Constant 
CMTM  =  2.54D-5  when  multiplied  by  a  dimension  given  in  mils  (.001  inch) 
converts  that  dimension  to  meters.  The  quantity  RADEG  =  ONE80/PI  converts 
radians  to  degrees,  and  V  =  2.997925D8  is  the  velocity  of  light  in  free  space. 
ERTE  is  the  permittivity  of  free  space. 

Cards  136- 342  read  the  input  data  concerning  dimensions  and  properties 
of  the  line  conductors  as  well  as  lengths  of  pigtail  sections  and  properties 
of  the  shields.  The  reader  is  referred  to  the  next  chapter  -  user's  manual  - 
for  a  definition  of  these  input  parameters.  Cards  344-517  read  the  terminal 
admittance,  current  source  and  shield  grounding  data. 

Cards  519-561  read  the  frequency  (F),  the  reference  conductor  impedance 
(ZG)  and  compute  the  self  impedances  of  the  unshielded  wires  (array  ZV.W),  the 
self  impedances  of  the  pigtail  wires  on  the  left  end  (array  ZWPLV)  and  right 
end  (array  ZWPRV)  of  the  line,  the  self  impedances  of  the  shields  (array  ZSV), 
the  diffusion  impedances  of  the  shields  (array  ZDV),  the  transfer  inductances 
of  the  shields  (array  LTV)  and  the  transfer  elastances  of  the  shields  (array 
STV) . 

Cards  563-633  compute  the  overall  chain  parameter  matrix  of  the  pigtail 
sections  on  the  right  end  of  the  line.  Card  570  stores  the  lengths  of  the 
right  pigtail  sections  (array  LPR)  into  array  LPT.  Cards  573-576  initialize 
the  overall  chain  parameter  matrix  of  the  right  end  to  the  identity  matrix. 
Cards  578-588  determine  the  length  of  the  first  (rightmost)  pigtail  section 


(variable  LMTNR)  as  well  as  the  longest  pigtail  section,  LMAXR,  and  initialize 

the  entries  in  array  KEY  to  indicate  whether  a  shielded  cable  is  shielded  over 

this  section  (KEY=1)  or  has  a  pigtail  over  this  section  (KEY=2).  Array  KEY 

is  the  key  to  keeping  this  information  sorted  out.  Card  589  calls  subroutine 

INDUCT  which  computes  the  per-unit-length  inductance  matrix,  L.  Card  593 

computes  Z  =  jwL  +  zqU  as  in  (2-17).  Card  596  adds  the  self  impedances  of 

the  unshielded  wires.  (See  equation  (2-11).  Matrix  Z  is  added  to  the  above.) 

~w 

Subroutine  IMP ADD  is  called  in  card  598  to  add  the  remaining  self  and  diffusion 
impedances,  Zg  and  Z^,  the  self  impedances  of  the  shielded  wires,  Z^,  as  well 
as  the  transfer  inductances  LT  to  complete  the  formation  of  (2-11)  and  thus  Z. 
Card  599  calls  subroutine  SCAP  which  forms  the  per-unit-length  elastance 
matrix  of  the  unshielded  wire,  shield  system;  that  is,  equation  (2-42)  without 
ST  and  Subroutine  ADMADD  in  card  600  adds  ST  and  S  to  this  result 

to  complete  the  formation  of  S  in  (2-42).  Cards  601-605  compute  S  ^  (and  stores 
S  1  in  array  TI).  Cards  607-610  compute  Y  =  jwC  =  jwS  ^  and  Y  \  Subroutine 
MULTC  forms  the  product  YZ  in  card  611  and  stores  in  array  TI,  and  subroutine 

EIGCC  computes  the  eigenvectors  of  YZ  (stores  in  array  T)  and  eigenvalues, 

2  -1 
Y  ,  (stores  in  array  GAM).  Cards  613-618  compute  T  and  cards  619-620 

compute  y  and  stores  in  array  GAM.  In  card  621,  subroutine  PHI  computes  the 

chain  parameter  matrix  of  this  section  (stored  in  array  PHIT)  and  MULTC  forms 

the  product  of  PHIT  and  the  previously  stored  (accumulated)  chain  parameter 

matrix  of  the  right  pigtail  section,  PHIR,  and  stores  in  array  TPHI.  Cards 

623-625  store  TPHI  in  array  PHIR.  Array  PHIR  thus  is  a  "running"  accumulation 

of  products  of  the  chain  parameter  matrices  of  these  individual  sections  of 

the  right  pigtails.  Cards  626-627  subtract  the  length  of  the  present  pigtail 

section,  LMINR,  from  each  of  the  right  pigtail  section  lengths  and  redefines 


them  in  array  LPT.  The  code  returns  to  statement  117  (card  577)  and  continues 
again  through  this  section  unless  there  are  no  nonzero  pigtail  section  lengths 
for  this  pass.  Card  581  checks  for  this  condition  and  when  no  nonzero  pigtail 
section  lengths  remain,  the  code  goes  to  statement  132  (card  630)  and  begins 
processing  the  left  pigtail  sections  in  an  identical  manner.  At  this  point 
the  overall  chain  parameter  matrix  of  the  entire  right  pigtail  section  is 
stored  in  PHIR. 

Cards  634-704  process  the  left  pigtail  sections  in  an  identical  fashion 
as  for  the  right  pigtail  sections  and  stores  the  overall  chain  parameter 
matrix  in  array  PHIL. 

Cards  705-747  compute  the  chain  parameter  matrix  for  the  shielded  section 
(a  section  with  no  pigtails  in  the  interior  of  the  line)  and  stores  it  in 
TPHI.  Card  748  computes  the  product  TPHI  *  PHIL  and  stores  in  PHIT  and  card 
749  then  computes  PHIR  *  PHIT  =  PHIR*TPHI*PHIL  and  stores  in  TPHI.  At  this 
stage,  TPHI  in  card  749  contains  the  overall  chain  parameter  matrix  of  the 
entire  line  and  it  remains  to  incorporate  the  terminal  conditions. 

If  the  line  should  contain  only  unshielded  wires,  the  code  bypasses  the 
above  pigtail  and  shield  calculations  and  computes  the  chain  parameter  matrix 
directly  in  cards  751-786  and  stores  in  TPHI. 

Cards  788-859  rearrange  the  overall  chain  parameter  matrix  rows  and 
columns  as  described  in  section  2.5.  Once  this  is  completed,  the  terminal 
conditions  are  incorporated  into  this  rearranged  chain  parameter  matrix  in 
cards  860-909,  and  equation  (2-62)  is  formed.  This  is  solved  in  card  910  and 
equations  (2-61)  are  formed  in  cards  911-924.  Array  EP  contains 


from  equation  (2-62)  and  array  EN  contains 


V  <*) 

£,<*> 


from  equation  (2-61). 

Cards  925-962  compute  the  magnitudes  and  angles  of  these  terminal 
voltages  of  the  unshielded  and  shielded  wires  and  print  the  results. 

3 . 2  Function  Subprogram  LSI 

Cards  FLS10001-FLS10014  contain  a  function  subprogram  for  computing  the 
per-unit-length  self  inductances  of  wires  above  a  ground  plane.  For  a  wire  of 
radius  r^.  at  a  height  of  above  ground,  the  result  is  [9,24] 

p  2y . 

LSI  =  an  (-— ~)  H/m 
2rr  r  . 

wi 


3 •  ^  Function  Subprogram  LM1 

Cards  FLM10001-FLM10017  contain  a  function  subprogram  for  computing  the 
per-unit-length  mutual  inductances  between  two  wires  above  a  ground  plane.  For 
two  wires  at  heights  y^  and  y  and  horizontal  coordinates  z^  and  z^ ,  the  result 
is  [1,9,24] 


LM1 


4yiyi 

+  V 


H/m 


where 


d.  . 
U 


+ 


3.4  Function  Subprogram  LS2 

Cards  FLS20001-FLS20014  contain  a  function  subprogram  for  computing  the 
per-unit-length  self  inductances  of  wires  within  an  overall,  circular, 


47 


cylindrical  shield.  For  an  overall  shield  of  interior  radius  R  and  a  wire 
of  radius  r^.  located  r^  from  the  axis  of  the  shield,  the  result  is  [9,24] 


LS2  in 

2n 


H/m 


3.5  Function  Subprogram  LM2 

Cards  FLM20001-FLM20019  contain  a  function  subprogram  for  computing  the 
per-unit-length  mutual  inductances  between  wires  within  an  overall,  circular, 
cylindrical  shield.  For  an  overall  shield  with  interior  radius  R  and  wires 
which  are  at  radii  r^  and  r  from  the  axis  of  the  shield  and  separated  by 
angle  0„,  the  result  is  [9,24] 


LM2  = 


2u 


Zn  < 


r . 

(-1) 
v  R  ' 


2  4  2  "I  1 

(r.r.)  +  R  -  2r.r.  R  cos(0..) 

JUU _ _ _ _ . .  .  .  >J. 


(r.r.)^  +  r J*  -  2r.r.^  cos(0..) 
i  3  J  i  J  ij 


H/m 


3.6  Function  Subprogram  STB 

Cards  FSTB0001-FSTB0073  contain  a  function  subprogram  for  computing  the 
per-unit-length  transfer  elastances  for  braided  shields  which  are  entires  in 
S^,  in  equation  (2-42).  A  table  of  values  of  this  parameter  for  discrete 
values  of  braid  weave  angle  for  angles  less  than  or  equal  to  45°  is  given  in 
[22].  An  interpolation  routine  described  in  [23]  is  used  to  interpolate  this 
table. 

3. 7  Function  Subprogram  LTB 

Cards  FLTB0001-FLTB0073  is  a  function  subprogram  for  computing  the  per- 
unit-length  transfer  inductances  for  braided  shields  which  are  entries  in  L 
in  equation  (2-11).  A  table  of  values  of  this  parameter  for  discrete  values 
of  braid  weave  angle  for  angles  less  than  or  equal  to  45°  is  given  in  [22]. 

An  interpolation  routine  described  in  [23]  is  used  to  interpolate  this  table. 
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IT— TT- 

Lc 

-- 

•  ~ 

■ 

3.8  Function  Subprogram  ZWW 

Cards  FZWW0001-FZWW0027  contain  a  function 

subprogram  for  computing  the 

-1 

per-unit-length  self  impedances  of  stranded  wires.  The  per-unit-length  self 

. 

J 

impedance  of  one  strand  is  determined,  zgT»  and 

the  NST  strands  are  treated 

s 

as  being  connected  in  parallel  so  that  the  total  per-unit-length  self 

impedance  of  the  stranded  wire  is  ZWW  =  z^/NST 

ill m.  The  programmed  equations 

are  described  in  [1,27].  For  a  solid  cylinder  of  radius  r  and  conductivity 

w 

1 

a,  the  skin  depth  is 

< 

-■ 

6  =  1 

/rrfu  a 

V 

N 

i 

2tt  /af  x  10  7 

'#  1 

■  ■} 

the  D-C  per-unit-length  resistance  is 

.  ’  * 

r 

rDC  '  2  a/" 

ire  r 

w 

■- 

-•  ■ 
•  J 

" 

and  the  D-C  per-unit-length  internal  inductance 

is 

5  Uv 

DC  8  it 

V  * 

■  'A 

=  .5  x  10-7  H/m 

Includirg  skin  effect,  the  programmed  equations 

are  [1,27] 

• 

-  A 

(!)  rw  <  6 

t  -  r0(.  S/m 

*  ■  «DC  H/m 

\ 

• 

-  _  tl 

•  ' 

-  ■  A 
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" 

... 

•  _ 

a.-*—,  a— a— 

(II) 


6  <  r  <36 

w 


r  =  i  +  3>roc  u/m 


a  =  [1.15  -  .15  (y)]  £dc  H/m 


(III) 


r  >  36 

w  — 


r  ■  2?  rDC  a/m 


1  ‘  f  lDC  H/” 
w 


The  per-unit-length  self  impedance  of  each  strand  of  radius  is 


zgT  =  r  +  j«£  n/m 


3 . 9  Function  Subprogram  2DB 

Cards  FZDB0001-FZDB0023  contain  a  function  subprogram  for  computing  the 
per-unit-length  diffusion  impedances  of  braided  shields.  For  a  braided  shield 


having  B  belts,  WPB  wires  per  belt  and  weave  angle  0  ,  the  result  is  [1] 


y  2ru 


ZDB  =  r 


DC  sinh(2Yr,  ) 

D 


n/ 


m 


where  the  braid  wires  have  radius  y. 


y  -  (1  +  jl)/6 


6  = 


/irfp  a 

v 


and  r  is  the  per-unit-length  D-C  resistance  of  the  braid  [1] 


nr  2 

tt  or*  B  WPB  co s  (G  ) 

D  W 


n/ 


m 


(All  braid  wires  are  treated  as  being  connected  in  parallel.) 
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3.10  Function  Subprogram  ZSB 


Cards  FZSB0001-FZSB0029  contain  a  function  subprogram  for  computing  the 
per-unit-length  self  impedance  of  a  braided  shield.  The  per-unit-length 
impedance  of  each  isolated  braid  wire,  z^,  is  computed  as  described  in 
section  3.7  and  the  result  is  [1] 


ZSB  = 


B  WPB  Cos (0  ) 
w 


ft/m 


3.11  Function  Subprogram  ZDS 

Cards  FZDS0001-FZDS0020  contain  a  function  subprogram  for  computing  the 
per-unit-length  diffusion  impedances  of  solid  shields.  The  shield  is  treated 
as  a  thin-walled  (relative  to  the  shield  radius)  cylinder  and  the  equations 
are  [17] 

Y  t 


ZDS  =  r 


DC  sinh  (yt  ) 
'  s 


ft/m 


where  the  shield  is  of  conductivity  a,  interior  radius  r  ,  thickness  t  and 

s  s 


DC  rr  a  t  (2r  +  t  ) 

s  s  s 


ft/m 


5  = 


V 


y  =  (l+jl)/6 


3.12  Function  Subprogram  ZSS 

Cards  FZSS0001-FZSS0024  contain  a  function  subprogram  for  computing  the 
per-unit-length  self  impedances  of  solid  shields.  The  shield  is,  as  in  ZDS, 
treated  as  a  thin-walled  cylinder  and  the  equations  are  described  in  [27]: 

ZSS  =  r  +  jcoJ.  ft/m 


51 


and  6  is  the  skin  depth. 

3.13  Subroutine  MULTC 

Cards  SMUL0001-SMUL0017  contain  a  subroutine  for  multiplying  two  complex 
matrices  B  and  C  as 

A  =  B  C 

A  is  of  dimension  NL  x  NR, 

B  is  of  dimension  NL  x  NM,  and 

C  is  of  dimension  NM  x  NR. 

3 . 14  Subroutine  SCAP 

Cards  SCAP0001-SCAP0024  contain  a  subroutine  for  computing  a  portion  of 
the  per-unit-length  elastance  matrix,  S,  described  in  section  2.1.  The 
routine  computes  equation  (2-42)  for  only  the  system  of  U  unshielded  wires  and 
S  shields,  that  is,  S  without  submatrices  ST  and  SIND  Those  submatrices  are 

added  to  this  S  in  subroutine  ADMADD.  The  matrix  to  be  computed  is 


Here  w  denotes  unshielded  wires,  s  denotes  shields,  p  denotes  pigtail  wires 
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and  w  denotes  shielded  wires#  First  the  per-unit-length  elastance  matrix  is 
computed  from  the  per-unit-length  inductance  matrix  L  as 


Then  a  check  is  made  using  array  KEY  to  determine  if,  in  fact,  a  shield  is 
actually  present  over  this  section  rather  than  a  pigtail  wire.  If  so,  those 
entries  in  S~~  need  to  be  changed  to  the  corresponding  ones  in  S  as 

WW  ^  S  S  £  p  p 

indicated  in  equation  (2-42).  Entries  in  S  A  and  S  A  ~  are  already  computed 

^ ww  ^ sw ^ pw 

correctly  in  L  by  subroutine  INDUCT. 

3.15  Subroutine  INDUCT 

Cards  SIND0001-SIND0169  contain  a  subroutine  for  computing  the  per-unit- 
length  inductance  matrix  L  given  in  equation  (2-14)  -  (2-17): 


L 

L 

L  - 

.ww 

~ws,wp 

_ww 

L  1 

L 

L  -  - 

_ws ,wp 

~ss,pp 

_sw,pw 

L  * 

L  - 

„ww 

_ sw , pw 

.ww 

subscripts  w,  p,  s,  w  have  the  same  meaning  here  as  in  the  previous  section. 

If  over  this  section,  no  shields  are  involved,  i.e.,  each  shield  is  in  fact  a 

pigtail  wire,  then  L  is  computed  in  a  straightforward  fashion  as  that  of  a 

system  of  U  unshielded  wires,  S  shielded  wires  and  S  pigtail  wires.  If,  over 

this  section,  a  shield  is  present  (determined  by  array  KEY),  then  appropriate 

entries  in  L  -  are  equal  to  the  corresponding  entries  in  L  since  the 

_ww  1  „ws,wp 

mutual  inductance  between  an  unshielded  wire  and  a  shielded  wire  is  the  same  as 
the  mutual  inductance  between  the  unshielded  wire  and  that  shield  [1]. 
Similarly,  certain  entries  in  ^  are  equal  to  the  corresponding  entries  in 
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Ci 


L  since  the  mutual  inductance  between  a  shield  and  its  shielded  wire  is 

~ss,pp 

the  same  as  the  self  inductance  of  that  shield  [1], 

3.16  Subroutine  PHI 

Cards  SPHI0001-SPHI0054  contain  a  subroutine  for  computing  the  chain 
parameter  matrix  of  a  section  of  line  as  described  in  section  2.3.  To  save 
array  space  several  dummy  arrays  are  used  temporarily  and  written  over. 

3.17  Subroutine  ADMADD 

Cards  SADM0001-SADM0030  contain  a  subroutine  for  adding  S„  and  STX_  to 

..T  _ IND 

the  S  matrix  generated  by  subroutine  SCAP  to  yield  the  per-unit-length 
elastance  matrix  S  in  equation  (2-42). 

3. 18  Subroutine  IMPADD 

Cards  SIMP0001-SIMP0039  contain  a  subroutine  to  add  the  self  impedances, 
diffusion  impedances  and  transfer  inductances  (Zg,  Z^,  L^,)  to  the  per-unit- 
length  impedance  matrix  to  yield  Z  in  equation  (2-11).  The  self  impedances 
of  the  shielded  wires,  Z- ,  are  also  added 

'  T.T  ' 
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IV.  User's  Manual 


In  this  section  we  will  provide  a  user's  manual  showing  the  format  and 
description  of  the  required  input  data.  The  program  requires  that  the  input 
data  groups  be  organized  as  shown  in  Table  1  where  the  number  of  unshielded 
wires  is  U  and  the  number  of  shielded  cables  is  S.  These  card  groups  must 
follow  the  main  program,  function  subprograms  and  subroutines  in  the  order 
shown  in  Table  1.  The  data  entries  on  the  cards  are  either  in  Integer  (I) 
format,  e.g.,  35,  or  in  Exponential  (E)  format,  e.g.,  12.6E-3.  All  data 
entries  must  be  right- j u^tif ied  in  the  assigned  card  column  block.  These  data 
entries  are  printed  out  by  the  program.  It  is  highly  recommended  that  the 
user  check  this  printout  of  input  data  to  insure  that  the  input  data  are  as 
intended. 

In  the  main  program,  the  user  must  appropriately  dimension  all  arrays  for 
each  problem.  Comment  cards  are  provided  at  the  beginning  of  the  main  program 
to  assist  the  user  in  providing  proper  dimensions. 

Program  SHIELD  considers  a  transmission  line  consisting  of  U  unshielded 
wires  and  S  shielded  cables  and  an  associated  reference  conductor.  The 
reference  conductor  type  is  determined  by: 

TYPE  =  1:  The  reference  conductor  is  a  ground  plane. 

TYPE  =  2:  The  reference  conductor  is  an  overall,  circular, 
cylindrical  shield. 

Cross-sectional  views  of  these  two  structures  are  shown  in  Fig.  4-1  and  Fig. 
4-2. 

4 . 1  Transmission  Line  Characterization  Cards,  Group  I 

The  purpose  of  this  group  is  to  define  the  type  of  reference  conductor, 
the  total  line  length,  the  number  of  unshielded  wires  and  the  number  of 
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TABLE  1 


3 


Card  Group  Organization  Data 

(Note:  the  card  groups  must  be  arranged  in  the  following  order.; 


L- 


N 


1)  Line  Characterization 

2)  Unshielded  Wires  (Omit  if  no  unshielded  wires) 

Wire  1 
Wire  2 

i 

i 

I 

Wire  U 

3)  Shielded  Cables  (Omit  if  no  shielded  cables) 

Cable  1 


Cable  2 


Cable  S 


Card  Group 
I 


II 

II 

i 

i 

i 

II 


III(c) 


Table 

2 

3 

3 

i 

l 

l 

3 

4 

5 

6 

4 

5 

6 

4 

5 

6 
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Card  Group 


Table 


m 

J.; 

i 


TABLE  1  (Cont  *d) 

4)  Terminal  Constraints 

(a)  Unshielded  Wires  (Omit  if  no  unshielded  wires) 
Wire  1 
Wire  2 


Wire  U 

(b)  Shielded  Cables  (Omit  if  no  shielded  cables) 

Cable  1 

Cable  2 

i 

i 

I 

Cable  S 

(c)  Mutual  Admittances 

(See  Table  9  for  order  information) 

5 )  Frequency  and  Reference  Conductor  Impedance 
frequency  1 
frequency  2 
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IV(a) 
IV  (a) 

i 

i 


IV  (a) 


IV(b) 

IV(b) 


IV(b) 


IV  (c) 


V 

V 


10 

10 


6* 


Reference  Conductor 
(Ground  Plane) 


shielded  wires.  In  addition,  if  the  reference  conductor  is  an  overall, 

circular,  cylindrical  shield  (TYPE  2)  as  shown  in  Fig.  4-2,  then  the  interior 

radius  of  this  shield,  R  is  defined. 

s 

The  format  of  these  cards  and  their  order  is  shown  in  Table  2. 

4 . 2  Unshielded  Wire  Characteristics  Cards,  Group  II 

The  purpose  of  this  card  group  is  to  define  the  relative  positions  and 
the  characteristics  of  the  U  unshielded  wires.  The  format  and  order  of  these 
cards  is  shown  in  Table  3. 

4 . 3  Shielded  Cable  Characteristics  Cards,  Group  III  (a) 

The  purpose  of  these  cards  is  to  define  the  relative  positions  and  the 
characteristics  of  the  shielded  cable  shields  and  the  characteristics  of  their 
shielded  wires.  The  format  and  ordering  of  these  cards  is  given  in  Table  4. 

4 . 4  Pigtail  Wire  Characteristics  Cards,  Groups  111(b)  and  III(c) 

The  purpose  of  these  groups  is  to  define  the  relative  locations  and  the 
characteristics  of  the  pigtail  wires.  The  cards  for  the  pigtails  at  the  left 
(x  =  0)  end  are  Group  111(b),  and  their  format  and  ordering  are  given  in 
Table  5.  The  similar  cards  for  the  pigtails  at  the  right  (x  =  ^  )  end  are 
Group  III(c)  and  are  described  in  Table  6. 

The  angular  positions  of  the  pigtail  wires  are  measured  in  the  counter¬ 
clockwise  direction  when  looking  from  x=0  towards  x=  . 

4 . 5  The  Terminal  Characteristics  Cards,  Group  IV 

The  general  structure  of  the  terminal  networks  is  illustrated  in  Fig.  2-5. 
At  x=0  and  x=^T  an  admittance  in  parallel  with  a  current  source  is  connected 
between  each  unshielded  wire  and  the  reference  conductor  and  between  each 
shielded  wire  and  the  reference  conductor.  Card  Group  IV (a)  (Table  7)  defines 
these  for  the  unshielded  wires,  and  Card  Group  IV (b)  (Table  8)  defines  these 
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TABLE  2 


Format  of  the  Line  Characterization  Cards,  Group  I 


Card  Column 


Card  #1 

TYPE  (1,2)  1 

Card  If 2 

£  (total  line  length  in  meters)  1-10 


Card  #3 

U  (number  of  unshielded  wires)  1-2 

Card  #4 

S  (number  of  shielded  cables)  1-2 


Format 


I 


E 


I 


I 


Card  it 5  (Omit  if  TYPE*!) 


TABLE  3 


Format  of  the  Unshielded  Wire  Characteristics  Cards,  Group  II 


(Omit  if  no  unshielded  wires) 

Card  Column 

Card 

*1 

(a) 

Wi 

(wire  number) 

1-2 

TYPE= 

-1: 

(b) 

Y  . 

Wl 

(height  of  wire  above  ground  in  meters) 

11-20 

(c) 

Z  . 

Wl 

(horizontal  coordinate  of  wire  in  meters) 

31-40 

TYPE= 

■2: 

(b) 

R  . 

Wl 

(radial  distance  of  wire  from  axis  of 

11-20 

overall  shield  in  meters) 

(c) 

0  . 
Wl 

(angular  position  of  wires  in  decrees) 

31-40 

Card 

J2 

(a) 

r  . 

Wl 

(radius  of  wire  in  mils) 

11-20 

(b) 

(number  of  strands) 

34-35 

(c) 

(radius  of  strands  in  mils) 

51-60 

(d) 

(conductivity  of  strands  relative  to  copper) 

71-80 

Format 

I 

E 

E 

E 


E 

I 

E 

E 
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TABLE  4 


Format  of  the  Shielded  Cable  Characteristics  Cards,  Group  III (a) 
(Omit  if  no  shielded  cables) 


Card  #1 

(a)  (shielded  cable  number) 

(b)  (shield  type:  l(solid),  2 (braided)) 
TYPE=1: 


Card  Column  Format 

1-2  I 

10  I 


(c)  (height  of  shield  above  ground  in  meters)  21-30  E 

(d)  (horizontal  coordinate  of  shield  in  meters)  41-50  E 

TYPE=2: 

(c)  R  .  (radial  distance  of  shield  from  axis  of  21-30  E 

si 

overall  shield  in  meters) 

(d)  0gi  (angular  position  of  shield  in  degrees)  41-50  E 


Card  #2 

(a)  r  .  (interior  radius  of  shield  in  meters)  11-20 

si 

(b)  (relative  permittivity  of  dielectric)  21-30 

(c)  r~ .  (radius  of  shielded  wire  in  mils)  31-40 

wi  — — 

(d)  (number  of  strands  in  shielded  wire)  55-56 

(e)  (radius  of  strands  in  shielded  wire  in  mils)  61-70 

(f)  (conductivity  of  shielded  wire  strands  71-80 

relative  to  copper) 


E 

E 

E 

I 

E 

E 


Card  #3 

Solid  Shields: 

(a)  tg^  (shield  thickness  in  mils)  11-20  E 
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TABLE  4  (Continued) 


Card  Column 


Format 


Card  // 3  (Cont'd) 


(b)  (conductivity  of 

Braided  Shields: 

shield  relative  to  copper) 

21-30 

E 

(a) 

(radius  of  braid 

wires  in  mils) 

11-20 

E 

(b) 

(conductivity  of 

braid  wires  relative  to  copper) 

21-30 

E 

(c) 

(weave  angle  of  braid  in  degrees) 

31-40 

E 

(d) 

(number  of  belts 

in  braid) 

44-45 

I 

(e) 

(number  of  wires 

per  belt) 

54-55 

I 

TABLE  5 


TABLE  6 


Format  of  Pigtail  Wire  Characteristics  Cards  for 


1ABLE  7 


ira 


1 


Format  of  Terminal  Characteristics  Cards,  Group  IV 
Self  Terms  for  Unshielded  Wires,  Group  IV (a) 
(Omit  if  no  unshielded  wires) 


Card  Column 


Card  #1 


(a) 


SOi  (current  source  between 
unshielded  wire  and 
reference  conductor  at 
x=0  in  amps) 

(b)  (admittance  between 

unshielded  wire  and 
reference  conductor  at 
x=0  in  Siemens) 
(current  source 
between  unshielded 
wire  and  reference 
conductor  at  x=jf  in 
amps) 

(admittance  between 
unshielded  wire  and 
reference  conductor 
at  x=  X  in  Siemens) 


real  part 


imaginary  part 


(■ 


real  part 


(c) 


f 


real  part 


(d)  Y 


1-10 


11-20 


21-30 


imaginary  part  31-40 


41-50 


imaginary  part  51-60 


real  part  61-70 

imaginary  part  71-80 
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Format 

E 

E 

E 

E 

E 

E 

E 

E 


TABLE  8 


Format  of  Terminal  Characteristics  Cards,  Group  IV 
Self  Terms  for  Shielded  Wires,  Group  IV (b) 


(Omit  if  no  shielded  cables) 


Card  Column  Format 


Card  //I 


Igo i  (current  source  between 
shielded  wire  and 


reference  conductor  at 


x=0  in  amps) 

Yrt  . .  (admittance  between 
Oil 

shielded  wire  and 


reference  conductor 


real  part 


imaginary  part  11-20 


real  part 


21-30 


at  x=0  in  Siemens) 


(current  source  be¬ 


tween  shielded  wire 


and  reference  conduc- 


imaginary  part  31-40 


real  part 


41-30 


tor  at  x=  X  in  amps ) 
Y  - . .  (admittance  between 

i  11 

shielded  wire  and 


imaginary  part  51-60 


real  part 


61-70 


reference  conductor 


at  x=  X  in  Siemens) 


imaginary  part  71-80 


Card  #2 


IGNDL  (grounding  code  for  left  end:  1  =  short  10 

2  =  open  ) 

IGNDR  (grounding  code  for  right  end:  1  =  short  40 


TABLE  9 


Ti 


Format  of  Terminal  Characteristics  Cards.  Group  IV 


•v 

Mutual  Admittances  Group  IV(c) 

Card  Column 

Format 

3 

Card  #1 

r* 

l'. 

(a) 

Y  .  (admittance  between 

r 

IL 

OiJ 

21-30 

r 

conductor  i  and  con¬ 

real  part 

E 

S* 

S 

ductor  j  at  x=0) 

imaginary  part 

31-40 

E 

(b) 

Y  -  (admittance  between 

i  ^ 

f 

r- 

real  part 

61-70 

E 

conductor  i  and  con¬ 

t* 

* 

ductor  j  at  x*J) 

imaginary  part 

71-80 

E 

Note:  Cards  must  be  arranged  in  order  (w  denotes  unshieded  wire  and 


s  denotes  shielded  wire)  so  as  to  fill  the  portion  of  the 
following  matrix  above  the  main  diagonal  b rows: 


wl  w2  w3  - wU  si  s2  - sS 
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1 


TABLE  9  (Continued) 

Thus  the  cards  are  arranged  sequentially  in  the  order  of  the  admittance 
connected  between 


wl  — ► 

w2, 

w3,  — 

- ,  wU ,  si. 

s2, 

- ,  sS 

w2  — ► 

w3. 

w4,  — 

- ,  wU,  si, 

s2, 

- ,  sS 

i 

i 

i 

i 


i 


w(U-l)  *• 

wU, 

si 

wU  - •> 

si. 

s2 

si  — 

s2, 

s3 

i 

i 

i 

s(S-2)  — *► 

s(S- 

-1) 

s(S-l)  — » 

sS 

s2, - ,  sS 

- ,  sS 

- ,  sS 


sS 


A  total  of 


(U+S)  (U+S-l) 
2 


cards  are  in  this  group 


TABLE  10 


Format  of  the  Frequency  and  Reference  Conductor  Impedance 
(Per-Unit-Length)  Cards,  Group  V 

Card  column 

Card  #1 

F  (frequency  in  Hz)  1-10 

Card  #2 

Zq  (reference  conductor  impedance  J  real  part  1-10 

per  unit  length  in  ohms/meter)  (  imaginary  part  11-20 

Note:  Group  V  may  be  repeated  at  the  end  of  the  input  deck  in  order 

for  the  program  to  analyze  the  structure  for  more  than  one 
frequency  without  having  to  redefine  the  previous  structure  and 
terminal  data  for  each  new  frequency.  However  the  previous 
data  Groups  I,  II,  III,  IV  are  assumed  to  apply  to  all  these 
frequencies. 


Format 

E 

E 

E 


for  the  shielded  wires.  The  current  sources  at  x=0  and  x=jL  are  directed 


from  the  reference  conductor  _to  the  end  of  the  wire.  In  addition,  the 
grounding  of  the  shields  at  x=0  and  x= X  are  defined  in  Card  Gro>  p  IV(b) 

(Table  8).  Each  shield  is  connected  via  a  pigtail  wire  (which  may  be  of 
zero  length)  to  the  reference  conductor  by  either  a  short  circuit  or  an  open 
circuit . 

A  mutual  admittance  (which  may  be  zero,  i.e.,  absent)  is  connected 
between  the  ends  of  all  pairs  of  these  wires.  The  format  and  ordering  of 
these  cards  in  Group  IV (c)  is  given  in  Table  9.  The  user  must  be  careful  to 
adhere  to  the  ordering  of  these  cards  as  described  in  Table  9.  The  ordering 
scheme  is  very  simple.  First  describe  the  admittances  connecting  unshielded 
wire  #1  and  all  other  wires  sequentially  in  the  order  unshielded  wires  first 
and  then  the  shielded  wires.  Next  describe  the  admittances  connecting  un¬ 
shielded  wire  #2  and  all  other  wires  sequentially  in  the  order  unshielded 
wires  first  and  then  the  shielded  wires.  Continue  in  this  fashion  until  the 
admittances  connecting  the  last  unshielded  wire,  wire  //U,  and  all  other  wires 
(in  this  case  only  the  shielded  wires)  have  been  described.  Then  describe  the 
admittances  connecting  the  first  shielded  wire  and  the  remaining  shielded  wires, 
sequentially.  Then  describe  the  admittances  connecting  the  second  shielded  wire 
and  the  remaining  shielded  wires,  sequentially.  Continue  until  the  admittance 
connecting  the  next  to  last  shielded  wire  (S-l)  and  the  last  shielded  wire  (S) 

Is  described  and  the  process  is  complete.  The  total  number  of  cards  in  this 
group  is 

(U+S)  (IJ+S-1) 

2 

If  the  wires  are  numbered  1,  2,  - ,  U  and  the  shielded  wires  are  re- 
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numbered  U+l ,  U+2 


U+S,  then  the  ordering  is 


1  (U+S) 


‘2  (U+S) 


(U-l)  u 


(U-l)  (U+S) 


U  (U+l) 


U  (U+S) 


(U+l)  (U+2) 


(u+l)  (U+S) 


(U+S-l)  (U+S) 


4 . 6  The  Frequency  and  Reference  Conductor  Impedance  Cards,  Group  V 

The  final  card  group  defines  (1)  the  frequency  for  the  analysis  and  (2) 
the  per-unit- length  impedance  of  the  reference  conductor  at  this  frequency. 
This  card  group  may  be  repeated  at  the  end  of  the  deck  so  that  the  program 
will  analyze  the  structure  at  more  than  one  frequency  without  the  need  to  re¬ 
peat  the  previous  line  and  termination  characterization  groups.  The  program 
will  process  the  data  in  Groups  I,  II,  III  and  IV  and  compute  the  terminal 
voltages  at  the  frequency  on  the  first  frequency  card  it  encounters.  It  will 
then  recompute  the  line  response  at  each  frequency  on  the  remaining  frequency 
cards  (and  use  the  reference  conductor  impedance  on  the  following  card).  The 
program  assumes  that  the  data  on  card  groups  1,  II,  III,  IV  are  to  be  used  for 
all  the  remaining  frequencies.  If  this  is  not  intended  by  the  user,  in 
particular  if  the  terminal  characterizations  are  frequency  dependent,  then  one 
may  only  run  the  program  for  one  frequency.  This  feature,  however,  can  be 
quite  useful.  If  the  termination  networks  are  purely  resistive,  i.e.,  fre¬ 
quency  independent  over  the  frequency  range  of  interest,  then  one  may  use  as 
many  repetitions  of  this  group  as  desired  and  the  program  will  compute  the 
terminal  responses  at  each  frequency  without  the  necessity  for  the  user  to 
input  the  data  in  Groups  I,  II,  III,  IV  for  each  frequency. 

4.7  Range  of  the  Values  of  the  Input  Data 

The  only  restrictions  on  the  values  of  the  input  data  are  that  only 
certain  input  parameters  can  be  input  as  zero  values .  These  are: 

(1)  NU  (the  number  of  unshielded  wires) 
or 

NS  (the  number  of  shielded  wires) 

(Note:  if  NU=0  or  NS=0  certain  arrays  must  be  dimensioned  of 
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size  1  not  size  0.  See  cards  SHLD0085  and  SHLD0097.) 

(2)  Y  .  or  Z  .  and  0  .  in  Table  3.  (Note:  No  two  wir*3  locations 

W1  W1  Wl 

may  be  the  same,  i.e.,  (Y  .,Z  .)  and  Y  .,Z  .)  or  (R  .,0  .)  and 

Wl’  Wl  WJ  W j  Wl*  Wl 


(R  .,0  .). 

wj  Wj 


(3)  Y  .  or  Z  .  and  0  .  in  Table  4.  (Note:  No  two  shield  locations 
si  —  si  si 


may  be  the  same,  i.e.,  (Y  .,Z  .)  and  (Y  .,Z  )  or  (R  . ,0  )  and 

si  si  sj’  sj  si’  si' 


(R  .,0  .).) 

sj  sj 


(4)  ^  (The  pigtail  lengths  for  left  or  right  or  both  pigtails. 


See  Table  5  and  Table  6). 


(5)  (The  angular  positions  of  the  pigtail  wires.  See  Tables  5 


and  Table  6. ) 

(6)  Any  of  the  terminal  current  source  or  admittance  values.  (see 
Table  7,  Table  8  and  Table  9.) 

(7)  Zq  (The  reference  conductor  impedance  at  a  frequency. 

See  Table  10.) 

All  of  the  other  data  entries  should  be  nonzero  (and  represent  realistic 


values) . 


V.  Program  Checkout 


In  order  to  check  the  proper  functioning  of  the  program,  several  examples 
were  run.  Examples  consisting  of  2  unshielded  wires  and  no  shit  Ided  wires 
were  run  for  TYPE  1  and  TYPE  2  structures.  The  outputs  of  SHIELD  were  com¬ 
pared  to  those  of  program  XTALK2  described  in  [27]  and  were  found  to  be 
identical  as  they  should  be  since  XTALK2  performs  the  same  calculations  for 
unshielded  wires  as  does  SHIELD.  SHIELD  was  also  run  for  2  shielded  wires 
each  having  8cm  pigtail  sections  on  both  ends  and  separated  by  1.5cm  and  1.5cm 
above  ground.  These  data  were  measured  experimentally  and  reported  in  [1]. 

In  addition  a  computer  program  was  written  for  [1]  in  a  "brute  force"  manner 
to  show  predictions  of  these  experimental  data.  Those  predictions  are 
reported  in  [1].  SHIELD  was  run  for  this  configuration  and  the  predictions 
compared  favorably  with  the  experimental  data  in  [1]  and  compared  almost 
exactly  with  predictions  given  in  [1]  by  the  "brute  force"  program. 

As  an  integrated  test,  we  performed  an  experiment  to  illustrate  the  pre¬ 
diction  accuracies  of  SHIELD.  The  experimental  configuration  is  shown  in 
Fig.  5-1.  Photographs  of  the  experimental  setup  are  given  in  Fig.  5-2.  The 
unshielded  wires  are  #22  gauge  stranded  (radius  =  15  mils)  with  7  strands  of 
#30  gauge  wire  (radius  =  5  mils).  The  insulation  is  pvc  (neglected  by  the 
program).  The  shielded  cables  are  identical  to  those  used  in  [1],  They  are 
braided  with  inner  shield  radius  of  35  mils  =  8.89E-4  m  and  an  inner,  shielded 
wire  which  is  #22  gauge  (radius  =  15  mils)  composed  of  7  strands  of  #30 

gauge  wire  (radius  =  5  mils).  The  inner  dielectric  is  teflon  with  =  2.1. 
The  braid  wires  have  radii  of  2.5  mils,  and  the  braid  consists  of  16  belts 
with  4  wires  per  belt  and  a  weave  angle  of  30  degrees.  The  pigtail  wires  are 
#20  gauge  (radius  =  16  mils)  solid  bare  copper  wires.  A  listing  of  the  input 
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data  is  given  in  Fig.  5-3.  A  sample  printout  of  some  of  the  results  for 
lOOkllz  — *■  900kHz  is  given  in  Fig.  5-4. 

The  unshielded  wires,  shielded  cables  and  shielded  wires  are  numbered  in 
Fig.  5-1  as  ul,  u2,  si,  s2,  s3,  wl,  w2,  w3.  In  the  experiment,  an  oscillator 
was  attached  to  the  right  end  of  the  second  unshielded  wire  (u2),  and  the 
voltage  of  the  oscillator  was  increased  until  the  input  voltage  to  that  line 
was  1  volt.  Then  the  voltages  coupled  to  the  other  ends  of  the  other  wires 
ul,  wl,  w2,  w3  were  measured  resulting  in  voltage  transfer  ratios.  A  zero 
source  impedance  source  is  not  allowed  in  the  SHIELD  input  data.  Thus  to 
simulate  this  situation,  the  source  at  the  right  end  of  the  second  unshielded 
wire  is  50V  with  source  resistance  of  50ft  or,  equivalently,  a  1A  source  in 
parallel  with  a  resistance  of  50ft  (2.0E-2  Siemens  admittance).  (See  the  list¬ 
ing  in  Fig.  5-3.)  Once  SHIELD  was  run  for  this  configuration,  the  voltage 
transfer  ratios  corresponding  to  the  experimental  data  are  obtained  by 
dividing  the  voltages  at  the  ends  of  each  line  by  the  voltage  at  the  right  end 
of  the  second  unshielded  wire. 

The  voltage  transfer  ratios  (magnitude  and  angle)  predicted  by  SHIELD  and 
experimentally  measured  at  the  right  end  (x  =Ji  )  of  each  line  are  shown  in 
Fig.  5-5  through  Fig.  5-8.  Fig.  5-5(a)  and  Fig.  5—5 (b )  show  the  magnitude  and 
angle,  respectively,  of  the  measured  voltage  at  the  right  end  of  the  first 
unshielded  wire,  V  ,  versus  the  predicted  voltage  transfer  ratio  of  SHIELD: 

vwl(*> 

Vw2(#> 

The  predictions  are  quite  good  and  well  within  ldB  below  10  MHz.  Above  10  MHz 
the  line  is  becoming  electrically  long,  and,  consequently,  rapid  variations 
in  both  the  magnitude  and  the  phase  of  the  voltage  transfer  ratio  are  observed. 
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Figure  5-3.  Input  data  for  the  experiment-line  structure  characterization  data 
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Figure  5-3  (Cont'd)...  Input  data  for  the  experiment-unshielded  wire  data. 
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Figure  5-3  (Coat'd)...  Input  data  for  the  experimental-shielded  cable  data. 
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Figure  5-3  (Cont'd)...  Input  data  for  the  experiment-terminal  characterization 
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Figure  5-3  (Cont'd)...  Input  data  for  the  experiment-frequency  and  common 

impedance  data. 
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Figure  5-4.  Printout  of  computed  results  for  the  experiment 
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Figure  5-4  (Cont’d)...  Printout  of  computed  results  for  the  experiment. 
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Figure  5-4  (Cont'd)...  Printout  of  computed  results  for  the  experiment 
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Figure  5-4  (Cont'd)...  Printout  of  computed  results  for  the  experiment 
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Figure  5-4  (Cont'd)...  Printout  of  computed  results  for  the  experiment. 
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Figure  5-4  (Cont'd)...  Printout  of  computed  results  for  the  experiment 
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Figure  5-4  (Cont'd)...  Printout  of  computed  results  for  the  experiment. 
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Figure  5-7 (a).  Magnitude  of  voltage  transfer  ratio  to  right  end  of  s2,  w2.  0 
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Nevertheless,  the  program  predicts  these  results  rather  well. 


A  similar  comparison  of  the  voltage  transfer  ratio  to  the  right  end  of 
the  first  shielded  wire,  V  and  the  predictions  of  SHIELD, 

Vi» 

V„2<*> 

are  shown  in  Fig.  5-6.  Note  that  above  1  MHz,  the  coupling  which  had  become 
"flat"  for  lower  frequencies  tends  to  rise.  This  is  a  result  of  the  coupling 
via  the  two  3cm  pigtail  sections  dominating  the  coupling  to  the  (much  longer) 
shielded  section  as  is  explained  in  [1,2].  SHIELD  predicts  these  results 
rather  well. 

The  comparison  of  the  voltage  transfer  ratio  to  the  right  end  of  the 
second  shielded  wire,  V^,  and  the  predictions  of  SHIELD, 

V02 

Vw2(  *> 

are  shown  in  Fig.  5-7.  Here  the  left  pigtail  (3cm)  is  grounded  but  the  right 
pigtail  (zero  length)  is  open.  The  predictions  of  SHIELD  are  again  quite 
accurate.  The  comparison  of  the  voltage  transfer  ratio  to  the  right  end  of 
the  third  shielded  wire,  V and  the  predictions  of  SHIELD, 

V03(*> 

V„2<*> 

are  shown  in  Fig,  5-8  and  are  quite  good.  For  this  shield,  the  left  pigtail 
(8cm)  is  ungrounded  (open)  and  the  right  pigtail  (8cm)  is  grounded  (short) . 
Note  that  even  though  the  left  pigtail  wire  is  ungrounded,  the  program 
assumes  the  pigtail  wire  is  present. 

Note  that  for  the  unshielded  wire  in  Fig.  5-5 (a),  there  is  a  slight 


change  in  the  response  from  a  20dB/decade  (linear  with  frequency)  response 
between  3  kHz  and  10  kHz.  Prior  to  this  experiment  the  author  has  observed 
the  coupling  to  similar  unshielded  wires  to  decrease  uniformly  at  20  dB/decade 
at  the  lower  frequencies  where  the  line  is  electrically  very  short.  This 
observed  behavior  is,  as  yet,  unexplainable  in  simple  terms  but  SHIELD 
predicts  this  result  very  accurately. 

This  experiment  tends  to  confirm  that  the  SHIELD  program  (1)  is  working 
as  intended  and  has  no  major  programming  errors,  and  (2)  provides  realistic 
predictions  of  crosstalk  for  shielded  cables.  The  user  is  cautioned,  however, 
that  if  the  cross-sectional  dimensions  (wire  spacing  and  height  above  ground) 
are  not  (1)  known  and/ or  (2)  controlled  such  as  occur  in  many  types  of  wiring 
harnesses  (which  we  refer  to  as  random  bundles)  then  the  SHIELD  program 
cannot  be  expected  to  provide  adequate  predictions  [28].  However,  in  con¬ 
trolled  characteristic  cables  such  as  the  flatpack  type,  accurate  predictions 
are  achievable.  Controlled  characteristic  flatpack  type  cables  are  presently 
being  marketed  which  contain  shielded  cables,  twisted  pairs,  etc.  For  these 
types  of  cables,  the  SHIELD  program  can  be  expected  to  provide  accurate 
predictions. 
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VI.  Summary 

This  report  contains  the  description  and  verification  of  a  digital 
computer  program,  SHIELD,  to  be  used  in  the  prediction  of  crosstalk  in  trans¬ 
mission  lines  consisting  of  U  unshielded  wires  and/or  S  shielded  cables.  The 
line  may  be  above  a  ground  plane  (TYPE  1)  or  within  an  overall,  circular, 
cylindrical  shield  (TYPE  2).  Each  shielded  cable  consists  of  a  circular, 
cylindrical  shield  which  may  be  solid  or  braided  and  a  wire  (the  shielded 
wire)  located  concentrically  on  the  axis  of  the  shield.  All  wires  may  be 
stranded  and  all  conductors  are  treated  as  imperfect  conductors;  that  is, 
their  per-unit-length  impedances  are  nonzero.  Through-braid  coupling  for 
braided  shields  as  well  as  diffusion  for  both  types  are  included  in 
the  model.  The  shielded  cables  may  have  exposed  sections  at  either  end 
(pigtail  sections)  in  which  the  shielded  wire  is  not  covered  by  the  shield. 
Over  these  pigtail  sections,  a  pigtail  wire,  parallel  to  the  shielded  wire, 
connects  the  shield  to  the  reference  conductor  at  that  end  via  either  a 
short  circuit  or  an  open  circuit.  These  pigtail  sections  are  included  in 
the  representation  to  simulate  the  common  practice  of  terminating  a  shielded 
cable  in  a  connector  via  these  pigtail  wires. 

The  multiconductor  transmission  line  equations  are  solved  for  sinusoidal, 
steady-state  excitation  of  the  line.  The  terminal  networks  at  the  two  ends 
of  the  line  are  modeled  by  a  generalized  Norton  equivalent  representation  in 
which  a  current  source  and  an  admittance  are  connected  between  the  reference 
conductor  and  the  ends  of  the  U  unshielded  wires  and  S  shielded  wires.  An 
admittance  is  connected  between  the  ends  of  each  pair  of  these  U+S  wires.  Any 
of  these  admittances  may  be  absent,  i.e.,  set  equal  to  zero. 


The  program  is  rather  expensive  to  run.  For  example  when  implemented  on 


an  IBM  370/165  digital  computer  using  the  WATFIV  compiler,  the  program 
required  3.82  sec.  compilation  time  and  84.35  sec.  execution  time  to  compute 
the  response  at  11  frequencies  for  the  example  in  Section  V  consisting  of  2 
unshielded  wires  and  3  braided-shield  cables.  This  represents  an  approximate 
per-f requency  computation  time  of  8  sec.  The  required  storage  for  this 
problem  was  115192  bytes  for  the  object  deck  and  32400  bytes  for  the  arrays. 

One  reason  for  this  rather  lengthy  run  time  is  a  result  of  the  decision 
to  include  the  capability  for  modeling  pigtails.  The  program  divides  the 
overall  line  into  cross-sectionally  uniform  sections.  Then  the  chain  para¬ 
meter  matrices  for  each  of  these  uniform  sections  of  line  are  determined  and 
the  overall  chain  parameter  matrix  of  the  line  is  determined  as  the  product 
(in  the  appropriate  order)  of  these  chain  parameter  matrices  of  the  individual, 
uniform  sections.  For  a  line  consisting  of  U  unshielded  wires  and  S  shielded 
cables,  each  chain  parameter  matrix  is  2(U+2S)  x  2(U+2S).  For  the  example  in 
section  V  consisting  of  2  unshielded  wires  and  3  shielded  cables,  each  chain 
parameter  matrix  was  16  x  16.  Thus  for  a  line  consisting  of  N  sections  we 
must,  at  the  least,  multiply  N  such  matrices  together.  Once  the  overall  chain 
parameter  matrix  of  the  entire  line  is  determined  in  this  fashion  we  must  in¬ 
corporate  the  terminal  conditions  which  requires  several  matrix  multiplications 
and  the  solution  of  2(U+S)  simultaneous,  complex  equations.  This  process  must 
be  repeated  for  each  frequency. 

Not  only  must  these  chain  parameter  matrices  for  each  section  be  multi¬ 
plied  together  but  they  must  be  computed.  In  determining  these  chain  parameter 
matrices  of  each  section  we  must  compute  the  eigenvalues  and  eigenvectors  of 
a  (U+2S)  complex  matrix.  This  is  a  result  of  the  need  to  consider  all  con¬ 
ductors  as  being  imperfect.  If  we  had  assumed  all  conductors  as  being 


108 


perfect  conductors,  the  eigenvectors  and  eigenvalues  would  not  need  to  be 
computed.  But  if  the  finite  conductivity  of  the  shields  are  not  considered, 
i.e.,  the  shields  are  considered  as  perfect  conductors,  no  coupling  would 
occur  to  the  interior  shielded  wires  for  solid  shields  -  clearly  an  unrealistic 
result.  For  braided  shields,  only  through-braid  coupling  would  occur  but  this 
is  almost  insignificant  for  many  typical  cables.  In  other  situations  one  can 
assume  perfect  conductors.  For  shielded  cables,  one  cannot  if  one  wishes  to 
obtain  predictions  approaching  any  degree  of  realism.  Once  the  decision  to 
include  shield  conductivity  is  made  we  may  as  well  include  the  finite  con¬ 
ductivity  of  all  other  conductors  -  very  little  additional  processing  is 
required . 

Consequently,  even  though  the  program  consumes  quite  a  bit  of  computation 
time  the  processing  time  probably  cannot  be  reduced  substantially  and  still 
have  the  program  handle  all  the  intended  functions  and  structures  that  it  now 
handles. 

The  program  (the  single  precision  version)  should  be  easily  implementable 
on  other  machines  since  it  is  written  in  standard  Fortran  IV  language.  For 
lines  with  known  and  controlled,  uniform  cross  sections,  the  program  provides 
accurate  predictions.  For  nonuniform  or  random  type  cable  bundles  such  as 
many  typical  wire  harnesses  where  the  relative  wire  positions  are  unknown  and 
uncontrolled  along  the  cable  length,  one  cannot  necessarily  expect  the  program 
to  provide  accurate  predictions. 
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PROGRAM  SHIELD 

(FORTRAN  TV,  DOUBLE  PRECISION) 

WPITT-N  BY 

CLAYTCN  R.  PAUL 

DEPAh  TMENT  OF  ELECTRICAL  ENGINEERING 
UNIVERSITY  OF  KENTUCKY 
LEXINGTON,  KENTUCKY  40506 

A  DIGITAL  COMPUTER  PRO  GF A  M  TO  COMPUTE  THE  TERMINAL  VC  L  T  AGES  OF 
A  ML LTICCNDUCTCR  TRANSMISSION  LINE  WHICH  NAY  CONSIST  OF 
UNSHIELDED  WIRES  AND  SHIELDED  WIRES.  THE  SHIELDED  WIPES 
NAY  HAVE  PIGTAIL  SECTIONS  AND  THE  PIGTAIL  SFCTICt:  LFNGTH3 
NEED  NOT  13 E  THF  SANE. 

THE  DISTRIBUTED  PARAMETER,  MU LTIC CNDUC1CR  TRANSMISSION  LINE 
EQUATIONS  ARE  SOLVED  FCR  STEADY-STATE,  SINUSOIDAL  EXCITATION 
OF  THE  LINE. 

THE  LINE  CONSISTS  CF  NU  UNSHIELDED  WIRES  AND  NS  SHIELDED  WIRES. 
THE  WIRES  MAY  EE  ABOVE  A  GROUND  PLANE  OR  WITHIN  AN  OVERALL, 
CIRCULAR,  CYLINDRICAL  SHIELD.  THE  IMEEDANCES  CF  ALL  CONDUCTORS 
ARE  INCORPORATED  INTC  THE  MODEL.  THE  SHIELDED  WIRES  N*Y  hAVE 
SOLID  OR  BRAIDED  SHIELDS.  THROUGH  BRAID  COUPLING  FOR  BRAIDFD 
SHIELDS  IS  INCORPORATE!  INTO  THE  MODEl. 

FUNCTION  SUBPROGRAMS  USED  (AND  SUPPLIED) : 

1 5  1 ,  L.N  1,LS2,LM2,STB,LTB,  ZWW,  ZDE  ,  2SB  ,ZD  S,  2  SS 

SUBROUTINES  USED  (AND  SUPPLIED): 

MULTC,SCAP,INDUCT,PHI,ADMADD,IMPADD 

SUBROUTINES  FROM  THE  IBSL  PACKAGE  USED: 

LEQT1C,  EIGCC 

*************************  *******************************  ************ 

IMPLICIT  REAL*8  (A-H,0-Z) 

REA  1*8  L, MU, MU02FI, MU04 PI , L MAXR , LMIN F, LMAXL , l.MINL  ,  LS  1 ,  LM  1 ,  LS2  , 
1LM2,LTRAN,LT3,LTR,LTL,L3HLD,I0R,I0I, ILR.ILI 
COMPLEX* 16  ZG,JOMEGA,ONEC,ZEROC,SUMC,SUMC1,SUMC2  ,XJ  ,ZWW  , 
1ZSS,ZDS,ZSE,ZDB,ZSELF,ZCIF,VWL,VWR,VSWL,VSWR 

************************************************************** *c 

c 


c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

INTEGER  TNT,TNUPS 

C 

IND  (NT, NT)  ,S  (NT, NT)  ,WK  (2*NT  (NT ♦ 1 )  )  C 

C 


THE  FOLLOWING  MATRICES  AND  VECTORS  SHOULD  HAVE  THE 
INDICATED  DIMENSIONS  WHERE 

NU=  NUM E  F R  OF  UNSHIELDED  WIRES 
NS=NUMEER  OF  SHIELDED  WIRES 
NUP  S  =  NU  *NS 

NT -N U*  NS*NS=TOT  AL  NUMBER  OP  CONDUCTORS 

TNT=2*NT 

IN  UP  S- 2  *  ( N U  ♦  N  3 ) 


♦**SHLDC001 
JHLD0002 
SHLD0003 
SHLDC004 
S  MLD0005 
SHLD0006 
SIILD  0007 
S  HL  0000  B 
3HLD0009 
SHLD001 0 
SHLDOO 1 1 
StiLDOO  12 
SH  LOCO  1  3 
S  IILDOO  1  4 
SHLDOO  16 
SHLD0016 
SIILD0017 
ShLDOO  1  3 
SHLD0019 
SHLD0020 
SHLD002  1 
ShLD0022 
SHLD0023 
SHLDOO 24 
SHLD0025 
SHLD0026 
SHLDOO  2  7 
SHLD0028 
SHLD0J29 
ShLDOO  30 
SHLDC031 
SHLDOO  3  2 
SHLDOO  33 
SHLD0034 
SHLDOO  36 
SHLDOO 36 
*  **SHLD  00  37 
SHLDOO 38 
SHLDOO  39 
SHLD0040 
SIILD004  1 
SHLD0042 
SHLD0043 
SHLD0044 
SHLD0046 
SHLD0046 
SHLD0047 
5HLD0048 
SHLD0049 
SHLD0050 
S1ILD006  1 
SHLD0052 
SHLD0053 
SHLD0054 
SHLDC055 
S  HLD0066 
SHLD0067 
SULD0056 
3  HLD0059 
SHLD0060 
SIILD  006 1 
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c* **************************************************** **********C 
REAl*8  IND  {  8,  8)  ,  S3  (  8,  8),HK(  144) 

****************************************************** *********Q 

c 

c 

Z  (Ml  ,NT)  ,PHIB  (TNT, TNT)  ,  HA  (NT)  ,TZ(NT,NI)  ,T(NT,NT)  ,  GAM  (NT)  ,  C 


A (NT, NT)  ,TI(NT,NT) ,  Y (NT, NT) ,YINV  (NT, NT) , B  (NT  ,  NT) ,  C 

Phil  (TNT, TNT) ,C(NT,NT)  ,D(NT,NT) ,EP (NT) ,  C 

EN(HT) , PHIL (TNT, TNT) ,TPHI (TNT, TNT) ,10 (HOPS) ,IL(NUPS) ,  C 

TO ( NUPS, HOPS)  , YL  (NUPS, NOPS)  C 

C 


C 

***************************************************************C 

COMPLEX*  16  Z(  8,  8)  ,  PHIB  ( 16,  16)  ,  MA  (  8),YZ(  8,  8), 

IT  (  8,  8),GAH(  8)  ,  A  (  8,  8),TI(  8,  8),Y(  8,  8),YINV(  8,  8), 

2B(  8,  8)  ,PHIT(16,16)  ,C  (  8,  8)  ,D  (  8,  8),EP(  8), 

3EN(  8)  ,PHIl(16,16)  ,TPH1(16,16),I0(  5),IL(  5), 

4Y0  (  5,  5) , YL  (  5,  5) 

***************************************************** **********C 

C 

C 

THE  FOLLOWING  ARB  OF  DIHENS ION  NO  C 

C 

NOTE:  IF  NU=0,  SET  ALL  DIMENSIONS*!  C 

C 

C 

***************************************************** **********Q 

INTEGER  W(2)  ,NWST(2) 

REAL *8  YH  (2) ,ZB(2)  , BN  (2)  , RNST  (2) , SIGNST  (2) 

COMPLEX*16  ZHHV  (  2) 

*********  ************************************************  ******c 

C 

c 

THE  FOLLOHING  ARE  OF  DIMENSION  NS  C 

C 

NOTE:  IF  NS=0,  SET  ALL  DIMENSIONS*!  C 

C 

C 

*********************************************** ****************c 

INTEGER  S(  3)  ,  STIPE  (  3)  ,R«H  (  3),BELT(  3),NPB(  3),NPl(  3), 

1 NPB (  3) ,IGNDL (  3) ,  IGNDB  (  3) ,KEY  (  3) 

R  E  AL*8  YS  (  3 )  ,  ZS  (  3)  ,RS  (  3),ERS(  3),RNH(  3),RWHST(  3),SIGHHS( 
ITS  (  3)  ,  SIGS  (  3)  ,  RB  (  3)  ,  SIGB  (  3)  ,THB  (  3),LFL(  3),RPL(  3),THPI( 
2RPB  l  (  3)  ,SIGPHL  (  3)  ,  RPHSTL  (  3),LPR(  3),RPB(  3),THPB(  3),RPWR( 
3BPBSTR (  3)  , SIGPHR (  3),LPT(  3),LTY(  3),STY(  3), 

4Y  PL  (  3)  ,  ZPL  (  3)  ,  YPR  (  3)  ,ZPR  (  3) 

COMPIE  X*  16  ZSV(  3)  ,  ZDV  (  3),Z*PLV{  3),ZWPBV(  3),Z»NHV(  3) 


END  DIMENSION  STATEMENTS 


COMMON  /RCON/  ZERO, ONE ,TNC, THREE , FOUR , HO, B002 PI , M004 PI , 
1PI,RADEG, ERTE, SIGCOP , R , F, V2 
COMMON  /C CON/  Z EROC,ON EC, XJ, JOHEG A 
ZERO-O. DO 
ONE= 1 , DO 
T»0  =  2.  DO 
TBREE=3. DO 
P00B=4#  DO 
CHTM=2, 84D-5 


S  HLD0062 
SHLD0063 
SHLD0064 
SHLD0065 
SHLD0066 
SHLD0067 
S  HLD0068 
SHLD0069 
SHLD007C 
SHLD0071 
SHLD0072 
SHLD0073 
SHLD0074 
SHLD0075 
SHLD  0076 
SHLD0077 
SHLD00  78 
SHLD0079 
SHLD0080 
SHLD008 1 
SHLD0082 
SHLD0083 
SHLD0084 
SHLD008S 
SHLD0086 
SHLD0087 
SHLD0088 
SHLD0089 
SHLDOO  90 
SHLDC091 
SHLD0092 
SHLD0093 
SHLD0094 
SHLD0095 
SHLD0096 
SHLD0097 
SHLD0098 
SHLD0099 
SHLD01 00 
SHLD0101 
SHLD0102 
3)  ,  SHLD  01  03 
3) ,  SHLD0104 
3),  SHLD0105 
SHLD 01 06 
SHLD0107 
SHLD0108 
SHLD01 09 
SHLDO 1 1 0 
SHLD01 1 1 
SHLD01 12 
SHLD01 13 
S  SHLDO 114 
SHLD0115 
SHLD0116 
SHLD01 17 
SHLD01 18 
SHLD01 19 
SHLD0120 
SKLD0121 
SHLD0122 


HU02PI=2.D-7 
BOO 4 FI* 1. D-7 
ONE  8  0- 1 80. DO 
V=2.997925D8 
SIGCOP=5. 8D7 
ONEC  =  DCH  PLX (ONE, ZERO) 

Z  E8 CC= DCS PLX (ZEBO, ZEBO) 
X J=DCHPIX (ZEBO, ONE) 

PI=  FOII R*DATAN  (ONE) 
RADEG=ONE80/PI 
B0=EUO2PI*T»O*PI 
V  2=V*V 

EST  E=ONF/ (NU  * V  2) 


REAC  INPUT  DATA 


REAC(5,1)  NTYPE,  L,NU,  NS 

1  rCR**AT  (I1/E10.3/I2/I2) 

■  BITE  (6,2) 

2  FORNAT ( '  1  •) 

IP  (NTYPE. EQ. 1)  GO  TO  3 
IF  (NTYPE. EQ. 2)  GO  TO  3 
GO  TO  220 

3  CONTINUE 

IF  (NU. EQ.O.AND.NS.EQ.O)  GO  TO  220 

NT=NU+NS*NS 

NUPS=NU*NS 

TNUPS=NUPS*NUPS 

TNT  =  NT  *-NT 

N  NU  =  NU 

IF (NU. EQ. 0)  NN  U=  1 
NNS=NS 

IF(NS.EQ.O)  NNS  =  1 

■  BITE  (6,4)  N U , NS , L 

4  PORH AT  (50X,* PROGRAM  SHIELD* /////4 7 X, 12 , 

1*  UNSHIELDED  BIRES • ,//48X , 12, 

2*  SHIELDED  HIBES*// 

342X, 'LINE  LENGTH  (METERS)*  ' ,  1PE 1  1 . 3///) 

IF (NTYFE. EQ. 1)  GO  TO  7 
READ  (5,5)  R 

5  PC8HAT(10X,E10.3) 

HRI TE  (6,6)  R 

6  PORHAT (10X,* BEFERENCE  CONDUCTOR  IS  AN  OVER ALL , Cl RCULAR , CYL 
1  SHIELD  WITH  INTERIOR  RADIUS*  ',1EE10.3,*  (HETEBS) »///) 

GO  TO  9 

7  CONTINUE 

■  RITE  (6, 8) 

8  FORNAT (40X,* REFERENCE  CONDUCTOR  IS  A  GROUND  PLANE*///) 

9  CONTINUE 

IF(NU. EQ. 0)  GO  TO  18 
HRI T  E (6 , 1 0) 

10  FORNAT  (/////45X,* DATA  FOR  THE  UNSHIELDED  HIBES*////) 

■  RITE(6.  11) 

11  FORNAT  (1 X, *  YH=  HIRE  HEIGHT  ABOVB  GROUND  (METERS)  OR  RADIAL 
IE  PBCH  CENTER  OF  SHIELD  (HETERS) '/ 

21 X,  •  ZB*  HORIZONTAL  COORDINATE  (HETEBS)  OR  ANGULAR  POSITION 
JS) * / IX , • HB*  HIRE  RADIUS  (NILS)*/ 

4  1 X, *  NH ST=  NUHBER  OF  HIRE  STRANDS*/ 

5 IX, • RH ST=  RADIUS  OP  HIRE  STBANDS  (HILS)  ’/ 


SHLD0123 
SHLD0 124 
SHLD0125 
SHLDO 1 2b 
SHLDO 1 27 
SHLD  0128 
SHLD0129 
SHLDO  1 30 
SHLD  0 1 31 
SHLDO 132 
SHLDO 1 3  3 
SHLD01 34 
SHLDO 1 35 
SHLDO 1 3b 
SHLD0137 
SHLDO 1 38 
SHLDO 1 39 
SHLD0140 
SHLDO 14 1 
SHLDO 142 
SHLDO 1 43 
SHLDO  144 
SHLDO 1 45 
SliLDOl  46 
SHLDO 147 
SHLDO 1 4  8 
SHLD0149 
SHLDO 1 50 
SHLDO 151 
SHLDO 1 52 
SHLDO 153 
SHLDO 1 54 
SHLDO 1 55 
SHLDO 156 
SHLDO  1 57 
SHLD0158 
SHLDO 159 
SHLDO 160 
SIILD01 6 1 
SHLDO 162 
SHLDO 163 
SHLD0164 
SHLDO  165 
SHLDO  166 
INDRICALSHLD0167 
SHLDO  168 
SHLDO 169 
SHLD017C 
SHLD0171 
SHLDO 1 72 
SHLDO 1 73 
SHLDO  1 74 
SHLDO 1 75 
SHLD  0 1 76 
SHLD0 1 7  7 
DISTANCSHLD0178 
SHLD  0 1 79 
(DLGREESHLD0180 
SHLDO 1 8 1 
SHLDO 1 82 
sHt.nniAT 


t- 


61X, 'SIGWST=  CONDUCTIVITY  OF  HIRE  STRANDS  (RELATIVE  TC  COPPER 

7) 

HBITE  (6,  12) 

12  FOB  NAT  (IX, 'WIRE' ,  6  X  , '  Y  W  ' ,  16X , ' Z W ' , 16X , • RH' , 1 2X, • N WST • , 3X , • H W 
116X, 'SIGWST'/) 

DO  17  1=1, NU 

READ  (5,13)  H (I) ,YB (I) , ZH  (I) 

13  FOBS  AT  (I2,8X,E10.3,10X,E10. 3) 

RFAC  (5,  14)  BN  ( I)  ,NNST(I)  ,RNST(I)  ,SIGNST(I) 

1*  FOHNAT(10X,E10.3,13X,I2,15X,E10.3, 10X,E10.3) 

N BITE  (6,15)  W (I) , YW (I) ,ZH  (I) ,RW  (I ) , NHST  (I) ,RNST (I) , SIGWST (I) 

15  FORNAT  (2X,I2,3X,  1PE10. 3, 8X, 1  PE  1 0. 3 , 8 X  ,  1  PE  1 0. 3 ,9 X , 12  , 
15X,1PE10.3,11X.1PE10.3) 

RH  (I) =RH  (I) *CHTM 
RNST  (1)  =RBST  ( I)  *C(1TH 
IF  (NTYPE. EQ. 1)  GO  TO  16 
ZH  ( I ) =ZW (I)/RADEG 

16  CONTINUE 

17  CONTINUE 

18  CONTINUE 

IF (NS. EQ. 0)  GO  TO  58 
WRITE (6 , 1 9) 

19  FORNAT  (/////45X,' DATA  FOR  THE  SHIELDED  HIRES'////) 

WRITE  (6, 20) 

20  FOR  S AT  (1 X, • YS=  SHIELD  HEIGHT  ABOVE  GROUND  (METERS)  OR  RADIAL 
1NCE  FRCN  CENTER  OF  SHIELD  (NETERS)'/ 

2 1 X,  •  ZS=  HORIZONTAL  COORDINATE  (HETERS)  OR  ANGULAR  POSITION  ( 
3S)'/1X,'RS=  INTERIOR  SHIELD  RADIUS  (NETEBS) •/ 

4  IX,  '  ERS=  RELATIVE  PER (1ITTI VITY  OF  INTERIOR  DIELECTRIC'/ 

51X, *RH=  SHIELDED  HIRE  RADIOS  (NILS) •/ 

61 X, ' NHST=  NUMBER  OF  SHIELDED  HIRE  STRANDS'/ 

7 1 X, ' RWST=  RADIUS  OF  SHIELDED  HI  BE  STRANDS  (NILS)'/ 

8 IX, • SIGWST=  CONDUCTIVITY  OF  SHIELDED  HIRE  STRANDS  (RELATIVE 
9PEB) •///) 

H  RITE  (6, 2 1 ) 

21  FORM  AT (2X,« SHIELD • ,8X,'YS* ,13X,«ZS',13X,'RS',13X,'ERS',12X,’ 
11  OX,  'NHST',9X,'RHST'  ,  1 1 X , ' S IGHST* ,//) 

NTHER=0 
N  ST  Y  P=  0 
DO  32  1=1, NS 

READ  (5, 22)  S < I)  , STYPE ( I) , YS (I) , ZS  (I) 

22  FCHPAT(I2»7X,I1, 10X, E 10. 3, 1  OX ,E  1  0. 3) 

IF  (STYPE  (I) . EQ. 1)  GO  TC  23 

IF (STYPE ( I) . EQ. 2)  GO  TO  23 
N  ST  Y  P=  1 

23  CONTINUE 

READ  (5,24)  RS(I)  , ERS  (I) ,RWH  (I) , NHH  (I) ,RHHSI (I) , SIGHHS (I) 

24  FORNAT  (10X,3(E10.3) , 10X.4X, I2,4X,2  (E  10.3) ) 

ITYEE=STYPE (I) 

GO  TO  (25,27) ,  ITYPE 

25  CONTINUE 

READ  (5 , 26 )  TS  (I) ,S IGS (I) 

26  PORNAT(10X,2(E10.3) ) 

GO  TO  29 

27  CONTINUE 

READ (5, 28)  RB(1)  , SIGB  (I)  , THB  (I)  ,BELT  (I)  , HPB(I) 

28  FORNAT  ( 10X, 3 (E 10. 3)  , 3X,I2,8X,I2) 

IF  (THB  (I) .LT. ZERO. OR. THE (I)  .GT.4  5.D0)  NTHER= 1 

29  CONTINUE 

READ  (5,30)  LPL(I) ,BPL(I)  ,1HPL(I)  , FPWl(I)  ,NPl(I)  ,P PHST l  (I) , 

1 SIGF8L  (I) 


) •///  SHLD0184 
SHLD0185 
SHLDO 18b 
SI',  SHLDO 187 
SHLD01 88 
SHLDO 1 89 
SHLD0190 
SHLDC191 
SHLDO 1 92 
SHLD0 1 9  3 
SHLDC194 
SHL  D0 1 95 
SHLD0196 
SHLD0197 
SHLDO 1 98 
SHLDO 199 
SHLD0200 
SHLD020  1 
SHLD0202 
SHLD0203 
SHLD0204 
SHLD0205 
SHLD0206 
SHLD0207 
DISTASHLD0208 
SHLD0209 
DEGREES HLD0210 
SHLD021  1 
SHLD0212 
SHLD02  13 
SHLDO  2  1 4 
SHLD0215 
TC  COPSHL DO 2 16 
SHLD0217 
S1ILD.02  1  8 
RH',  SHLD0219 
SHLD0220 
SHLD0221 
SHLD0222 
SHLD0223 
SHLD0224 
SHLD0225 
SHLD0226 
SHLD0227 
SHLD0228 
SHLD0229 
SHLD0230 
SHLDO  231 
SHLD0232 
SHLD0233 
SHLD0234 
SHLD0235 
SHLD0236 
SHLD0237 
3HLD0238 
SHLD0239 
SHLD0240 
SHLD024  1 
SHLD0242 
SHLD0243 
SH  LD  02  4  4 
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»- 

►V 

• 

30  PORNAT  (10X,3(E10.3)/10X,E10„3,3X,I2,5X,2(E10.3)) 

SHLD0245 

• 

RFAE(5,31)  LPR(I)  ,RPR(I)  ,THPR(I)  ,BPWR(I)  ,NPR(I)  ,RPWSTR(I)  , 

3HLD0246 

1SIGPHR  (I) 

SHLD024  7 

31  PORN  AT  ( 1  OX, 3  (El  0. 3) /10X,E10.3,JX,I2,5X,2(E10. 3)  ) 

5HLD0248 

•  // 

32  CONTINOE 

SHLD0249 

IF  (NTHER. EQ. 1 )  GO  TO  220 

SHLD0250 

* 

IP  (NSTYP. EQ.  1 )  GO  TO  220 

SHLD0251 

*  •  •  ■  r 

v  m  . 

DO  35  1=1, NS 

SHLD0252 

• 

WRITE  (6,33)  S  (I)  , IS  (I)  ,ZS  (I)  ,  BS  (I)  ,EBS  (I)  ,  RW H  (I)  ,  NBU  (I)  , 

RWSHLD0253 

»’  - . 

1 BST  (I)  , SIGHHS  (I) 

SHLD0254 

•  - 

3  3  PORN AT ( 3X, 12, 7X,  1PE 1 0.  3, 5X,  1  PEI  0 . 3 , 5  X  ,  1  PE  10 . 3  ,5  X,  1  PE  10.  3  , 

SHLD0255 

1 5  X,  1PE10.3,7X,I2,7X,1PE10.3,5X, 1PE10.3) 

SHLD0256 

HUH  (I) =  RHH  (I) *CNTN 

SHLD0257 

RHHST(I) =  RH  HST ( I) *CNTH 

SHLD0258 

IP  (NTYPE.EQ.1)  GO  TO  34 

SHLD0259 

'■<§  ' 

ZS(I)=ZS(I)/RADEG 

SHLD0260 

34  CONTINUE 

SHLD0261 

35  CONTINUE 

SHLD0262 

*  *  -  * 

■  BITE  (6,36) 

SHLD02  6  3 

36  PORNAT  (///// 35 X, 'SHIELD  CHARACTERISTICS'//) 

SIILD0264 

■  81 TE  (6 ,37) 

SHLD0265 

- 

37  PORNAT  (IX,  'TS  =  SOLID  SHELD  THICKNESS  (NILS)'/ 

SHLD0266 

1 IX, • SIGS=  SOLID  SHIELD  CONDUCTIVITY  (RELATIVE  TO  COPPFR) •/ 

SHLD0267 

21 X, • 8B=  RADIUS  CP  BRAIL  HIRES  (NILS)'/ 

SHLD0268 

3 IX, ' SIGB=  BRAID  HIRE  CONDUCTIVITY  (RELATIVE  TC  COPPER)'/ 

SHLD0269 

4 IX, *THB=  BRAIL  HEAVE  ANGLE  (DEGREES)'/ 

SHLD0270 

51 X, ' BELTS=  NUNBER  OF  BELTS  IN  BRAID'/ 

3HLD0271 

6 IX, *  HPB=  NUNBER  OF  HIRES  PER  BELT'///) 

SHLD0272 

HRITF (6, 38) 

SHLD0273 

38  FCRHAT (2X, ■ S HI  ELD • , 4 X, 'TYPE* , 8X , ' TS* , 13X, ' SIGS' , 1 4X , • RB' , 1 6 X , 

SHLD0274 

! 

1'SIGE* , 1 1X, 'THB',1 OX, 'BELTS', 10X,  'HPB  *//) 

SHLD0275 

i 

DC  43  1=1, NS 

SHLD0276 

IF(STYPE  (I) . EQ.2)  GO  TC  40 

SHLD0277 

”  j 

■  RITE  (6,  39)  S(I)  ,TS(1)  ,SIGS  (I) 

SHLD0278 

39  PORNAT  (3X, 12, 6X,' SOLID* ,4X, 1PE10. 3,5X  ,1  PE  10. 3) 

SHLD027 9 

TS  (I) =TS  (I) *CNTN 

SHLDO20O 

RS(I)=RS(I)*TS(I) 

3HLD0281 

GC  TO  42 

3HLD0282 

40  CONTINUE 

SHLD0283 

•  —  / 

■  RITE  (6,41)  S  (I)  ,RB(I)  ,SIGB  (I)  , TUB  (I)  ,BELT  (I)  ,HPB(I) 

SHLD0284 

-  s 

4  1  PORNAT(3X,I2,5X,'BBAIDED',35X,1PE10.3,8X, 1PE10. 3 , 5 X , 1  PE  10 . 3 , 9X 

,  SHLD0285 

112, 12X,I2) 

SHLD02B6 

RE(I)=RB(I)*CNTN 

SHLDO207 

TS(I)=THO*RB(I) 

SHLDO280 

BS(I)=RS(I)*TS(I) 

SHLD0289 

• 

THB  (I) =THE  (I) /RADEG 

SHLDC290 

•  -1 

42  CONTINUE 

3BLU0291 

43  CONTINUE 

SHLD0292 

■  BITE  (6,44) 

SHLD029J 

44  PORNAT  (////, 39X,' DATA  POR  THE  PIGTAILS'//) 

SHLD0294 

■RITE (6,45) 

SHLD0295 

45  PORNAT  (421, 'LEFT  PIGTAILS'//) 

SHLD029b 

■BITE (6,46) 

SHLD0297 

•  ’ 

46  PORNAT  (IX, *LPL=  LENGTH  CP  LEFT  PIGTAIL  (NETERS)'/ 

SHLD0298 

1 1 X, • BPL=  RADIAL  SEPARATION  OP  PIGTAIL  HIRE  PRCH  SHIELDED  HIRE 

(NETSHLD0299 

2ERS)  •/ 

SHLD0300 

3 1 X, *TH PL*  ANGULAR  POSITION  OP  PIGTA1I  HIRE  (DEGREES)'/ 

SHLD0301 

4 1 X, ' FPHL=  RADIUS  OF  PIGTAIL  HIRE  (NILS)'/ 

SHLD0302 

SIX,  »NP'L=  NUNBER  OP  STRANDS  IN  PIGTAIL  HIRE'/ 

SHLD0303 

6 1 X, • BPHSTL=  RADIUS  OP  PIGTAIL  HIRE  STRANDS  (NILS)'/ 

SHLD0304 

71X, • SIGPHL=  CONDUCTIVITY  OP  PIGTAIL  HIRE  STRANDS  (RELATIVE  TO 

COPPSH  LDO  305 

-  -  -4 

119 

■ 

• 

■ 

#  ' 

~  .. 

-  ^  ^ . -J 

8ER) •//) 

WRITE(6,47) 

4  7  FOR  FAT  (IX, ‘SHIELD* ,9X, • LPL' , 14X, ' BFL' ,  12X, 'THPL* , 10X, 'RPWL' , 9X 
1  *  NPL ' , 9X, 'RPWSTL' , 12X, 'SIGPWL'//) 

DC  49  1=1, NS 

IMin  E  (6,4  8)  S  (I)  ,  LPL  (I)  ,  RPL  (I)  ,  T  H  I  L  (I)  ,  RPW  L  (I)  ,  SPL  (I)  , 

1WSTL  (I) , S IGP  WL (I) 

4  8  F0RNAT(3X,I2,7X,  1t>El0.3,7X,1PE10.3,5X,1PE10.3,4X,1PF10.3,7X, 

1I2,7X, 1PE10.3,8X, 1  PE  10. 3) 

RPWL  (I)  =R?WL  (I)  *CNTN 
R  PWSTL  (I)  =RPWSTL ( I)  *CNIN 
THPL  (I) =THPL (I) /RADEG 

49  CONTINUE 
MBI T  E  (  6 , 50) 

50  FCRNAT  (////, 42X,' RIGHT  PIGTAILS'//) 

WRITE  (6,51) 

51  FCRNAT (1X,'LPR=  IENGTH  CF  RIGHT  PIGTAIL  (NFTE  PS)  1 / 

1 1 X, ' RPR=  RADIAL  SEPARATION  OF  PIGTAIL  WIRE  FRCN  SHIELDED  WIRE 
2ERS)  •/ 

3 1 X, • TH  PR=  ANGULAR  POSITION  OF  PIGTAIl  WIRE  (DEGREES)'/ 

4 1 X, ' BPWR=  RADIUS  OF  PIGTAIL  WIRE  (HIIS) '/ 

5 IX, • NPR=  HUN  D  ER  OF  STRANDS  IN  PIGTAIL  WIRE'/ 

6 1 X,  •  EPWSTR=  RADIUS  OF  FIGTAIL  WIRE  STRANDS  (NILS)'/ 

71 X, ' SIGPWF=  CONDUCTIVITY  OF  PIGTAIL  WIRE  STRANDS  RELATIVE  TO 
8R)  •//) 

WBITE{6,52) 

52  FOR  RAT  (IX, 'SHI  ELD' , 9 X, ' LPR* , 1 4X , • EPB • , 1 2X, • THPR • , 1  OX , 

1' RPWB' ,9X, 'NPH' ,9X,'RPWSTR* , 12X, • SIGPWR'//) 

DC  54  1=1, NS 

WRITE (6,53)  S  (I) ,LPR  (I) ,RPH  (I)  , T  HFR  ( I)  , RPWR  ( I) , NPR(I)  , 

1WSTH  (I)  , SIGPWR  (I) 

5  3  FCR"AT(3X,I2,7X, 1PE10. 3,7X, 1PE10. 3 , 5X , 1PE 1 0 . 3 ,4 X , 1  PE  10. 3, 

17X,I2,7X,1PE10.3,3X, 1PE10.3) 

RPWR  (I) =RPWR (I)  *CMTN 
RPWSTR(I) =RPWSTH (I) *CNTN 
THPR  (I) =THPR (I) /RADEG 

54  CONTINUE 

DC  57  1=1, NS 

I F  ( NTYFE. EQ. 2)  GO  TO  55 

TPL(I)  =YS(I)  ♦RPL(I)*DSIN  (THPL(I)) 

ZPL(I)  =ZS(I)  ♦RPL  (I)  *DCOS  (THPL  (I)  ) 

YPR  (I)  =YS  (I)  ♦RPR  (!)  *DSIN  (THPR  (I)  ) 

ZPR  (I)  =ZS  (I)  *RPR  (I)  *DCOS  (THPR  (I)  ) 

GO  TO  56 

55  CONTINUE 

YPL(I)  =DSQRT  (YS(I)  **2  +  RPL  (I)  **2-TW0*YS(I)  *RPL  (I)  *DCC3  (THPL  (I)  )  ) 
ZPL(I)  =ZS(I)  ♦DATAN2(  (RPL(I)  +D3IN  (THPL  (I)  )  )  ,  (YS(I)  ♦ 

1RFL  (I) *DCOS (THPL (I) ) ) ) 

YPR  (I) =DSyRT (YS(I) **2»HPR  (I) ** 2-T KC* IS (I) *RP8 (I) *DCCS  (THPR (I) ) ) 
ZPR  (I) =ZS  (I)  ♦DAT AN 2 (  (RFR  (I)  ♦DSIN (T H PR  (I) ) ) , ( Y S  ( I) ♦ 

1RPR (I) *DCOS  (THPR  (I) ) ) ) 

56  CONTINUE 

57  CONTINUE 

58  CONTINUE 
WRITE (6,59) 

59  FORNAT (////, 33X, 'TERNINAL  SOURCE  AND  IHPECANCE  DATA'////) 
WRITE(6,60) 

60  FORNAT (IX, *IOR=  REAL  PART  OF  CURRENT  SOURCE  EFTWEEN  WIRE  AND 
1ENCE  CONDUCTOR  AT  X=0'/ 


SUL  DC  3  06 
SHLD0307 
,  3ULD0308 

SHLD0309 
SHLD031C 
R  PS11LD0  3 1 1 
SHLDC312 
SHLDO 3  13 
SHLDO 3  1  4 
SHLDO.)  15 
SHLDO  316 
SHLDO 31 7 
SHLD0318 
SHLDO  319 
SOLDO  320 
SHLD032  1 
SHLD0322 
(NETSHLD0323 
SHLD0324 
SHLD0325 
SHLD0326 
3HLD0327 
SIILD0326 
COPPE3HLD0326 
SHLDC330 
SHLD033  1 
SHLDO  3  32 
SHLD0333 
SHLDO 334 
PPSHLD0J35 
SHLD0336 
SHLDO  3  3  7 
SOLD  03  38 
SHLD0339 
SHLD0340 
SHLDO 34  1 
SHLDC342 
SIILD0343 
SliLDO  344 
SHLDC345 
SHLD0346 
SHLDO  34  7 
SHLDC343 
SliLDO  349 
SIILD035  0 
SHLDO  351 
SHLDO  3  5  2 
SHLD035 J 
SHLDO 354 
SHLD0355 
SHLD035b 
SHLDC357 
SHLD0356 
SHLD0359 
3HLD0360 
SHLD036 1 
SHLD0362 
REFERSHLD0363 
SHLD0364 


21 X, ' 101=  I  HAG  PART 
3NCUCT0R  AT  X=0'/ 


OF  CURRENT  SOURCE  E EX WEEN  WIRE  AND  REFERENCE 


COShLD  0  365 
SHLDC366 
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61 


62 


63 


66 


67 


68 


69 


4  IX,  ' YOR  =  SEAL  PART  OF  ADMITTANCE  EEISEEN  WIHE  AND  REFERENCE 
5TCR  AT  X=0'/ 

6 1 X, • Y0I  =  I  NAG  PART  OF  ADMITTANCE  E  ETSEEN  WIRE  AND 
7TCR  AT  X=0'/ 

81X, ' ILR=  REAL  PART  OF  CURRENT  SOURCE  BETNEEN  WIRE 
9NCUCTOR  AT  X=L'/ 

A 1 X, * ILI=  IMAG  PART  OF  CURRENT  SOURCE  EETWEEN  HIRE 
ENCUCTOR  AT  X  =  L'/ 

C IX, ' YLR=  BEAL  PART  OF  ADMITTANCE  EEISEEN  HIRE  AND 
DTOR  AT  X  =  L'  / 

BIX, • TL 1=  IMAG  PART  OF  ADMITTANCE  EETSEEN  HIRE  AND 
FTOR  AT  X=L'///) 

WRITE (6,61) 

FOR MAT (IX, 'WIRE', 8X, • IOR ' , 1 2X , 'IOI ' , 12X , * YOR • , 1 2X ,  •  YO I • , 1 2X , 
1'ILB',12X,'ILI',12X, 'YiR',12X,'YLI'//) 

IF(NO.EC.O)  GO  TC  65 
DO  64  1=1, NU 

READ  (5,62)  IOR ,  IOI ,  YOR  ,  YOI , ILR  ,1 II  , HR  , YU 
FOR  RAT  (8  (El  0.  3)  ) 

WRITE (6, 6 3)  I, IOR, IOI, YOR , IOI , ILR ,ILI ,  YLR , YLI 
FORMAT (2X, 12, 5X,1PR10. 3, 5X, 1PB10. 3, 5X  ,  1  PE  10. 3  ,5 X, 1 P E  10. 3 ,5X , 
11  PE  10. 3,5X,1PE10.3,5X, 1PE10.3,5X,  IP E 10. 3) 

10(1) =IOR*XJ*IOI 
II(I)=ILR*XJ*ILI 
YO (I ,1) =YOR*XJ* YOI 
YL(I,I)=YLR*XJ*YLI 
CCN7INUE 

IF(SS.EC.O)  GC  TO  73 

NISR=0 

NVSH=0 

NIS1=0 

NVSL=0 

WRI T  E (6 , 66) 

FORMAT (//, IX, 'SHIELDED  W IRE • , 3X , • IOR • , 12X, • IOI' , 1 2X , • YOR • , 1 2X 
1, 'YOI* ,12X,'ILR',12X,'ILI',12X,'YLR',12X,*ILI»//) 

NGNDR=0 
NGNDI=0 
DC  72  1=1, NS 

BEAD  (5,6  7)  I  OB  ,101 ,  YOR  ,  YOI  ,ILH  ,1  LI ,  UR  ,  III 
FORMAT (8  (El  0.3) ) 

WRITE (6, 68)  I, IOR, IOI, YCR , IOI , ILR , I II , YLR , III 
FORMAT (3X, 12, 8X,1PE10. 3,5X, 1 PE10 . 3, 5X , 1 PE10. 3 ,5 X , 1  PE  10. 3  ,5X , 
11  PE  10. 3,5X,1PE10.3,SX,1PE10.3,5X,  1PE10. 3) 

I0(I+NU)=I0R*XJ*I0I 
II (I *NU)  =ILB«-X  J*ILI 
YO (I»NU,I»NU) =Y0R>XJ* YOI 
IL  ( I *NU ,  I  +  N'J)  =YLR»XJ*YLI 
READ  (5,69)  IGNDL (I)  , IGNDB  (I) 

FORMAT  (9X.I1 ,29X,I1) 


CCNDUCSHLDQ367 
SHLD0368 

REFERENCE  COSDU CS HLDO 369 
SHLD0370 

AND  REFERENCE  C03HLD0371 
SHLD0372 

AND  REFERENCE  COSHLD0373 
SHLDC374 

REFERENCE  CCNDUCSHLD0375 
SHLD0376 

REFERENCE  CCNDU CSHL D0377 
SHLD0378 


70 


71 


IF  ( IGNDL  (I) 
IF ( IGNDL (I) 
NGNDl=1 
CONTINUE 
IF  (IGNDR(I) , 
IF  (IGNDR  (I) 
NGN  CR=  1 
CCNTINUE 
IF  (IGNDB (I) 
IF  ( IGNDR  (I)  , 
IF (IGNDL  (I) 
IF  ( IGNDI  (I) 


EQ 

EC 


D 

2) 


EQ.  1) 
EC. 2) 


EC.1) 
EQ.  2) 
EQ.  1) 
EC. 2) 


GO 

GO 


GO 

GO 


TO 

TC 


TC 

TC 


70 

70 


71 

71 


NISR=NISR ♦ 1 
N  V  SR  =N  VSR  ♦  1 
NI SL=N ISL ♦ 1 
NVSL=NVSL*1 


SHLD0379 
SHLD0380 
SHLD038  1 
SHLD0382 
SHLD0383 
SHLD0384 
SHLD0385 
SHLD0386 
SHLD0387 
SHLD0388 
SHLD0389 
SHLD0390 
SHLD039 1 
SHLD0392 
SHLD0393 
SHLD0394 
SHLD0395 
SHLD0396 
SHLDC397 
SHLDC398 
SHLD0399 
SHLD0400 
SHLD0401 
SHLD0402 
SHLD0403 
SHLD0404 
SHLD0405 
SHLD0406 
SHLD04  07 
SHLD0408 
SHLD0409 
SIILD0410 
SHLD041 1 
SHLD04 12 
SHLD04 1 3 
SHLD04 14 
SHLD0415 
SHLD0416 
SHLD04 1 7 
SHLD04 18 
SHLD04 1 9 
SHLD0420 
SHLD042 1 
SHLD0422 
SHLD0423 
SHLD0424 
SHLD0425 
SHLD0426 
SHLD0427 
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72  CONTINUE 

IF(NGNDL.EU.  1.0B.NGNDR.EQ.  1)  GO  TC  220 

73  CONTINUE 
K  1  =  NUPS-1 

DC  75  I =  1  #  K  1 
IP1=I* 1 

DO  75  J=IP1 , NUPS 

BEAD  (5, 74)  TOR ,YOl , YLR , YLI 

74  P OH  MAT  (2 OX, 2  (E10.3)  , 20 X , 2  (E 10. 3) ) 

TO  (I,J) =Y0R*XJ*Y0I 

Y  L  ( I ,  J  ) =YLR»XJ*YLI 
TO  (J,I) =  Y0  (I, J) 

75  Yl  (J,I) =  YL (I ,J) 

8BITE(6, 76) 

76  FOR  NAT  (// ,  10X, 'IMPEDANCES  BETWEEN  HIRES'//) 

IF  (NU. EC. 0)  GO  TO  80 

WRITE(6, 77) 

77  FOR  NAT (IX, 'HI  PE' ,2X, 'HIRE* , 18X, • YCR'  ,12X,' YOI' , 12 X, • YLB* ,12X, 

1 ' YLI'//) 

MKN  =  NU -  1 
DC  79  1=1, HKN 
IF1 =1*1 

DC  79  J  =  IP1,NII 

WRITE (6, 78)  I,J,YO(I,J) ,YL(I,J) 

7  8  FOB  NAT  (2 X, 12 ,4 x , 12 , 1 5X , 1  PE  1 0. 3, 5X , 1  PE  10. 3, 5X , 1PE 1 0 . 3 , 5X , 1PE 1 0 . 3) 

79  CONTINUE 

80  CONTINUE 

IF  (NU. EC.O.OR. NS. EQ.O)  GO  TC  84 
WBITE<6, 81) 

8 1  FOR  FAT  (///, IX, 'HIRE* ,2X,' SHIELDED  8 1 R E ' , 1 2X , ' YO R •  ,  1 2 X  ,  • YO I • ,  1  2X , 
1' YLB* , 12X,' YLR'//) 

DC  83  1=1, NU 
DC  83  J=1 , NS 

HRITE(6,02)  I,  J, YO  (I, J'NU)  ,YL(I,J»NU) 

82  FORMAT (2X,I2,7X,I2,  15X,1PE10.3,5X,1PE10.3,5X,1PE10.3,5X, 

1 1  PE  10.  3) 

83  CONTINUE 

84  CONTINUE 

IF  (  NS. EC. 0)  GC  TO  100 
8RITE(6,85) 

85  FORM  AT (///, IX, ' SHIELDED  8 IR E ', 2X ,' SHIELDED  HIRE’,7X, 

1' YOB  ',  19X, 'YOI' , 19X, 'YLB', 19X, • YLI'//) 

MKS=NS-1 
DC  87  1=1, MRS 
IP1=I»  1 

DC  87  J=IP1 , NS 

HRITE(6, 86)  I , J , YO  (  I* NU , J ♦ NU)  , YL  (I* NU, J*  NU) 

86  FORMAT  (5 X ,12 , 1 3X, 12,  10X , 1PE 10. 3, 1 2X , 1PE 10. 3 , 12X , 1  PE  1 0. 3 , 

1 12X,1PE10.3) 

87  CONTINUE 
8RI TE (6 , 88) 

88  FORMAT  (////, 20 X, ' SHIE  ID  TERMINATION  DATA'//) 

WRITE (6,  89) 

89  FCRNAT(1X, 'SHIELD' ,5X, 'TERMINATION  AT  X=0' , 

110X, 'TERMINATION  AT  X=L'//) 

DO  99  1=1, NS 

IF(IGNDL (I) . EQ.  1 . AN D. IGNDR ( I)  .EQ.  1)  GC  TO  90 
IF  (IGNDL  (I) . EQ. 1. AND. IGNDR(I) .EQ. 2)  GO  TO  92 
IF  ( lGNDL  (I) • EQ. 2. AND. IGNDR (I) . EQ. 1)  GO  TO  94 
IF(IGNDL(I).EQ.2.AND.IGNCR(I) .EQ.2)  GO  TO  96 

90  CONTINUE 


SHLDC428 
SHLD0429 
SHLD04  30 
SHLD0431 
SHLD0432 
SHLD0433 
SHLD04  34 
SHLD0435 
SHLD0436 
SHLDC437 
SHLD04  38 
SHLD04  39 
SHLD0440 
SHLD0441 
S.HLD04  4  2 
SHLD0443 
SHLD04  4  4 
SHLD0445 
SHLD0446 
SHLD0447 
SHLD0448 
SHLD0449 
SHLD0450 
SHLD0451 
SHLD  C4  52 
S HLD0453 
SHLD0454 
SHLD0455 
SHLD045b 
SHLD04S7 
SHLD0458 
SHLD0459 
SHLD0460 
SHLDC461 
SHLD0482 
SHLD046  3 
SHLD0464 
SHLD0465 
SHLD0466 
SHLD0467 
ShLD0468 
SHLD0469 
SHLDC470 
SHLD0471 
SHLD04  7  2 
SHLD0473 
3HLD04  74 
SHLD0475 
SHLD 04 76 
SHLD0477 
SHLD04  78 
SHLDC4  79 
SHLD0480 
SHLD048  1 
SHLD0432 
SHLDO403 
SHLD0484 
SHLDC485 
SHLD04  86 
3HLD0487 
SHLD04  88 


• 

' 

isB 

H  BITE (6,91)  I 

SHLD0489 

•  ^ 

91  FORN  AT  (3X,  12, 12X,’  SHORT  •  ,23X,»  SHORT) 

SHLD0490 

,  *. 

GC  TO  98 

SHLD049  1 

92  CONTINUE 

SHLDC492 

WRITE(6,93)  I 

3HLD049  3 

i  *  ■  * 

9  3  FORBAT  (3 X,I2 , 12 X, ’ SHORT • ,24X , • OPEN* ) 

SHLD0494 

■  "  -  *  4 

GO  TO  98 

SHLD0495 

*.  % 

94  CONTINUE 

SIILD0496 

“ 

HBITE(6,95)  I 

SHLD0497 

.  i  • 

95  FOBS  AT (3X, 12 , 1 3X , • OPEN  * ,23X , • SHOBT ' ) 

SHLD0498 

.  •  /, 

GC  TO  98 

SHLD0499 

v 

96  CONTINUE 

SHLD0500 

W  HI T  E  (6  , 97)  I 

SHLD050 1 

*-  - 

97  FORBAT (3X.I2, 13X,’OPEN* ,24X,’OPEN’) 

SHLD0502 

98  CONTINUE 

SHLD050  3 

1  ■Hi 

99  CCNIINUE 

SHLD0504 

♦ 

100  CONTINUE 

SHLD0505 

' 

DO  101  1  =  1  , N  UPS 

SHLD0506 

DO  1C1  J  =  1 ,N  UPS 

SHLD0507 

r  * 

YO {I ,J) =- YO ( I, J) 

SHLD0508 

';•% 

101  Yl (I,J)=-YL(I,J) 

SHLD0509 

DC  103  1= 1 , N  UPS 

SHLD0510 

m 

SUHC  1=ZEBOC 

SHLD0511 

i 

SUBC2=ZEBOC 

SHLD05  12 

DC  102  •!=  1  ,N  UPS 

SHLD051 3 

SUBC 1=SUBC 1* YO (I, J) 

SHLD0514 

102  SUBC2  =  SUBC2* YL  (I,J) 

SHLD05  15 

YO (1,1) =-SUBCl 

SHLD  05  1  6 

103  YL (1,1) =— SUNC2 

SHLD0517 

104  CONTINUE 

SHLD0518 

1 

C 

SHLDC519 

r 

C 

SHLD0520 

C  REA  C  FREQUENCY  AND  CONFUTE  CONDUCTOR  ISPEDANCES 

SHLD0521 

C 

SHLD0522 

• 

i 

c 

5  ML  CO  5  2  3 

, 

READ (5,105, END=225)  F 

SHLD0524 

>  ' 

• 

105  FORBAT  (E10. 3) 

SHLD0525 

OBEGA=TWO*PI*F 

SHLD0526 

JCNEGA=XJ*OMEGA 

SHLD0527 

Mi  4 

READ  (5, 106)  ZGR.ZGI 

SHLD0328 

106  PORBAT  (  2  ( E  10 .  3)  ) 

SMLD0529 

ZG=ECNPLX (ZGR.ZGI) 

SIILD05  30 

WRITE  (6,1 07)  F, ZGR.ZGI 

SHLD0531 

-  . 

107  FORBAT (////, IX, ’FREQUENCY (HZ)  =  ',1PE10.3// 

SHLD0532 

1 1 X , ' FEFERENCE  CONDUCTOF  IBPEDANCE=  ’,1PE10.3,’  ’, 

SIILD0533 

2 1  PE  1 0. 3 , '  OHBS  PER  HETEB’) 

SHLD0534 

I E ( NU. EQ.O)  GO  TO  109 

SHLD0535 

DC  108  1=1, NU 

SHLD0536 

108  ZWWV(I)=ZHW  (R  WST  (I)  ,  NWST(I)  ,SIGWST(I)) 

SHLD0537 

109  CONTINUE 

SHLD0538 

•.  ■ 

IF(NS.EQ.O)  GO  TO  113 

S1ILD0539 

DC  112  1=1, NS 

SHLD0540 

ZWWHV(T) =ZWW (RHHST (I) , NWH(I) , SIGNHS(I)) 

SMLD054  1 

• 

ZWPLV(I) =ZWW  (RPWSTL(I)  , NPL  ( I)  ,S IG EWL  (I) ) 

SHLD0542 

•1  -4 

ZWPRV  (I) =ZWW  (RPWSTR (I) ,NPR (I) ,S IG PH B (I)  ) 

SHLD0543 

IF (STYPE (I) . EQ. 2)  GO  TO  110 

SHLD0544 

ZSV  (I)  =ZSS(RS  (I)  ,TS  (I)  ,5  IGS  (I)  ) 

SHLD0545 

ZDV  (I)  =ZDS  (RS  (I)  ,TS  (I)  ,S IGS  (I) ) 

SHLD0546 

LTV  (I)  =ZE80 

3HLD0547 

*  * 

STV (I) =ZERO 

SIILD054H 

GO  TO  111 

SHLD0549 

'1 

123 

' 

110  CONTINUE 

ZSV (X) =  ZSB  (RS  (I)  , BB  (I)  ,S  IGB  (I) , T H  E  (I) , 

1BELT  (I) ,HPB  (I) ) 

ZCV(I)  =ZDB(RS(I)  ,RB  (I)  ,SIGB(I)  ,THB(I)  , 

1BELT  (I)  , WPB  (I) ) 

LTV  (I)  =LTB  (BS  (I)  ,RB(I)  ,THB  (I)  , 

1  BELT  (I)  ,  WPB  ( I) ) 

STV  (I)  =  STB  (RS  (I)  ,  RB  ( I)  ,THB(I)  , 

1BELT  (I) ,»PB  (X) , ERS  (I) ) 

111  CONTINUE 

112  CONTINUE 

113  CONTINUE 

IF(NS. EQ.O)  GO  TO  164 
C 
C 

C  COMPUTE  CHAIN  PARAMETER  MATRICES  FOR 

C  RIGHT  PIGTAIL  SECTIONS 

C 
C 

DC  114  1=1, NS 

114  LPT  (I) =LPH  (I) 

IPR  =  0 

LHAXR=ZERO 
DO  116  1=1, TNT 
DO  115  J=1,TNT 

115  PHI  R  (I  ,  J)  =ZEROC 

1 16  PHIB  (1,1) =0N  EC 

117  CONTINUE 

DC  118  1=1, NS 

118  IF  (LPT  (I) .GE.LHAXR)  LMAXB=LPT  (I) 

IF(IPR.EQ.O)  LTR=LM AXR 
IP(LMAXR.EQ.ZERO)  GO  TC  132 
IPR=IPR»1 

LMIN  R=LHAXR 
DO  119  1=1, NS 

119  IF  (LPT  (I)  -LE.LMINH.AND.LPT(I) .NE.ZERC)  IHINR  =  LPT(I) 

DO  120  1=1, NS 

KEY  (I) =1 

120  IF  (LPT  (I) .GE.LHINR)  KEY(I)=2 

CALL  INDUCT (IND,NU,NS,NT,NUPS,NTYPE,RB,YW,ZH,RS,YS,ZS, 
1RFMR,YPR,ZPR,RWH, KEY,NKU,NNS) 

DC  121  1=1, NT 
DO  121  J=  1, NT 

121  Z  (I,J) =JOMEGA*IND(I,J) *ZG 
IF(NU.EQ.O)  GO  TO  123 

DO  122  1=1, NU 

122  Z (I ,1) =Z  (1,1) +ZNNV  (I) 

123  CONTINUE 

CALL  IMPADD  (Z  , KEY  ,  NO ,  NS,  NT , NUP  S  ,ZS  V  ,ZD  V  ,  LTV  ,  ZW PR  V  ,  Z W  WH V  ,  N'!0  ,  H NS) 
CALL  SC AP (SS ,IND, NU, NS , NT, NU PS, KEY, AND, NNS) 

CALL  A  DM  ADD  (KEY, NU, NS , NT, MOPS , RS , TS, R»H,ERS, SS, STV,NNU,NNS) 

DO  125  1=1, NT 
DO  124  J= 1 , NT 
T  (I,J) =SS  (I, J) *ONEC 

124  TI  ( I , J) =Z EROC 

125  TI  (I , I) =ONEC 

CALL  LEQT1C  (T,NT,NT,TI,NT,NT,0,WA,IER) 

DO  126  1=1, NT 
DO  126  J= 1 , NT 

Y  (I, J) =JOMEG A*DREAL (TI (I,J)  ) 

126  YINV  (I, J) =SS  (I,J) /JOMEGA 


SIILD0550 
SIILD055 1 
SHLD0552 
SI1LD0553 
SHLD0554 
SHLD0555 
SHLD0556 
SHLD0557 
SHLD 05 5 8 
SHLD0559 
SHLD0560 
SHLD0561 
SHLD0562 
SHLD056  3 
SHLD  0564 
SHLD0565 
SHLD0566 
SHLD0567 
SHLD0568 
SHLD056  9 
SHLDC57C 
SHLD0571 
SHLD0572 
SHLD0573 
SHLD0574 
SHLD0575 
SHLD0576 
SHLD0577 
SHLD0578 
SHLD0579 
SHLD0580 
SHLD059  1 
SHLD0582 
SHLD0583 
SHLD0S84 
SHLD'05  85 
SHLD0586 
SHLD0587 
SHLD0588 
SHLD0589 
SHLD  05  90 
S1ILDC59  1 
SHLD0592 
SHLD0593 
SHLD0594 
SHLD0595 
SHLDC596 
SHLD0597 
SHLD0598 
SHLDC599 
SHLD0600 
SHLD060 1 
SHLD0602 
SHLDObO  3 
SHLD0604 
SHLD0605 
SHLD0606 
SHLD0607 
SHLD0608 
SHLD0609 
SHLD06  10 


124 


* 

*  k 

• 

3S 

CALI  HULTC(TI,  Y,  Z,NT,  NT,  NT) 

3HLD06  1 1 

CALI  EIGCC (TI,NT»NT,2,GAR,T,NT, WK  ,  It  B) 

S  li  L  f)  0  6  1  2 

DO  128  1=1,  NT 

S  HL  D  06  1  3 

DC  127  J=  1 , NT 

3HLD0614 

Z  II,  J)  =T  |I,  J) 

3HLD0615 

127 

T I ( I , J) =Z  EfiOC 

3  ')  L  D  0  6  1  6 

128 

TI (1,1) =0NEC 

SHLD06  1  7 

Cl 

CALI  LECT1C (Z,NT ,NT,TI ,NT, NT,0, WA, IER) 

SHLDO0I8 

• ' 

DC  129  1=1,  NT 

SHLDCb 1 9 

129 

GAR  (I)  =CD3QRT  (GAR  (  I)  ) 

3HLD0620 

CALI  PHI(PHIT,NT,TNT,L*<INR,GAN,YfYINV,T 

,TI,  EP,EN) 

SHLD062 1 

CALI  HULTC (TPHI,PHIR,PHIT,INT,TNTfTNT) 

SHLD0622 

*_ 

DC  130  1=1, TNT 

SHLD0623 

DC  130  J=1,TNT 

SHLD0624 

130 

PHIR  (I ,  J) =  TPHI ( I , J) 

SHLD0625 

lB 

DC  131  1=1, NS 

3HLD0626 

• 

131 

IF  (IPT  (I)  . NE.ZE8C)  LPT  (I)  =L PT  ( I)  - LR I N F 

SHLD0627 

L '  • 

Lf1AXR=ZERO 

SHLD0628 

GC  IC  117 

3HLD062S 

132 

CONTINOE 

SHLD0630 

•  • 

mm 

U  RIT  E (6  ,  1  33)  IPR 

SHLD0631 

*  . 

mm 

133 

FCR  R AT (// 10X , • THE  NUREER  CF  FIGHT  PIGTAIL  SECTIONS  = 

SHLD0632 

«H 

1.12//) 

SHLD0633 

c 

3HLD0634 

c 

SHLD0635 

c 

CCREUTE  CHAIN  FARARETEE  HAT RICES  FOR 

SHLD0G36 

c 

LEFT  PIGTAIL  SECTIONS 

S  HL  D  06  37 

c 

SHLD063E 

c 

SHLD0639 

~ 

DO  134  1=1, NS 

SHLDC640 

■ 

;i 

134 

LFT(I)  =IPL  (I) 

SHLD064  1 

-• 

► 

IPt  =  0 

SHLD0642 

r 

LHAXL=ZERO 

SHLD0643 

DC  136  1=1, TNT 

SIILD0644 

‘ 

DC  135  J=1 ,TNT 

8HLD0645 

135 

PHIL  (I, J) =ZEHOC 

SHLD0646 

136 

PHIl  (1,1)  =ONEC 

SHLD0647 

137 

CCNTINOE 

DO  138  1=1, NS 

SHLD0648 

SHLD0649 

138 

IF  (IPT  (I)  .GE.LNAXL)  IHAXL=LPT(I) 

S  HLD0650 

IF  (IEL.EQ.O)  LTL  =  LHAXL 

SHLD0651 

IE(LRAXL.EQ.  ZERO)  GO  TO  152 

SHLD0652 

IFL=IPL+ 1 

S  HLD0653 

» . 

l HI N 1= I HAXL 

SHLD0654 

DC  139  1=1, NS 

SHLD0655 

139 

IP  (LPT  (I) .LE.LRINL.AND.LPT(I)  .  NE.ZERC) 

l RINL= l PT (I) 

S  HLD0o56 

• 

DC  140  1=1, NS 

SI1LD0657 

*• 

KEY  (I)  =1 

SHLD0658 

140 

IF  (LPT(I)  .GE.LHINL)  KEY(I)=2 

SHLD0659 

CALI  INDOCT  (IND,NU,NS,ST,NUPS,NTYPE,HU, 

YU  ,ZH, R5,YS,ZS, 

SHLD0660 

1RPML,YPL,ZPL,RNH,KSY,NMJ,NNS) 

SHLD0661 

DC  141  1=1, NT 

SHLD0662 

DC  141  J=1 , NT 

SHLD0663 

A 

141 

Z  (I,J) =JOHEGA*IND(I,J)  *ZG 

SHLD0664 

w 

IF(NU.EQ.O)  GO  TO  143 

SHLD0665 

DC  142  1=1, NU 

SHLD0666 

142 

Z  (I,  I)  =Z  (1,1)  ♦ZUUV(I) 

SHI.DC667 

143 

CONTINUE 

SHLD0668 

CALL  INFADD(Z,KEY,NU,NS,NT,NOPS,ZSV,ZDV 

, LTV, ZHPLV,ZHUHV, NNU , NN  S) 

SHLD0669 

CALL  SCAP(SS,IND,NU,NS,NT,NUPS,KEY,NNU, 

NNS) 

SHLD067C 

CALL  ADRADD(KEY,NU,NS,NT,N0PS,RS,1S,RUH 

,EHS,SS,STV, M>U,NNS) 

SHLD067 1 

9 

125 


u  u  u  o  u  u 


DC  145  1=1, ST 
DO  144  J= 1 , NT 
T  (I,J)  =SS(I, J) *ONEC 

144  TI(I,J)=ZEROC 

145  T I  (1,1) =ONEC 

CALL  LEQT1C(T,NT,NT,TI,NT,NT,0,WA,IER) 

DC  146  1=1, NT 
DO  146  J=1 , NT 

Y  (I,  J) =JOMEGA*DREAL (TI  (I, J) ) 

146  Y1NS  (I,J)=SS  (I,J)/JC«EGA 
CALI  MULTC  (TI,Y,Z,NT,NT,NT) 

CALL  EIGCC(TI,NT,NT,2,GAM,T,NT,WK,IEB) 

DC  148  1=1, NT 
DC  147  J= 1 , N  T 
Z (I,J) =T (I, J) 

147  TI (I,J) =ZEROC 

148  TI  (1,1) =ONEC 

CALL  LEQT1C(Z,NT,ST,TI,NT,NT,0,WA,IER) 

DC  149  1=1, NT 

149  GAH (I) =CDSQRT  (GAM  (I) ) 

CALL  PHI(PHIT,NT,TNT,LMINL,GAH,Y,YINV,T,TI,EP,EN) 

CALI  MUITC (T PHI, PHIT, PHIL, TNT, TNT, TNT) 

DO  150  1=1, TNT 
DC  150  J=  1 , TNT 

150  PHI  L  (I  ,  J)  =TPHI  (I ,  J) 

DO  151  1=1, NS 

151  IF  (LPT  (I)  .NE. ZERO)  LPT (I) =LPT (I) - LMI NL 
I  FA  X  L=ZEFC 

GC  TO  137 

152  CONTINUE 

WRITE  (6,153)  IPL 

153  FORMAT  (//I  OX , ' THE  NUMBER  OF  LEFT  PIGTAIL  SECTIONS  =  • 

1,12//) 


CONFUTE  CHAIN  PARAMETER  MATRIX  FOB 
SHIELDED  SECTION 


ISHID=L-LTL-LTB 
IF(LSHLD.LE.ZEBO)  GO  TO  222 
DO  154  1=1, NS 

154  KEY  (I)  =1 

CALI  INDUCT (IND,NU, NS, NT, NUP5,NTYPE,RW,YW,ZW, RS,YS,ZS, 
1RPWfi,YPR,ZP8,RWH,KEY,NKU,NN5) 

DC  155  1=1, NT 
DO  155  J= 1 , NT 

155  Z (I, J) =JOHEGA*IND(I, J) *ZG 
IF(NU.EQ.O)  GC  TO  157 

DO  156  1=1, NU 

156  Z(I,I)  =Z(I,I)  *ZWWV(I) 

157  CONTINUE 

CALL  IHPADD(Z,KEY,NU,NS,NT,NUPS,ZSY,ZDV,LTV,ZWPLV,ZWNHV, NNU, NNS) 
CALL  SCAP(SS,IND,NU,NS,NT,NUPS,KEY,KNU,NNS) 

CALI  ADMADD(KEY,NU,NS , NT , NU PS , RS , TS , BWH,ERS, SS, ST  V , N NU, NNS) 

DO  159  1=1, NT 
DC  158  J=1,NT 
T  (I  ,  J)  =SS  (I,  J)  *0 NEC 

158  TI  (I,J) =ZSROC 

159  TI (1,1) =0N EC 

CALI  LECT1C (T,NT,NT,TI,NT,NT,0,WA,IER) 


SHLD0672 
SHLD0673 
SHLD0674 
SHLD0675 
SHLDC676 
SHLD0677 
5HLD0678 
SHLD0679 
SHLD0680 
SHLD0681 
SHLD0682 
SI1LD0683 
SHLD0684 
SHLD0685 
SHLD0686 
SHLD0687 
SHLDC688 
SHLD068S 
SHLD0690 
SHLDC691 
SHLD0692 
SHLD0693 
SHLDC694 
SHLD0695 
SHLD0696 
SHLDC697 
SHLD0698 
SHLD0699 
SHLD0700 
SHLD070 1 
SHLD0702 
SHLD0703 
SHLD0704 
SHLD0705 
SHLQ0706 
SHLD0707 
SHLD0708 
SHLD07  09 
SHLD0710 
SHLD071 1 
SHLD0712 
SHLD07  1  3 
SHLD07  14 
SHLD0715 
SHLD0716 
SHLD0717 
SHLD0718 
SHLD07 1 9 
SHLD0720 
SHLD07  2  1 
SHLD0722 
SHLD0723 
SHLD0724 
SHLD0725 
SHLD0726 
SHLD0727 
SHLD0728 
SKLD0729 
SHLD0730 
SHLD073 1 
SHLD0732 
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DO  160  1=1, NT 
DC  160  J  =  1 ,  N  T 

If  (I  ,J)  =JCMEGA*DREAL(TI  (I,J)  ) 

160  YIN  V  (I  ,  J) =SS  (I,J)/JCMEGA 
CALL  MULTC(TI,Y,Z,NT,NT,NT) 

CALI  E IGCC (XX, NT, NT, 2, GAM, T, NT, HK,IER) 

DC  162  1=1, ST 
DC  161  J  =  1 ,  N  T 
Z(I,J)=T(I,J) 

161  T1  (I  ,  J ) =ZEROC 

162  T I (I, I) =ON  EC 

CALL  LEQT1C (Z,NT,NT,TI,NI,NT,0,HA,IER) 

DC  163  1=1, ST 

163  GAM  (I)  =CDSQRT  (GAM(I)  ) 

CALL  PHI (TPHI,NT,TNT,LShLD,GAM,Y,YINV,T,TI,EP,EN) 
CALL  NULTC  (PHIT,TPiiI,  PHIL, T NT, T  NT, TNT) 

CALI  SULTC(TPHI,PHI8,PHIT,INT,TNT,TNT) 

GC  TC  173 


C  CCMFUTE  CHAIN  PARAMETEF  MATRICES  FOR 

C  LINE  HAVING  NO  SHIELDED  HIRES 

C 
C 

164  CONTINUE 

CALL  INEUCT(IND,NU,NS,NT,NUPS,NTYFE,RH, YH,ZW,RS,YS,2S, 
1RPHF  ,YFP,ZPR,RWH, KEY, NFU.NNS) 

DO  165  1=1, NT 
DC  165  J=  1 , N T 

165  Z  (I  ,J) =JGMEGA*IND(I,J)  *ZG 
DO  166  1=1, NU 

166  Z  (1,1)  =Z(I.I)  ♦ZHHV(I) 

CALI  SCAP(SS,IND,NU,NS,NT,NUPS,KEY,NNU,NNS) 

DO  168  1=1, NT 
DC  167  J=  1 , N  T 
T  <I,J)  =SS  (I,  J)  *0  NE 

167  TI  (I,J) =Z EROC 

168  TI  (1,1) =0N  EC 

CALI  LECT  1C(T,NT,NT,TI,NT,NT,0,HA,IER) 

DC  169  1=1, NT 
DC  169  J  =  1 ,  N  T 

Y  (I  ,J)  =JCMEGA*DREAL(TI (I,J)  ) 

169  YINV  (I,J)=SS  (I,J)/JCf1EGA 
CALI  MULTC(TI,Y,Z,NT,NT,NT) 

CALI  EIGCC(TI,NT,NT, 2 , GA M,T , NT , HK , I E R) 

DC  171  1=1, NT 
DO  170  J=  1 , N  T 
Z(I,J)  =T  (I ,  J) 

170  TI  (I ,J) =ZEBOC 

171  TI  (1,1) =ONEC 

CALL  LECT1C(Z,NT,NT,TI,NT,NT,0,HA,IER) 

DC  172  1=1, NT 

172  GAM  (I)  =CDSQR  I(GAM  (I)  ) 

CALL  PHI (TPHI , NT, TNT, L  ,G AM , Y , YIN V ,T , TI , E P , EN) 

173  CONTINUE 
C 

C 

C  REARRANGE  PHI  MATRIX 

C 

C 

C  REARRANGE  ROHS 
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n  n  n  n  n  no 


IF(NU.EC.O)  GO  TO  175 
DC  174  1=  1 , N  U 
DO  174  J=1,TNT 
PHI 1 <I,J) =TPHI  (I,J) 

174  PHIT  (I+NUPS, J)=TPHI  (I*NT,J) 

175  CONTINUE 
IF(NS.EC.O)  GO  TO  182 
DC  17b  1=1, NS 

DO  176  J=1,TNT 

PHIT (I ♦ NU , J) =TPHI  (I*NUPS, J) 

176  PHIT  (I*NUPS*NU,J)  =TPHI  ( I*NT *NUPS , J) 

K  =  0 

H  =  0 

DC  181  1=1, NS 

IF  (IGNDB  (I) . EC. 1)  GO  TC  178 
K  =  K  ♦  1 

DC  177  J=  1, TNT 

PHIT  (K*TNUPS,J)=TPHI(I*NU,J) 

177  PHIT  (K*TNUPS*NS*NISB, J) =TPHI  (I*TNUPS,J) 
GC  TO  180 

178  CONTINUE 
N  =  H*1 

DC  179  J=  1, TNT 

PHIT  (M  +  TNDPS  +  NVSR, J) =T PHI  (  I *T NUPS ,J) 

179  PHIT  (M*TNUPS*NS,J)=TPHI  (I*NU,J) 

180  CONTINUE 

181  CONTINUE 
132  CONTINUE 

DC  183  1=1, TNT 
DC  183  J  =  1 ,TNT 
183  TPH  I  (I  ,  J)  =PH  IT  (I ,  J  J 


BEARBANGE  COLUMNS 


IF(NU.EQ.O)  GO  TO  185 
DO  184  J  =  1 , N  U 
DC  184  1=1, TNT 
PHIT  (I, J) =TPHI  (I, J) 

184  PHIT (I, J*NUPS) =IPHI(I,J*NT) 

185  CONTINUE 

IF  (NS. EQ.O)  GO  TO  192 

DC  166  J=  1 , M  S 

DO  186  1=1, TNT 

PBTT  (I, J*NU)  =TPHI  (I, J* NUPS) 

186  PHIT (I, J*NUPS*NU) =TP H I ( I , J* NT* NUPS) 

K=0 

n=o 

DC  191  J=1,NS 

IF  (IGNDl  (J) . EQ. 1)  30  TC  188 
K  =K  ♦  1 

DC  187  1=1, TNT 

PHIT  (I , K*TNU PS) =TPHI  (I,J*NO) 

187  PHIT (I,K*TNUPS*NS*NISL)=TPHI  (I, J*TNUPS) 
GC  TO  190 

188  CONTINUE 
H  =  H  ♦  1 

DC  169  1=1, TNT 


SIILD  0794 
3HLD0795 
SIILD0796 
3HLDC797 
SHLD0798 
SIILDC799 
SHLD0800 
SHLD030 1 
SHLD0802 
SHLD0803 
3IILD0304 
SIILD0805 
SHLD08C6 
S  HLD0807 
SHLD0808 
SHLD0809 
SHLD03  10 
SHLD081 1 
SHLD08 12 
SHLD081 J 
SHLDC814 
SHLD0815 
SIILD08  16 
5HLD0817 
S  HLDOd  1 8 
SHLD0319 
SHLDC320 
SHLD082  1 
SHLD0822 
SHLD0823 
3HLD0824 
SHLD0325 
SHLD0826 
5 HLD0827 
SHLD0828 
SHLD0329 
SHLD0830 
SHLD0331 
SI1LD0832 
S  HLDOd  3  3 
SHLD08 34 
SULD03 35 
SHLD08  36 
SHLD08  37 
SHLD0J38 
SHLD0339 
SHLD0340 
SHLD094 1 
3HLD0d42 
SHLD0843 
SHLD0844 
SHLD0345 
SHLD0346 
3HLDC847 
SHLD0348 
S3LD0849 
SHLDC850 
SHLD0851 
SHLD0852 
SH  LD  03  33 
3HLD0864 


o  n  o  n  n 


PHIT (I,B*TNUPS*NVSL) =TPHI(I,J*IH0PS) 

3HLD0855 

189 

PHI 7  (I , fl+TNU  PS*NS) =TPHI  (I,J*HU) 

SHLDC856 

190 

CONTINUE 

SHLDG357 

191 

CONTINUE 

SHLD0358 

192 

CONTINUE 

SHLD0859 
SHLD0360 
SHLD036  1 

I NCC  SPOBATE  TEBHINAL  CONDITIONS 

SHLD0862 

SHLDC363 

SKLD0964 

DC  194  1*1, HUPS 

SHLD0865 

DC  194  J=  1 ,  N  OPS 

SHLD0366 

SOBC=ZEBOC 

SHLD0867 

SUHC 1=ZER0C 

SHLD0868 

SUBC2=ZEB0C 

SULD0869 

DC  193  K=1,NUPS 

SHLD0870 

SUHC=SOaC*PHIT  (I ,  K*NU  PS)  *Y0  (K,J) 

ShLD  0871 

suhc  i= subc i *phit (I*n ups, k*nups>  *yc  <k,j) 

SHLDC372 

193 

SUBC2=SUBC2*YL(I,K) *PHIT(K,J*NUPS) 

SHLD0373 

B  <I,J)  =-subc*phit  41,  J) 

SHLDC374 

C  (I,J) =  SUBC 1-PHIT (I*NUPS,J) 

SHLD0375 

194 

D  (I  ,J) =-SUBC2*PHIT(I*NUPS, J*NUPS) 

SHLD0876 

IF(BS.EC.O)  GO  10  199 

SHLD0377 

DO  196  1= 1, N  UPS 

SHLDC878 

DC  196  J=1,NS 

SHLD0379 

SOBC=ZEHOC 

SHLD0880 

DC  195  K=  1, N  UPS 

SHLD0881 

195 

suhc=subc*yl(I,k)  *phit {k,j*tnups) 

SHLD0882 

196 

A (I, J*NUPS) =SUBC-PHIT  (I*NUPS, J*TNUPS) 

SHLD0383 

DO  198  1=  1, NS 

SHLDC384 

DO  198  J-  1, HUPS 

SI1LD0885 

SOBC=ZEBOC 

SHLD0386 

DO  197  K=  1 , M  UPS 

SHLD0887 

197 

SUBC=SOBC*PHIT(I*TNUPS*NS,K*NUES) *10  (K,J) 

SHLD0388 

198 

A  |I*NUPS,J)=SUBC-PHI1  (I*TNUPS*NS,  J) 

SHLD0889 

199 

CONTINUE 

SHLDC890 

DC  202  1=1, HUPS 

SHLD089  1 

SUBC=ZEBOC 

SHLDC892 

DO  201  J=1,NUPS 

SHLD0893 

SUBC1=ZER0C 

S  HLD0894 

DO  200  K=1 , HUPS 

SHLDC895 

200 

S0BC1=SUHC1*YL(I,K)  *B  <K,  J) 

SHLDC896 

A  (I,J)  =C(I,J)  *30801 

SHLDC897 

201 

SUBC=SUBC*D (I,J) *10 (J) 

SHLD0898 

202 

EP(I)  =  IL(I)  *SUBC 

SHLDC399 

IF(BS.EQ.O)  GO  TO  206 

SHLD0900 

DC  203  1*1, NS 

SHLD 0901 

DO  203  J=1,NS 

SHLDC902 

203 

A  (I ♦ HUPS, J* HUPS) =-PHIT (I*TSUPS*HS, J+TNUPS) 

SHLD090  3 

DO  205  1*1, HS 

SHLD0904 

SUBC=Z  EBOC 

SHLDC905 

DC  204  J=  1 , HUPS 

SIILD0906 

204 

SUBC  =  SUBC*PHIT  (I *TNUPS *NS , J *NUPS ) *10 (J) 

SHLD0907 

205 

EP(I*NUPS)=SUBC 

SHLDC903 

206 

CONTINUE 

SHLD0909 

CALI  LECT1C(A,NT,NT,KP,1,NT,0,WA,IER) 

SHLD  09 1 0 

DO  208  1=1, HUPS 

SHLD091 1 

SU8C*ZER0C 

5HLD0912 

SOSC 1=ZEBOC 

ShLD 091 3 

DO  207  J* 1,M0PS 

SHLD0914 

S0BC=S0HC*B(I,J) *EP(J) 

SHLD09  15 
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207 

208 


209 

210 
211 


212 


213 

214 

215 


216 


217 

218 

219 

220 
221 


222 

223 

224 

225 


soil  i=sumci ♦phit  ( i ,  j ♦  n ') t :;)  *io  (j) 

EN (1)=SUMC*3UMC1 
I F (NS. EQ. 0)  00  TO  211 

DC  210  1=1 ,  N  U PS 
SUMC=Z  EROC 
DC  209  J  =  1 ,  N  S 

SUMC  =  SUMC*PHIT (I .J  +  TNUES)  *EP  (J'NUFS) 

EN ( I) =EN  (I) *SUMC 
CONTINUE 

IF(NU.EC.O)  GC  TO  215 
WRITE (6,212) 

FORMAT  (/////, IX. 'WIRE* ,14X, • VON (VOLTS)  •  ,15X, • VO  A ( 
128X, • VIM (VOLTS)  • ,15X,  • VL A  (DEGREES ) •//) 

DO  214  1=1, NU 
VWL=EP  (I) 

VWLP=CDA3S  ( V  WL) 

VWLR=DREAL  ( V  W L) 

VWLI  =  DIMAG(V«IL) 

VWLA=DATAN2 (VWLI  ,  V  W  LR)  *R ADEG 

V  WR  =  EN (I) 

VWRM=CDA3S  (VWR) 

VWR  R  =D  RE AL  (VWR) 

VWR  I  =  DXHAG (VWR) 

VWRA=DATAN2( VWRI ,VWRR)  *R ADEG 
WRITE  (6,213)  I,VWLM,VWLA,VWRM,VWRA 

FORMAT  (IX, 12,2  (15X, IP  £  10 . 3)  ,  1  5X , 2  (  15X , 1  PE  1 0. 3) ) 

CONTINUE 

CONTINUE 

IF(NS.EQ.O)  GC  TO  219 
WRITE(6,  216) 

FORMAT  (///, IX, 'SHI  ELDED  WIR E' ,5X , • VOM  (VOLTS)  * ,15X 
128X,«VLM  (VOLTS)  * ,  15X  ,  ' VLA  (DEGREES) *//) 

DC  218  1=1, NS 
VSWI=EP (I +NU) 

VSWLM=CDABS (VSWL) 

VSWLH=DREAL( VSWL) 

VSWII=DI MAG (VSWL) 

V  SW  L  A=  D AT AN  2 (VSWL I, VSW IR) *RADEG 
V SW  B=EN (I+NU) 

VSWB  M=CDABS (VSWR) 

VSW  PR=DREAL (VSWB) 

VSWBI=DIMAG(VSWR) 

VSWRA=DATAN2 (VSWBI, VSWBB)  *R ADEG 

WRITE  (6,217)  I, VSW LM, VSWL A, VSWB M, VSW BA 

FCRMAT(1X,I2,2(15X, 1PE10.3)  ,  1 5X , 2  ( 15 X  ,1  PE  1 0. 3) ) 

CONTINUE 

CONTINUE 

GC  TO  104 

CONTINUE 

WRITE (6,221) 

FORMAT(20X,' INPUT  DATA  ERROR'/ 

120X, 'CONSULT  USERS  MANUAL  FOR  ALLOWEE  RANGE  CP  FO 
GO  TO  224 
CONTINUE 
WRITE (6,223) 

FORMAT (10X, ' THERE  IS  A  ZERO  OR  NEGATIVE  LENGTH  SH 
ID  WIBE'/IOX, 'CHECK  PIGTAIL  LENGTHS  AGAINST  IOTAL 
CCNTINUE 
STOP 
END 


DEGREES) 


SHLDO?  16 
SIILDC7  17 
SHLDO  3 1 6 
3I1LD0J  19 
SH  LD  092  0 
SHLD0921 
SHLD0922 
SHLDC923 
SHLD092U 
SHLD0925 
SHLDC926 
SIILD0927 
SHLD0928 
SHLDC929 
3  IiL  DO  9  30 
SHLDO  3  3 1 
SH  LD  0932 
SHLD0933 
SHLD0934 
SHLD0935 
S1ILD0936 
SHLD0937 
SHLD0930 
SHLD0939 
SHLD0940 
SHLD094 1 
SHLD0942 
SHLD0943 
SHLD0944 
SHLD0945 

,' VOA  (DEGREES) • ,  SHLD0946 
SHLDC947 
SHLD0948 
SHLDO 949 
SHLDC95C 
SHLD095  1 
SHLD0952 
SHLD0953 
SHLD0954 
SHLD0955 
SHLDC956 
SHLD0957 
SHLD0956 
SHLDC959 
SHLD0960 
SHLD0961 
SHLD0962 
SHLD0963 
SHLD0964 
SHLDC965 
SHLD0966 
SHLD0967 
SHLDC968 
SHLD0969 
SHLD0970 
ON  A  SHIELDESHLD0971 
LENGTH')  SHLD0972 
SHLD0973 
SHLD0974 
SHLD0975 
FLS  10001 


PM  AT  ) 


I  ELD 
LINE 
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C 

c 

C 

C 

c 

c 


FUNCTION  LSI  FOB  COMPUTING  THE  SELF  INDUCTANCES  CF  CCHDUCTOES 
ABOVE  GBCUND 


BEAI  FUNCTION  LS1*8(RW,H) 

IMPLICIT  REAL*8  (A-H.O-Z) 

REA  I *3  HU,MU02PI,HU04PI 

CCMBCN  /RCON/  ZE RO, ONE  ,TRC, THREE, TO U R, MU, MD02PI, M U04PI, 
1PI,BADEG, ERT  E , S IGCOP ,R ,  F  ,  V2 
LS1=BU02PI*DLCG(TH0*H/BN) 

RETURN 

END 


FUNCTION  LM1  FOR  COMPUTING  MUTUAL  INDUCTANCES  OF  CONDUCTORS 
ABOVE  GROUND 


REAL  FUNCTION  LM 1 *3 ( II  ,Y J,Z I, ZJ) 

IMPLICIT  R E AL*8  (A-H.O-Z) 

REAl*8  MU,NU02PI,MU04PI 

CCMBCN  /RCON/  ZERO,  ONE ,THC, THBEE, FOUR, HU, BU0  2 PI , MU04PI , 
1PI, RADEG,ERTE,S IGCOP, R,F,V2 
ZC=ZI-ZJ 
YD=YI-YJ 

DIJ2=ZC*ZD»Y0*YD 

LM1=NUO4PI*DLCG(ONE«T0UR*YI*YJ/DIJ2) 

BETURN 

END 


FUNCTION  LS2  FCR  COMPUTING  THE  SELF  INDUCTANCES  CF  CONDUCTORS 
IN  AN  CVERALL  SHIELD 


REAL  FUNCTION  LS2*8  (BHI,RI) 

IMPLICIT  REAL*8  (A-H,0-Z) 

BEAl*8  MU,HU02PI,MU04PI 

CCHBCN  /RCON/  ZERO .ONE, THC,THBEE, FOUR, MU, BU02PI, M U04PI, 
1PI.RACEG, ERTE, S IGCOP  ,  B  ,  P , V2 
LS2=HU02PI*DL0G( (B*8-RI*RI) /(R*RNI)J 
RETURN 
END 


FUNCTION  LH2  FOR  COMPUTING  THE  MUTUAL  INDUCTANCES  OF  CONDUCTORS 
IN  AN  OVERALL  SHIELD 


REAL  FUNCTION  LM 2*8 (R I ,R J, THI , TH J ) 

IMPLICIT  BE AL*8  (A-H.O-Z) 

RE Al*8  MU,HU02PI,MU04PI 

COMMON  /RCON/  ZERO , ONE , TWO , THREE , FOUR , MU, BU02 PI , MUC4 PI , 
1PI.BADEG, ERTE, S IGCOP, R  ,  F , V2 
RI2=RI*RI 
RJ2=BJ*RJ 
THETA=THI-THJ 
RS2=B*R 

LM2=HU04PI*DL0G( (RJ2/RS2) « (RI 2*R J 2* R S2*RS2-T«C* RI * RJ *FS 2  * 
1DCOS  (THETA) ) / { 812 *R J 2» B J2*R J2-T NO *RI *RJ *RJ 2* CCOS  (THETA) ) ) 


FLS10002 
ELS  10003 
FLS  10004 
FLS10005 
FLS 10006 
FLS 10007 
FLS 1 0008 
FLS 1 0009 
FLS  10010 
FLS 10011 
FLS  10012 
FLS  1 00 1 3 
FLS 10014 
FLM10001 
FLM10002 
FLM10003 
PLH  10004 
FLM10005 
FLM 10006 
FLM10007 
FLM 100C8 
FLM  10009 
FLM  100  10 
FLM  1 00 1 1 
FLM  10012 
FLM  100  1 3 
FLM  10014 
FLM  100  15 
FLM  100  16 
FLM 1 00 1 7 
FLS2000  1 
FLS20002 
FLS  20003 
FLS  20004 
FLS2&005 
FLS  20006 
FLS20007 
FLS20008 
FLS 20009 
FLS200  10 
FLS2001 1 
FLS 20012 
FLS20013 
FLS20014 
FLM  2000  1 
FLM20002 
FLM20003 
FLM  20004 
FLM  20005 
FLM2C006 
FLH20007 
FLM20008 
FLM2C009 
FLS200  10 
FLM200  11 
FLB20012 
FLM  200  1 3 
FLM200  14 
FLM20015 
FLM  200  1 6 
FLM20017 
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RETURN 

END 


FUNCTION  STB  FOR  COMPUTING  THE  TRANS FEB  ELASTANCES  FCP 
BRAIDED  SHIELDS 


REA l  FUNCTION  STB*8 ( BS , R E , TUN , E , NP E , E RS) 

IMPLICIT  RE AL*8  (A-H,0-Z) 

INTEGER  E,MP3 
REA L *4  TABLH  (5,9) 

REAl*S  NU,nU02PI,MU04PI 

CCMflCN  /RCON/  ZERO,  ONE,  1HO,  THREE,  FOUR,  NU,  NU02EI , MU04 PI, 
1PI,PADEG,ERTE, S IGCOP , R , F , V2 

DATA  TABLH/.  521, 1,0 36, 1.  545, 2. 049, 2. 54 7,. 502,  1.001,  1.496, 

1  1.96  7,  2.4  75,  .491,. 987,  1.463, 1.  944, 2. 423,. 4 80,  .956,1.431  , 
21.902,2.  371,.  471,.  939,  1.  404 , 1 . 866 ,2. 326 ,. 463, .922, 1 . 379, 

3 1. 834, 2. 2 85,. 4  58,. 9  12,  1.  364, 1.  81  3,2.259, .455, .9  06, 1  .  355, 

4  1.80  1,  2. 244,. 454,. 904,  1.  352,  1.  796,2.238/ 

THW=THH*RADEG 

XF=FLO AT (3) *PL0A1  (4PB) ♦THO*RB/ (FO UR*EI*RS*DCOS {I H4/R ADEG)  ) 
CPC=TNC*XF-XF»XE 

F AC  =  (  (CNE-OPC) **  (THREE/TNO)  ) *TWC* E FS/ (ERS »ON E) 

DX= 1.0D-2 
D Y=  5 . 0  DO 

IF  (FAC.LE.DX)  GO  TO  2 
IF(FAC.CE. 5.0D-2)  GO  TO  3 
DC  1  1=1,5 
FAC1=DX*PLCAT (I) 

X  =  F  AC-  F  AC  1 

IF ( X .UE. ZERO)  GO  TO  1 

IX 1 =1-  1 

1X2=1 

DIFX  =  FAC-DX*FLOAT  (1X1) 

GC  TC  4 

1  CONTINUE 

2  1X1=1 
1X2  =  2 
DIFX  =  Z  ERO 
GC  TC  4 

3  1X1=4 
1X2=5 
DIF  X  =  D  X 

4  CONTINUE 

IF(THW.LE.DY)  GO  TO  6 
IF(THH.GE.45.0D0)  GO  TC  7 
DC  5  J=1 ,9 
PSI1=DY*FLOAT(J) 

Y=TH4-PSI 1 

IF  (Y.GE.ZERO)  GC  TO  5 

JI1=J-1 

JY2  =  J 

DIFY=THN-DY* FLOAT  (JY1) 

GO  TO  8 

5  CONTINUE 

6  J  Y 1  =  1 
JY2=2- 
DIF  Y  =  ZERO 
GC  TO  8 

7  JY1=8 


FLM2C018 
FLM20019 
F3TE000 1 
FSTB0002 
FSTD0003 
FSTD0004 
FSTB0005 
FSTB0006 
PSTB0007 
FSTB0008 
FSTB00C9 
FSTBOO  10 
FSTB001 1 
FSTE0012 
FSTBOO  13 
FSTB001 4 
FSTB0015 
FSTBOO  16 
FSTBOO 1 7 
FSTBOO 18 
FSTB0019 
FSTB0020 
FSTBOO  2  1 
FSTB0022 
FSTB0023 
FSTB0024 
FSTB  002  5 
FSTB0026 
FSTB0027 
FSTB0028 
FSTB0029 
FSTB0030 
FSTB003 1 
FSTB0032 
FSTBOO 33 
FSTB0034 
FSTB0035 
FSTBOO  36 
FSTB0037 
FSTB0038 
FSTBOO  39 
FSTB0040 
FSTB004 1 
FSTBOO  42 
FSTB004  3 
FSTB0044 
FSTB0045 
FSTB0046 
FSTEC047 
FSTB0048 
FSTB004  9 
FSTBC050 
FSTB005 1 
F5TB0052 
FSTBC053 
FSTB0054 
FSTB  0055 
FSTE0056 
FSTB0057 
FSTB0058 
FST30059 
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c 

c 

c 

c 

c 

c 


JY2  =  9 
DIPY  =  D¥ 

8  CONTINUE 

P1=TABLH  (I  X 1  ,  J  T  1 ) 

F2=TABLH  (1X2, JY1) 

F3=TABLH(IX1,JY2) 

P4  =  TABIH  (1X2, JY2) 

H=F 1 *DIFX/DX*  (F2-F1) ♦DIFY/Dt*  (F3-F1) ♦DIPY*DIFX/ (D¥*BX) » (F4- 
1F3-F2  +  F1) 

H  =  H  *DX 

SIB  =  H/  (THO*FLOAT (3) *EHS*ERTE) 

T HH=THH/R ADEG 

RETURN 

END 


FUNCTION  LTB  FOR  COMPUTING  THE  TRANSFER  INDUCTANCES  FOR 
ERAICEE  SHIELDS 


REAL  FUNCTION  LTB*8 (RS , R B, TH8 , B, HPB) 

IMPLICIT  REAL*8  (A-H,0-Z) 

INTEGER  B, HP B 
R  EAL*4  TABLG (5, 9) 

REA  I *8  HU,HU02PI,aU04PI 

CCHMCN  /BCON/  Z ERG, ONE, THO, THREE, FOUR, MU, N002PI,M DO 4PI, 
1PI,8ADEG, ERTE, S IGCOP ,R , F , V2 

DATA  TABLG/. 539,  1. 09 1 ,  1. 657 , 2. 236 , 2. 829, . 555 ,  1 .  1 1 7,  1 . 6 96 , 

12. 26 1.2. 844..  584.1. 173.1.767.2. 366.2.569. .62 1 . 1. 246. 1 . 874. 2. 507, 

23. 143. . 667. 1.337.2. 010.2.687.3.366. .727.1.457.2. 190.2.928, 

33. 66e,. 807, 1.6 18, 2. 433,3.253,4.076, .910,1.825,2.745,3.671, 
44.603, 1.045, 2. C98, 3. 156,4.226,5. 303/ 

THH=THH*RADEG 

X  F=  F  LO AT (3)  *  FLOAT (HPB) *THO* RB/ (FOUR*PI*RS*DCOS (TH W/P ADEG) ) 

OFC=THO*XF-XP*XF 

PAC  =  (CNE-OPC) ♦*  (THREE/THC) 

DX= 1.0D-2 
DY=5.0D0 

IF(FAC.IE.DX)  GO  TO  2 
IF(FAC.GE.5.0D-2)  GO  TO  3 
DC  1  1=1,5 
FAC1=DX*FL0AT (I) 

X=F  AC- F  AC  1 

I F  (  X .G  E. ZERO)  GO  TO  1 

IX1=I-1 

1X2=1 

DIFX=FAC-DX*FLOAT(IX1) 

GC  TC  4 

1  CONTINUE 

2  1X1=1 
1X2  =  2 
DIFX=Z  ERO 
GC  TC  4 

3  1X1=4 
1X2  =  5 
DIF  X  =  D  X 

4  CONTINUE 

IF(THH.LE.DY)  GO  TO  6 
IF(THH.GE.45.0C0)  GO  TC  7 
DC  5  J=1 ,9 
PSI1=DY*FLOAT(J) 


FSTB006J 
FSTBC06 1 
FSTB0O6  2 
FSTBOOfi  3 
FSTB0064 
FSTB0065 
FSTB0066 
FSTB0067 
F3TB0068 
FSTB0065 
FSTB0970 
FSTE00  7  1 
FSTB0072 
FSTD007  3 
FLTB000 * 
FLTB0002 
FLTB000  3 
FLTB0004 
FLTB0005 
FLTB0006 
FLTE0007 
FLTB0008 
FLTB0009 
FLTBOO  10 
FLTBOO  1  1 
FLTBOO 1 2 
FLTEOO  1  i 
FLT300  14 
FLTBOO  1  5 
FLTEOO  1b 
FLTBOO 1 7 
FLTBOO  18 
FLTBOO 1 9 
FLTS0020 
FLTBOO/  1 
FLTB0022 
FLTBC333 
FLTB0024 
FLTB0025 
FLTBC026 
FLTB0027 
FLTB0028 
FLTBC029 
FLTE0030 
FLTBOO  3 1 
FLTE0032 
FLTBOO  33 
FLTBOO  34 
FLTE0035 
FLTDOO  3fc 
FLTB0037 
FLTCOO  38 
FLT  D  00  35 
FLTB0040 
F  LI E  004  1 
FLT BOO  4 2 
FLTB0043 
FLTB0044 
FLTB0045 
FLTB0046 
FLTEC047 


•  ' 
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c 

c 

c 

c 

c 

c 


Y=THW-PSI 1 

IF  (Y.GE.ZERC)  GC  TO  5 

J  Y 1  = J-  1 

JY2=J 

DIFY=THW-DY*FLCAT  (JY1) 

GC  TO  8 

5  CONTINUE 

6  J  Y 1  =  1 
JY2=2 
DIFY=ZERO 
GC  TC  8 

7  JY  1  =  8 
J  Y2  =  9 

D IF Y  =D Y 

8  CONTINUE 
F1=TABLG  (1X1 ,JY1) 

F2=TABLG(IX2, JY1) 

F 3=T AB LG  (IX  1,JY2) 

F4  =  T  ABLG  (1X2, JY2) 

XG=F1*DIFX/DX*(F2-F1)  ♦DIFY/DY*  (F3-F1)  *DIFY*DIFX/(DX*DY)  *(f*-F3- 
1F2«-F  1) 

XG=  XG*  D  X 

LTB=MU*XG/(TWO*FLOAT(B) ) 

T H  W  =  TH  W/RADEG 

RETUBN 

FND 


FUNCTION  ZWW  FOR  COMPUTING  THE  SEIF  IMPEDANCES  OF  STRANDED  WIRES 


COMPLEX  FUNCTION  ZWW* 16 <R ST , NST, SIG) 

IFPLICIT  REAL*8  (A-H.O-Z) 

REAl*8  MU,NU02PI,MU04PI,LDC 
CCM  FLEX  *  1 6  ZEROC,ONEC,XJ,JOMEGA 

DATA  P25/.25DO/.ONEP15/1. 15D0/, P5/.5C0/, PI  5/. 15  DO/ 

COMMON  /RCON/  ZERO , ON E ,T KC , THREE, FOUR , MO, MUO 2PI , M UO 4P I , 

IP I, BA DEG, ERTE, S IGCOP ,K  ,  E , V2 
CCMrCN  /CCON/  ZEROC,ONEC,XJ, JOBEGA 
DELTA  =  ONE/(THO*PI*CS08T (SIG*SIGC0P*F*MU04PI)  ) 
RCC=ONE/(PI*SIG*SIGCOP«RSI*RST) 

LDC  =  FUC4  PI/T WC 
IF(RST.LE. DELTA)  GO  TO  1 
IF(RST.GE.THRE E* DELTA)  GO  TO  2 

ZWW  =  (P25*  (RST/OELTA»THREE) ♦RDCtJOEEGA* (ONEP 1 5-P 1 5 *R ST/DELTA )  ♦ 
1 LEC) /FLOAT  (NST) 

GC  TO  3 

1  ZWW  =  (ROC* JCMEGA*LDC) /FLOAT  (NST) 

GO  TO  3 

2  ZNN  = (P5*BST*HCC/DELTA*JOMEGA*TWO*DEITA*lDC/RST) /FLOAT  (NS f ) 

3  CONTINUE 
RETURN 
END 


FUNCTION  ZD3  FOR  COMPUTING  THE  DIFFUSION  IMPEDANCES  OF  BRAIDED 
SHI E IDS 


COMPLEX  FUNCTION  ZDB* 1 6  (BS , RB, S IG E ,T HU , B, W PB) 
IMPLICIT  REAL*8  (A-H.O-Z) 


FLTZ004U 
FI7B0046 
FLT3C05C 
FLT3005  1 
FLTB0052 
FL:bC053 
FLT30054 
FLTUC355 
FLTB0056 
FLT30057 
FLTU0058 
FLTE0059 
FLTBC060 
FLTB006  1 
FLTB0J62 
FLT300G3 
FLT30064 
FLTE0065 
FLTB0066 
FLTB0067 
FLTB0066 
FLTB0C69 
FLTB007C 
FLTB007  1 
FLTB0072 
FLIBC073 
FZWW000  1 
FZWWC002 
FZWW0003 
FZWW0004 
FZWWC005 
FZ W  W  0006 
FZKW0007 
FZWW0008 
FZWW0009 
FZWWOO  10 
FZKWC01  1 
FZWWOO  12 
FZWWOO  13 
FZWW0014 
FZWWOO  15 
FZWWOO  16 
FZWWC017 
FZWWOO 18 
FZWWOO 19 
FZWWC02C 
FZWW002  1 
FZWV0022 
FZWWC023 
FZW  WOO  24 
FZWW0025 
FZWW002b 
FZWW0027 
FZDBOOO  1 
FZDB0002 
FZDBOOO 3 
FZ0B0004 
FZDB0003 
FZDB0006 
FZDB0007 
FZ3B0008 
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c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 


INTEGER  3, WPS 
R  EA l *8  MU,  "UC2PI ,  -1')04ri 

COMPLEX* 16  ZEROC,ONEC,XJ,JOMEGA,GA:i,SINH 

CCMMCN  /RCCN/  ZERO,  ON  E  ,TWC,  THREE  ,  F  C  OR  ,  KU,  MIJ02  P I  ,  X  004  M  , 
1PI,PADEG,E:>TE,SIGCC?,R,F,V2 
COMMON  /CCON/  Z ErtOC , OH  EC , X J , JOM EG  A 
DELTA=CNE/DSQRT{PI*F*MU*S1GB*SIGCCP) 

GAM  = (ONE*XJ) /DELTA 
D=TWC*  RB 

SINH-(CDEXP(GAM*C) -CDEXF(-GAM*D) ) /TWO 

RDC  =  CNE/ (PI*83*RE*3IGB*SIGCCP*PL0AT  (  E) *  I LO AT  (  WPD)  * 

1DC0S  (THH)  ) 

ZCE=RDC*GAM*D/SINH 

RETORN 

FND 


FUNCTION  ZSB  FOR  COMPUTING  THE  SELF  IMPEDANCES  CF  BRAIDED  SHIELDS 


COMPLEX  FUNCIION  ZSB* 1 6 (RS, R B , SI GE , T HS ,B , HPB) 

IMPLICIT  RE AL*8  (A-H,0-Z) 

INTEGER  B , HPB 

RE AL*8  MU,HU02PI,M004PI,LDC 
COMPLEX* 1 6  ZEROC, ONEC,XJ,JONEGA 

CCMMCN  /RCCN/  ZERO, ONE ,T WC,THRE E, f CU R, MU, HU02PI , MU04PI , 
1PI/RACEG, ERT  E, 3 IGCOP , R , F, V2 
CCMMCN  /CCON/  ZEROC, CN EC , X J , JOMEG A 
DATA  P25/.25D0/,CNEP15/1. 15 DO/ , P5/ . 5D0/, P 1 5/.  1500/ 
DFLTA=CN2/(TH0*PI*CSQRT (SIGB*SIGCCP*F*8004PT)  ) 
RDC=CNE/(PI*SIGE*3IGC0P*RB*RB) 

IDC=BUC4PI/TNC 

IE(RE. LE. DELTA)  GO  TO  1 

IF(RE.GE.THREE*DELTA)  GC  TO  2 

ZSB=  (P25* (R3/DELTA*THHEE) *H  CC*0  ./N  EGA*  (ONEP  1 5-PI  5*  RB/DELT  A)  * 
1LEC)/(FL0AT(B) * FLO AT  (HPB) *DCOS (THH) ) 

GC  TC  3 

1  ZSB=  (RDC*JCMEGA*LDC)/ (FLCAT (B) *FLCAT (HPB) *DCOS(THW)  ) 

GC  TO  3 

2  ZSB= (P5*RB*RDC/DELTA*JCHEGA*THO*DEITA*LDC/R3)  /(FLOAT ( D)  * 

IF  10 AT  (HPB) *DCCS (THW)  ) 

3  CONTINUE 
RETURN 
END 


PUNCTICN  ZDS  FCR  COMPUTING  THE  DIIFUSION  IMPEDANCES  OF 
SOLID  SHIELDS 


COMPLEX  PUNCTION  ZDS* 1 6 (RS, TS, SIG S) 

I  BP l IC IT  REAL*3  (A-H,0-Z) 

RE A l *8  MU,MU02PI,MU04PI 

COMPLEX* 16  ZSROC,ONEC,XJ,JOMEGA,GAM,SINH 

CCMMCN  /RCCN/  ZERO  ,ONE  ,THC,THRE E, ECU E , HO, MUO 2P I , M U04P I , 
1PI, RADEG, ERT  E, S IGCOP ,R , F , V2 
COMMON  /CCON/  ZEROC ,0N EC ,XJ , JOMEG A 
DELT A=CNB/DSQRT ( PI *F*HO *S 1GS *S IGCOP) 

GAH=  (ONE*XJ) /DELTA 

R CC  =  ON  E/ (THO*PI*SIGS*SIGCOP*TS* (RS-TS/THO) ) 

SI NH  = (CDEXP (GAB*TS) -CD EX P  (-G AM*TS) ) /THO 


FZDB0003 
FZDBOO 1 0 
FZ  DECO  1  1 
FZDbOO  12 
FZDB001 3 
FZ  DB  00 1  4 
FZDBOO  15 
F2D3001 6 
PZDBC01 7 
FZDBOO  18 
FZDBOO 1 9 
FZDCC020 
FZDB002  1 
FZD30022 
FZDC0023 
FZSBOOO  1 
FZSL0002 
FZSE0003 
FZSB0004 
FZSB0005 
FZSE0006 
FZSBO  07 
FZSB0008 
FZSBC009 
FZSB0010 
FZSB001 1 
FZSBC012 
F7SB00 13 
FZSBOO 1 4 
FZ3E0015 
FZSB0016 
FZSB001 7 
FZSBC018 
FZSBOO 19 
FZSB0020 
FZSB0021 
FZSB0022 
FZS30023 
FZ3B0024 
FZSB0025 
FZSB0026 
FZ330027 
FZSB0028 
FZSB0029 
FZDS0001 
FZDS0002 
FZDS0003 
FZDS0004 
FZDS0005 
FZDS0006 
FZDS0007 
PZDS0008 
FZDS0009 
FZDSOO  10 
FZDSOO  1  1 
FZD50012 
FZDSOO 13 
FZDSOO  14 
FZDSO  0 1 5 
FZDSOO  16 
FZDSOO  17 


*1 


m 
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ZDS=RDC*GAH*TS/SINH 

RETURN 

END 

C 

c 

C  FUNCTION  ZSS  FCB  COMPUTING  THE  SEIF  IHPEDANCES  OF 

C  SOLID  SHIELDS 

C 

c 

CCH F IE X  FUNCTION  ZSS*  1  6  (RS,TS,SIGS) 

IHPLICIT  BEAL*8  (A-H,0-Z) 

BEA I  *8  HU,HU02PI,NU04PI 
CCHFIEX*16  ZERCC,ONEC,XJ,JOHEGA 

CCHHON  /BOON/  ZERO, ONE, T»C, THREE, FO UB, MU , RUC2 PI , HU04PI, 
1PI,PADEG,ERTE,SIGCOP,R,F,V2 
CCMFCN  /CCON/  ZEROC,ONEC,XJ, JOMEGA 
DELTA=ONE/DSQRT(PI»F*MU*SIGS*SIGCOP) 

RSS=CNE/ (SIG3*SIGCOP*DElT A) 

CCSH=  (DEXP  (THO*TS/DELTA) ♦CEXP (~THO*TS/ DELTA)  )  /TWO 
SINE=(DEXP(THO*TS/DELTA)-DEXP  (- TiO*T S/ DELTA) )/THO 
RB=BSS* ( (3INH*DSIN(THO*TS/ DELTA) ) /  (COSH-DCOS  (THO*TS/ DELTA) ) )/ 

1  (THC*PI*  (BS-TS/THO) ) 

OLI=BSS* ( (SINH-DSIN (THO*TS/ DELTA) )/ (COSH-DCOS (T HO =TS/ DELTA) ) )/ 

1 (THC^PI* (BS-TS/THO) ) 

ZSS=RB»XJ*CLI 

PETURN 

END 

C 

C 

C  SUBROUTINE  HULTC  FOB  MULTIPLYING  THO  COHPLEX  1 ATR ICES 

C 

C 

SUBPOUTINE  NULTC  ( A  ,  B , C , NL, NH, NB) 

IHPLICIT  BEAL*8  (A-H.O-Z) 

C0HFIEX*16  A  (NL.NR)  ,  B(NL,NH)  ,C(NH,NB)  , SUMC, ZEROC, ON  EC ,XJ, JOMEGA 

CCHHCN  /CCON/ZEROC,ONEC,XJ, JOHEGA 

DC  2  1=1, NL 

DC  2  J=1 , NB 

SUHC=ZEBOC 

DC  1  K=  1, NH 

1  SUHC=SUNC*B  (I, K) *C (K,J) 

2  A  (I, J) =SUMC 
RETURN 

END 

C 

C 

C  SU3BCUTINE  SOAP  FOR  CONFUTING  THE  INVERSE  O'  THE  CAPACITANCE 

C  NATRIX 

C 
C 

SUBROUTINE  SCAP ( SS , IND , NU ,N S, NT ,N OPS  , BE Y , N NU , NNS) 

IFPIICIT  REAL*8  (A-H.O-Z) 

INTEGER  KEY  (NNS) 

REAL*8  IND(NT.NT) ,SS(NT,NT) 

V=2.99792SD8 
V2=V^V 
DO  1  1=1, NT 
DC  1  J= 1 , NT 
1  SS (I, J ) =V2* I ND  (I , J) 

I F ( NS. Eg. 0)  GO  TO  4 
DC  3  1=1, NS 


FZDS001  8 
FZDSC019 
FZDS0020 
FZSS000 1 
FZSS0002 
FZSSOOO  3 
FZSS0904 
FZSS0005 
FZSS0006 
FZSS0007 
FZSS0008 
FZSSOOOS 
FZSSOO 1 0 
FZSS001 1 
FZSSOO 1 2 
FZSS001 3 
FZSSC014 
FZSSOO  15 
FZSSO0 16 
FZSS0017 
FZSSOO  18 
FZSSOO 1 9 
FZS30020 
FZSSOO  2 1 
FZSS0022 
FZSS0023 
FZSS0024 
SHUL0001 
SHUL0002 
SHUL0003 
SHUL0004 
SHUL0005 
SNUL0006 
SHUL0007 
SHUL0008 
SHUL00C9 
SHULOO  10 
SNULOO  1 1 
SHUL0012 
SHU  LOO  1 3 
SHULOO  14 
SHUL0015 
SHULOO 16 
SHULOO 17 
SCAP0001 
SCAP0002 
SCAP0003 
SCAP0004 
SCAF0005 
SCAP0006 
SCAP0007 
SCAP0008 
SCAP0009 
SCAPOO  10 
SCAPOO  1 1 
SCAP0012 
SCAPOO  13 
SCAPOO  1 4 
SCAPO0 1 5 
SCAPOO  16 
SCAPOO  17 
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IF  (RET  (I) .EQ. 2)  GO  TO  2 

SCAPOO 1 8 

SS(I*NOPS,I*NUPS)=SS (I»H0,I»MU) 

SCAP001 9 

2 

CONTINUE 

3CAP0020 

3 

CONTINUE 

SCAP002 1 

4 

CONTINUE 

SCAP0022 

RETOBN 

SCAP00  2  3 

END 

SCAP0024 
SINDOOO 1 
SIND0002 

SUBBCUTINE  INDOCT  FOR  CCHEUTING  THE  INDUCTANCE  MATRIX 

SIND0003 

SIND0304 

SIND0005 

SUBROUTINE  I NDUCT (I ND , NU , NS, NT , NU PS, NTTPE, RN , YH ,Z H, BS,YS,ZS, 

SIND0006 

1RHP, YP,ZP,BHH,KEY,NNU,NNS) 

SIND0007 

IMPLICIT  REAL*8  (A-H.0-7.) 

SIND0008 

INTEGER  KET(NNS) 

SIND0009 

RFAL*8  HU,HU02PI,MU04PI 

SIND0010 

REAI*8  IND (NT, NT) ,RH(NNU)  ,  YH  (NNU)  ,RS  (NNS)  ,YS(NNS)  ,RHH  (NNS)  , 

SIND0011 

IZS(NNS)  ,  YP  (NNS)  ,ZH  (NNU)  ,ZP(NNS)  ,  RHP  (NNS)  ,  LS  1,  LM  1,  LS2,  Lfl2 

SIND0012 

COMMON  /RCOH/  ZERO , ONE ,THO, THREE , FOUR, HO , HU02PI , HU04 PI, 

SINDOO 13 

1  PI,BADEG,ERTE,SIGC0P,R,F,V2 

SIND0014 

GO  TO  (1,29),  NTYPE 

SINDO0 1 5 

1 

CONTINUE 

SINDOO  16 

IF  (NU.  EC. 0)  GO  TO  3 

SINDOO  17 

DC  2  1=1, NU 

SIND0018 

2 

IND  (J,  I)=LS1  (RH(I)  ,  YH  ( I)  ) 

SINDOO  19 

3 

CONTINUE 

SINDOO  20 

IF  (NS.  EQ.  0)  GO  TO  8 

SIND  002 1 

DC  7  1=1, NS 

SIND0022 

I  RE  Y  =  KEY  (I) 

SIND0023 

GO  TO  (4,5) ,  IKEY 

SIND0024 

4 

IND  (I* NO, I*NO) =LS1 (R S ( I)  , YS  (I) ) 

SIND0025 

GC  TC  6 

SIND0026 

5 

IND  (ItNU, I»N  U) =  LS  1  (RHP  (I) , YP  (I) ) 

SIND0027 

6 

CONTINUE 

SIND0028 

IND  (I»NDPS,I  +  NUPS)=LS1  (RHH (I) , Y S ( I) ) 

SIND0029 

7 

CONTINUE 

SIND0030 

8 

CONTINUE 

5IND003 1 

IF  (NU. EC. 0)  GO  TO  10 

SINDOO  32 

NUH 1=NU— 1 

SIND0033 

DO  9  1=1, NUH  1 

SINDOO  34 

IP1=I*1 

SIND0035 

DC  9  J=IP1,NU 

SIND0036 

9 

IND  (I,  J)  =LM  1  ( YH  (I)  ,  YH(J)  ,ZH(I)  ,  ZH(J)) 

SINDOO  37 

10 

CONTINUE 

SIND0038 

IF (NS. EQ . 0)  GO  TO  21 

SIND0039 

NSH 1=N  S-  1 

SIND0040 

DO  16  I=1,NSM1 

SIND004  1 

IP1=I*  1 

SIND0042 

DO  16  J=IP1,NS 

SIND0043 

IE (KEY (I)  .EQ. 1. AND. KEY  (J) .EQ. 1)  NK=1 

SIND0044 

IF (KEY (I) .EQ. 1. AND. KEY  (J) .EQ. 2)  NK=2 

SIND0045 

IF (KEY (I) .EQ. 2. AND. KEY  (J) .EQ.  1)  NK=3 

SIN  D0046 

IF (KEY (I)  .EQ. 2. AND. KEY  (JJ.EQ.2)  NK  =  4 

SIND004  7 

GO  TO  (11,12,13,14) ,  NK 

SIND0048 

11 

CONTINUE 

SIN  D0049 

IND  (I* NO, J*NU)=LH1 (YS (I) , YS (J) , ZS (I ) , ZS (J) ) 

SIND0050 

GO  10  -15 

SIND  005 1 

12 

IND(I  +  NU,J*NU)  =LM1  ( YS  ( I)  ,YP(J),ZS(I)  ,ZP(J)) 

SIND0052 

GC  TC  15 

SINDOO  5  3 

13 

IND (I»NU, J»N  U) =LM 1  (YP  (I) , YS  ( J) , ZP  (I ) , ZS (J) ) 

SIND0054 
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• 

GO  10  15 

SIND0055 

14  IND  (I*NU,J*NU)  =LN1(YP{I)  ,YP(J)  ,ZP(I)  ,ZP(J)) 

SIND0056 

• 

.  - 

15  CONTINUE 

SIND0057 

IND(I*NUPS,J*NUPS)  =  LH 1  (YS(I) ,YS(J),ZS(I)  ,ZS(J)) 

3IND0058 

7*. 

16  CCN1INUE 

SIND0059 

' 

IP  (NO. EC. 0)  GO  TO  20 

SIND0060 

DO  19  1=1, HU 

SINDOObl 

DC  19  J= 1 , NS 

SIND0Q62 

IF  (KEY  (J) .EQ.2)  GO  TO  17 

SIND0063 

• 

IND  (I,  J»NU)  =1.111  (YH  (I)  ,YS  (J)  ,ZH  (I)  ,ZS  ( J)  ) 

SINB0064 

GC  TC  18 

SIND0065 

17  CONTINUE 

SIND0066 

IND  (I, J+NU) =  LN 1 ( YH  (I) ,YP (J) ,ZU(I) ,ZP (J) ) 

SINDC067 

r. 

18  CONTINUE 

3IND0066 

IND (I, NUPS) =LN1 (Yfl (I) , YS ( J) ,ZH  ( I) , ZS  (J) ) 

SIND0069 

19  CONTINUE 

SINDC07C 

- * 

20  CONTINUE 

SIND007 1 

21  CONTINUE 

SI ND  0072 

IF (NS, EQ. 0)  GO  TO  28 

SIN  D0073 

DC  27  1=1, NS 

SIND0074 

IP  (KEY  (I) .EQ.2)  GO  TO  24 

SIND0075 

IND(I»NU,I*NUPS)=LS1  (RS(I)  ,  YS  Cl)  ) 

SIND0076 

DC  23  J= 1 , NS 

SIND0077 

IF(I.EQ.J)  GC  TO  22 

SIND0078 

fe  ' 

IND(I»NU,J*NUPS)=LN1  (YS(I) , YS (J) ,ZS(I)  ,ZS(J) ) 

SIND0079 

22  CONTINUE 

SIND0080 

23  CONTINUE 

SIND008  1 

GC  TO  26 

S IN  DC  082 

24  CONTINUE 

SIND0083 

DC  25  J=1 , NS 

SIND0084 

25  IND(I*NU,J*NUPS)  =LH  1  (YP  (I)  ,YS(J)  ,ZP(I)  ,ZS(J)  ) 

SIND0085 

m  ■ 

26  CONTINUE 

SIND0086 

W 

27  CONTINUE 

SIND0087 

28  CONTINUE 

SIND0088 

GO  TC  57 

SIND0C89 

,  ■> 

29  CONTINUE 

SINDC090 

*  '  ' 

I F ( NU. EQ. 0)  GO  TO  31 

SINDC091 

DO  30  1=1, NU 

S  TND0092 

30  IND  (I,I)=IS2  (RH(I)  ,YW(I)  ) 

SIND0093 

i 

31  CONTINUE 

SINDC094 

■  i 

IF ( MS. EQ. 0)  GC  TO  36 

SIND0J95 

DC  35  1=1, NS 

SI ND  0096 

TKEY=KEY (I) 

S1NDC097 

GC  TC  (32,33) ,  IKEY 

SIHD0098 

3  2  IND  (I*NU,I*NU)=LS2 (RS(I)  ,YS  (I)) 

SINDC099 

GO  TO  34 

•SIND0100 

-  *• 

3  3  IND (I»NU,I»NU) =LS2  (BNP  (I)  ,YP  (I) ) 

SIND0101 

34  CONTINUE 

SIND0102 

IND(I*NUPS,I*NUPS)  =LS2  (HUH  (I)  ,  Y  S  ( I)  ) 

SINDC103 

35  CONTINUE 

SIND0104 

36  CONTINUE 

SI ND  0 1 05 

IF (NO. EQ. 0)  GO  TO  38 

SIND01 Ob 

NUN  1  =  NU- 1 

SINDO 107 

DC  37  1=1, NUN1 

SI ND  0 1 08 

IP1=I*  1 

3IND01C9 

DO  37  J=IP 1 , NU 

STND01  10 

37  IND  (I, J)=LH2  ( Y W  (I ) , YU  ( J) ,ZH  (I) ,ZN  (J) ) 

SIND  0  1  1  1 

38  CONTINUE 

SIND01  12 

IF  ( NS. EQ. 0)  GO  TO  49 

SINDO 11 3 

NSN 1 =NS-1 

SIND01 1 4 

DO  44  1=1, NSN 1 

SIND01 15 

* 

138 

1 

1 

: J 

1 1 1  =  1  +  1 

DC  44  J=I PI ,  NS 

IF (KEY (I) .  EQ. 1. AMD. KEY  (J) .  EQ. 1)  NK=1 
IF  (KEY  (I)  .  EQ.  1.  AMD.  KEY  (J)  .EQ.  2)  NK=2 
IF (KEY (I) .EQ.2. AND. KEY  (J) . EQ.  1)  NK=3 
IF (KEY (I) .EQ. 2. AND. KEY  (J).EQ. 2)  NK=4 
GC  TC  (39,40,41,42) ,  NK 
39  CONTINUE 

IND  (I»NU, J*NU)=LM2  (YS  (I)  ,YS  (J)  ,2S  (I)  ,ZS  (J)  ) 

GO  TC  43 

4  0  IND  (1+  NU, J+  NU) =  LH2 ( YS ( I)  , YP  ( J) , ZS  (I) , ZP  (J ) ) 

GC  TO  43 

4  1  IND  (I*  NO,  J*NU)  =LN2  (  YP  (  I)  ,YS  (J)  ,ZP  (I)  ,ZS(J)  ) 

GC  TC  43 

42  IND  (I*NU,J*NU)=LM2 (YP (I) ,YP (J) ,ZP (I)  ,ZP(J)  ) 

43  CCNTINOE 

IND  (I*NUPS,J*NUPS) =  LM2  (YS  (I)  ,  YS  ( J)  ,  ZS  (I)  ,  ZS  (J  )  ) 

44  CCNTINOE 

IF ( NU. EQ. 0)  GO  TO  48 
DC  47  1=1, NU 
DO  47  J= 1 , NS 

IF  (KEY  (J)  •  EQ .  2)  GO  TO  45 

IND  (I, J  +  NU) =  LM2 (YH  (I)  ,YS (J)  ,Z« (I) ,ZS (J)  ) 

GC  TO  46 

45  CCNTINOE 

IND  (I,  J+  NU) =LH2  (YH  (I)  , YP  (J)  ,ZH (I) ,ZP  (J) ) 

46  CONTINUE 

IND (I, J*NUPS)  =  LH2(YH  (I)  , YS (J)  ,ZH (I) ,ZS(0)) 

47  CCNTINUE 

48  CONTINUE 

49  CCNTINUE 
IF(SS.EQ.O)  GC  TO  56 
DO  55  1=1, NS 

IF  (KEY  (I). EQ. 2)  GO  TO  52 

IND  (1+  NU, 1  +  NUPS) =LS2 (BS (I)  ,YS  (I) ) 

DO  51  J=1 ,NS 

IF ( I. EQ. J)  GO  TO  50 

IND (I*  NO, J+NUPS) =  LS2 ( YS ( I) , YS (J) , ZS  (I) ,ZS  (J) ) 

50  CONTINUE 

51  CONTINUE 
GC  TC  54 

52  CONTINUE 

DC  53  J=1,NS 

53  IND  (1+  NU,  J*  NUPS)  =  LN2 ( YP (I)  , YS (J) , ZP  (I) ,ZS  (J) ) 

54  CONTINUE 

55  CONTINUE 

56  CCNTINUE 

57  CONTINUE 

DO  58  1=2, NT 
I N 1 =1- 1 
DO  58  J= 1 , IN  1 

58  IND  (I,J)  =IND(J,  I) 

RETURN 

END 


SUBROUTINE  PHI  FOR  COMPUTING  THE  CHAIN  PARAMETER  MATRICES 


SUBROUTINE  PHI  (PH IT, NT  ,T  NT  ,  LH  IN  ,G  Afl,  Y ,  Y I N  V,  T,  TT ,  EP ,  EN) 
IMPLICIT  BEAL*8  (A-H.O-Z) 


S IN  DO  1  16 
SIND01  17 
3IND01  18 
SIND01  19 
SIND0120 
SIND0121 
S IN  DO  1 22 
SIND0123 
SIND0124 
3IND0125 
SIND0126 
SINDC1 27 
SIND0128 
SINDO 129 
SINDC130 
SINDO 131 
SINDO  132 
SIND01 33 
SIND0134 
SINDO  135 
SIND01 36 
SINDO 1 37 
SINDO  1  38 
SIND01 39 
SIND0140 
SIND0141 
SI NDO 1 42 
SIND0143 
SINDO  144 
SIND01 45 
3IND0146 
SINDO 147 
SINDC148 
3IND0149 
SIND0150 
SIND01 5 1 
SIND0152 
SINDO 153 
SI ND  C 1 54 
SIND0155 
SINDO 1 56 
SINDC157 
SIND0158 
SINDO  159 
SI ND  C 1 6  0 
5IND0161 
SINDO  162 
SIND0163 
SINDO 164 
SINDO 165 
SI ND  01 66 
S IN  DO  1 6  7 
SINDO 168 
SI ND  0 1 6  9 
S  PH  1000  1 
SPHI0002 
SPH10003 
S  PH  10004 
SPH 10005 
SPH I  0006 
SPH 10007 


Ti 

INTEGER  TNT 

spuioooe 

jj 

PEA  1*3  LfllN 

5PHI00C9 

CCN  ELEX* 1 6  PHIT  (TNT, TNT)  , GAN  (NT)  , Y  (Nl.NT)  ,  YIN  V (NT, NT)  ,T  (NT, NT)  , 

SPHI001 0 

*/;  - 

1TI ( NT, NT)  ,  EP (NT)  ,  EN (NT) 

3PHI00 1 1 

CCM  EIE  X* 1 6  ZEROC,ONEC,XJ,JCMEGA,SUNC,SUNC1,SUNC2 

SPHIC012 

CCHNGN  /CCON/  ZEROC,ON EC, XJ , JONEGA 

3  PHI  00 1 3 

T  NO  =  2.  DO 

SPHI00 14 

s 

DC  1  1=1, NT 

SPH 10015 

EP(  I)  =  CDEXP  (GAN  (I)  *LM  IN) 

S  PH  100 1 0 

i 

£N(I)=CDEXP(-GAN(I) *LNIN) 

SPIIIOO  17 

, 

DC  3  1=1, NT 

SPHIC01 8 

DC  3  J=  1 , NT 

3  PH  100 1 9 

“ 

SOHC=ZERCC 

SPH IOO  20 

SUHC 1=ZE80C 

SPH I 002 1 

SUHC  2=Z  EROC 

S  PHI0022 

DC  2  K= 1 , NT 

SPHI0023 

,i 

SUHC1=SUHC1+T  (I,  K)  *  (EP  (K)  +  EN  (K)  )  *TI  (K,  J)  /THO 

SPHI0024 

'-V- 

3UBC2=SOHC2*T(I,K)  *GAN  (K)  *  (EP  (K)-EN  (F) )  *TI  (K,  J)  /TWO 

SPHI0025 

2 

SUNC=SUHC>T(I,K) ♦ ( EP (K) -EN ( K) ) *TI (K , J) /  (T HO*G AH (K)  ) 

S  PH  100  26 

PHIT  (I , J) =  SU  HC 1 

SPHI0027 

PPIT (I, J»NT)  =— SONC2 

SPHI0028 

PHIT (I ♦ NT, J) =-SUMC 

SPHI0029 

"* 

3 

PHIT (I*  NT, J*  NT) =SUHC 1 

SPHI 0030 

DC  5  1=1, NT 

SPHI0031 

DO  5  J= 1 , NT 

SPHI0032 

SUNC=ZEBOC 

SPHI 0033 

SUHC 1=ZEROC 

SPHI0034 

DC  4  K= 1, ST 

SPH IOO  35 

*  *’ 

SUHC=SUNC*YINV (I,K)*PHIT  (K,J) 

SPHI 0036 

N-.‘ 

4 

SUHC 1= SUHC 1* YIN V (I,K) *PHIT  ( K , J*NT) 

SPHI0037 

m 

T  (I  ,  J)  =SUHC 

SPHI0038 

5 

TI (I,J) =SUNC1 

SPHI 0039 

>*;"  ■ 

DC  7  1=1, NT 

SPHI0040 

DO  7  J= 1 , NT 

SPHI004 1 

■  .  v 

SUNC=ZEROC 

SPHI0042 

SUHC 1=ZEBOC 

SPHI0043 

DO  6  K= 1 , NT 

SPHI0044 

SUNC=SUNC*T (I , K) *Y  (K, J) 

SPHI0045 

B 

6 

SUHC 1=SUHC 1 ♦PHIT (I*NT,K) *Y  (K,J) 

SPHI0046 

PHIT  (I  ,  J)  =SU  HC 

SPHI0047 

7 

YINV (I,J) =SUHC  1 

SPHI0048 

DC  8  1=1, NT 

5PHI0049 

/ 

DC  8  J= 1 , NT 

SPHI0050 

V*. 

PHIT (I , J*NT) =TI (I, J) 

SPH 1005 1 

8 

PHIT (I^NT.J)  =YINV (I, J) 

SPHI0052 

■■m 

BETUBN 

SPHI0053 

END 

SPHI 0054 

c 

SADH0001 

c 

SADH0002 

SUBBCUTINE  ADH ACC  FOR  ACOING  ELASTANCES  TO  THE  INVERSE 

SADHC003 

c 

OF  THE  CAPACITANCE  HATRIX 

S ADH0004 

c 

SADH0005 

:.m  c 

SADH0006 

1  m 

SUBROUTINE  A CN A  CD (KE Y, NU , NS ,NT, NUPS , HS , TS , RWH  ,ERS , SS  ,5T V  ,NNU,NNS) 

S ADH0007 

IEPIICIT  BEAL*8  (A-H.C-Z) 

SADN0008 

■ 

TNTICEB  KEY  (NNS) 

SADHC009 

REAl*8  HU,HU02PI,aU04PI 

SADH0010 

-  ■  . 

REA  1*8  RS  (NNS) ,TS  (NNS)  ,BNH(NNS)  , ERS  (NNS) , SS  (NT, NT) ,STV (NNS) 

SADHOO 1 1 

COHHCN  /BCON/  ZEBO, ONE  , THO, THREE, FOUR, HU, NU02PI,HU04PI, 

SADH0012 

f- 

1PI,BADEG,ERTE, SIGCOP , P , F , V  2 

SADHOO 1 3 

• 

IF(NS.EC.O)  GO  TO  4 

SADHOO  14 

140 

• 

r- 

* 

iMHMM 1 

DC  3  1=1, NS 

IF  (KEY  (I)  •  EQ  .  2)  GO  TO  1 
STB  A  N=ST V  (I) 

SBB=DLCG  ( (BS (I) -T3 (I) )/RBH (I) ) / (TBO* El* EBS  (I) *EBTE) 
GC  TC  2 

1  CONTINUE 
STR  A  N=ZEBO 
S8B=ZER0 

2  CONTINUE 

SS(I*N0,I*NUPS)=SS (I*NU,I+NUPS)  +  STB A  N 
SS(I*NUPS,I«-NU)  =SS(I*NUPS,I»NU)  ♦ SIB AN 
SS(I+NUPS,I*NUPS) =SS(I*NUPS,I*NOPS)  ♦SBB+STRAN 

3  CONTINUE 

4  CONTINUE 
BETUBN 
END 


SUBROUTINE  INPACE  FOR  ADDING  INPEIANCES  TO  THE  INPEEANCE 
3  ATE  IX 


SUBROUTINE  IMP ADD  (Z , KE Y , NU, NS, NT, NUES , ZS V ,ZD V , LT V , ZBPLV , ZBBHV , 

1 NNU  ,  NNS) 

I FPIICIT  BEAL*8  <A-H,0-Z) 

REAI*8  NU,HU02PI,NU04PI, LTB AN, LTB 
CCNFIEX*16  Z ERCC, ON EC, XJ , JON EGA 
INTEGER  KEY (NNS) 

REAI*8  LTV (NNS) 

CCNFIEX* 16  ZSELF,ZDIF,Z(NT,NT)  ,Z5V(NNS)  ,ZDV(NNS) ,ZWFIV(NNS)  , 
1ZBBHV  (NNS) 

CCHHCN  /RCON/  ZERO ,ONE  , TWO, THREE , FOUR, HU, BU02PI , HU04 PI , 
1PI,BADEG , EBTE, SIGCOP , B , F , V2 
CONNCN  /CCO N/  ZEBOC,ONEC,XJ, JONEGA 
IF (NS. EQ. 0)  GO  TO  4 
DO  3  1=1, NS 

IF  (KEY  (I) . EQ. 2)  GO  TO  1 
ZSEIF=ZSV  (I) 

ZDI F=ZDV (I) 

LTRAN=LTV  (I) 

GO  TO  2 

1  CONTINUE 

ZSE I F=ZBPLV (I) 

Z  El F=Z  SELF 
LTB A  N=ZERO 

2  CONTINUE 

Z(I*NU,I*NU)  =Z(I*NU,I*NU) »ZSELF 

Z  (I*NU,I*NUPS)=Z(I+NO,I*NUPS)  ♦ZSELF-2EIF- JONEGA* LTR AN 
Z  (I ♦ NUPS, It NU) =Z  (I* NUPS, I  +  NU) *ZS El F-ZEIF-JON E G A*LTP AN 
Z ( I *NU PS, I*NUPS) =Z  (I*NUPS,I*N0PS) *TBC*  (ZSELF- ZDIF-JCMEG A*LT R AN) 
1 *ZB BHV  (I) 

3  CONTINUE 

4  CONTINUE 
BETUBN 
END 


SADS0015 
SADN0016 
5ADM00 17 
SADH0018 
SADN0019 
SADN0020 
SADN0021 
S AON 00  22 
SADN0023 
SADHC024 
S ADH0025 
SADN0026 
SADH0027 
SADN0028 
SADH0029 
SADN0030 
SINP0001 
SINP0002 
SINP0003 
SIMP0004 
SIHP0005 
SI NP 00 06 
SIMP0007 
SIRP0008 
SIHP0009 
SIN  POO  1 0 
SIMPOO  1 1 
SIHP00'2 
SIM  POO  1 3 
SINPOO  14 
SIHP001 5 
SIN  PO0 1 6 
SINPOO 17 
SIHPC018 
S IN  PO0 1 9 
SINP0020 
SINP002 1 
SIHP0022 
SIHP0023 
SINP0024 
SINP0025 
3INP0026 
SISP0027 
SINP0028 
SIHP0029 
SINP0030 
SIMP0031 
SIM  POO  32 
SINP0033 
SIMP0034 
SIN  POO  35 
SIHP0036 
SINP0037 
SIMPOO  38 
SINP0039 


Conversion  of  SHIELD  to  Single  Precision 


I 


i 


Delete  0038 


Card 

Double  Precision 

Single  Precision 

0039 

REAL*8 

REAL 

0041 

COMPLEX* 16 

COMPLEX 

0063 

REAL* 8 

REAL 

0075 

COMPLEX* 16 

COMPLEX 

0090 

REAL *8 

REAL 

0091 

COMPLEX* 16 

COMPLEX 

0103 

REAL* 8 

REAL 

0108 

COMPLEX*16 

COMPLEX 

0117 

0.0D0 

O.OEO 

0118 

1.0DO 

l.OEO 

0119 

2.0D0 

2.0EO 

0120 

3.0D0 

3.0EO 

0121 

4.0DO 

4.0EO 

0122 

2.54D-5 

2.54E-5 

0123 

2.D-7 

2.E-7 

0124 

l.D-7 

l.E-7 

0125 

180 . DO 

180. E 0 

0126 

2.997925D8 

2.997925E8 

0127 

5.8D-7 

5.8E-7 

0128 

DCMPLX 

CMPLX 

0129 

DCMPLX 

CMPLX 

0130 

DCMPLX 

CMPLX 

0131 

DATAN 

ATAN 

_  0 


m- 
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Card 


Double  Precision 


Single  Precision 


f  1  '  '  •  • 

l, 

T-: 

:• 


[W 

0241 

45.  DO 

45. EO 

•  ■ 

0345 

DSIN 

SIN 

0346 

DCOS 

COS 

m 

vj 

0347 

DSIN 

SIN 

V  ' 

.  - 

0348 

DCOS 

COS 

0351 

DSQRT 

SQRT 

* 

0351 

DCOS 

COS 

A.  < 

• 

0352 

DATAN2 

ATAN2 

0352 

DSIN 

SIN 

r* 

0353 

DCOS 

COS 

ft 

- 

0354 

DSQRT 

SQRT 

i 

0354 

DCOS 

COS 

0355 

DATAN2 

ATAN2 

r 

- - 

9 

0355 

DSIN 

SIN 

•.  ■ 

i 

0356 

DCOS 

COS 

V  1 

0530 

DCMPLX 

CMP  LX 

V 

w-.f 

0609 

DREAL 

REAL 

■ 

0620 

CDSQRT 

CSQRT 

£ 

0680 

DREAL 

REAL 

0691 

CDSQRT 

CSQRT 

.  ■ 

0735 

DREAL 

REAL 

• 

0746 

CDSQRT 

CSQRT 

;• 

0774 

DREAL 

REAL 

■ 

0785 

CDSQRT 

CSQRT 

; 

i 

0931 

CDABS 

CABS 

j 

• 
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Card 


Double  Precision 


i 


Single  Precision 


0932 

DREAL 

REAL 

0933 

DIMAG 

AIMAG 

0934 

DATAN2 

AT  AN  2 

0936 

CDABS 

CABS 

0937 

DREAL 

REAL 

0938 

DIMAG 

AIMAG 

0939 

DATAN2 

ATAN2 

0950 

CDABS 

CABS 

0951 

DREAL 

REAL 

0952 

DIMAG 

AIMAG 

0953 

DATAN2 

ATAN2 

0955 

CDABS 

CABS 

0956 

DREAL 

REAL 

0957 

DIMAG 

AIMAG 

0958 

DATAN2 

ATAN2 

FUNCTION  LSI 

0007 

LS1*8 

LSI 

Delete  0008 

0009 

REAL* 8 

REAL 

0012 

DLOG 

ALOG 

FUNCTION  LM1 

0007 

LM1*8 

LM1 

Delete  0008 

0009 

REAL* 8 

REAL 

0015 

DLOG 

ALOG 
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FUNCTION  LS2 


:* 


7* 

L 

0007 

LS2*8 

LS2 

Delete  0008 

0009 

REAL*8 

REAL 

V 

H 

0012 

DLOG 

ALOG 

FUNCTION  LM2 

0007 

LM2*8 

LM2 

3 

Delete  0008 

0009 

REAL *8 

REAL 

1 

* 

1 

0016 

DLOG 

ALOG 

0017 

DCOS 

COS 

"  V 

* 

0017 

DCOS 

COS 

i 

FUNCTION  STB 

0007 

STB*8 

STB 

•  ■ 

v 

Delete  0008 

r> 

t~ . 

1 

0010 

REAL* 4 

REAL 

1 

0011 

REAL*8 

REAL 

's 

0020 

DCOS 

COS 

r— 

0023 

1.0D-2 

1.0E-2 

>  * 

0024 

5.0D0 

5.0E0 

0026 

5.0D-2 

5.0E-2 

* 

0045 

45.0DO 

45.0EO 

M 


B 
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FUNCTION  LTB 

0007 

LTB*8 

LTB 

Delete  0008 

0010 

REAL*4 

REAL 

0011 

REAL*8 

REAL 

0020 

DCOS 

COS 

0023 

1.0D-2 

1.0E-2 

0024 

5.0DO 

5.0EO 

0026 

5.0D-2 

5.0E-2 

0045 

45.0D0 

45.0EO 

FUNCTION  ZWW 

0006 

ZWW* 16 

ZWW 

Delete  0007 

0008 

REAL*8 

REAL 

0009 

COMPLEX* 16 

COMPLEX 

0010 

Change  all  D's  to  E’s 

0014 

DSQRT 

SQRT 

FUNCTION  ZDB 

0007 

ZDB*16 

ZDB 

Delete  0008 

0010 

REAL*8 

REAL 

0011 

C0MPLEX*16 

COMPLEX 

mwi; 


itIKi] 


3 


function  zsb 


i 

I 

1 

I  - 

*.i 


i 


0006 

ZSB*16 

ZSB 

Delete  0007 

0009 

REAL*8 

REAL 

0010 

C0MPLEX*16 

COMPLEX 

0014 

Change  all  D's  to  E's 

0015 

DSQRT 

SQRT 

0021 

DCOS 

COS 

0023 

DCOS 

COS 

0026 

DCOS 

COS 

FUNCTION  ZDS 

0007 

ZDS*16 

ZDS 

Delete  0008 

0009 

REAL*8 

REAL 

0010 

COMPLEX* 16 

COMPLEX 

0014 

DSQRT 

SQRT 

0017 

CDEXP 

CEXP 

0017 

CDEXP 

CEXP 

FUNCTION  ZSS 

0007 

ZSS*16 

ZSS 

Delete  0008 

0009 

REAL*8 

REAL 

0010 

COMPLEX* 16 

COMPLEX 

0014 

DSQRT 

SQRT 

0016 

DEXP 

EXP 

BE 


L6 


H 


& 

FUNCTION  ZSS  (Cont’d) 

sa 

0017 

DEXP 

EXP 

0017 

DEXP 

EXP 

r\\ 

t  . 

r.  •  "• 

0018 

DSIN 

SIN 

0018 

DCOS 

COS 

0020 

DSIN 

SIN 

»„  ■ 

j 

0020 

DCOS 

COS 

s 

»*. 

SUBROUTINE  MULTC 

/ 

Delete  0007 

'*% 

V 

0008 

COMPLEX*16 

COMPLEX 

. 

SUBROUTINE  SCAP 

Delete  0008 

1 

0010 

REAL*8 

REAL 

0011 

2.997925D8 

2.997925E8 

■ 

Delete  0008 

SUBROUTINE  INDUCT 

0010 

REAL*8 

REAL 

0011 

REAL* 8 

REAL 

SUBROUTINE  PHI 

Delete  0007 

0009 

REAL* 8 

REAL 

0010 

COMPLEX* 16 

COMPLEX 

0012 

COMPLEX* 16 

COMPLEX 

GH 


L4 


2 


lili] 


K 


SUBROUTINE  PHI  (Cont’d) 


0016 

CDEXP 

CEXP 

• 

0017 

CDEXP 

CEXP 

SUBROUTINE  ADMADD 

■ 

Delete  0008 

0010 

REAL* 8 

REAL 

0011 

REAL*8 

REAL 

'«■  '  « 

>• 

0018 

DLOG 

ALOG 

SUBROUTINE  IMPADD 

Delete  0009 

* 

0010 

REAL*8 

REAL 

0011 

C0MPLEX*16 

COMPLEX 

* 

0013 

REAL*8 

REAL 

' 

0014 

COMPLEX* 16 

COMPLEX 

ij 
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APPENDIX  C 


Flowchart  of  SHIELD 


Read  and  print  terminal  characteristics  data 
(including  shield  grounding  data) 


Determine  shortest  pigtail,  LMINR,  and  set  entries 
in  array  KEY  to  indicate  for  each  cable 
whether  a  shield  (1)  or  a  pigtail  (2) 
is  present  over  LMINR 


Compute  chain  parameter  matrix  of  this  right  pigtail 
section  of  length  LMINR  and  store  in  PHIT 


Multiply  this  section  chain  parameter  matrix  PHIT, 
and  the  previous,  accumulated  chain  parameter  matrix, 
PHIR,  as  TPHI  =  PHIR*PHIT  and  restore  in  PHIR 


Subtract  LMINR  from  each  previous  pigtail  length 
and  restore  in  array  LPT 


GO  TO  117 


Store  left  pigtail  lengths  in  LPT  and  initialize  chain 
parameter  matrix  of  left  pigtail  section, 

PHIL,  to  the  identity  matrix 


Determine  longest  pigtail, 
LMAXL,  in  array  LPT 


LMAXL=0 


LMAXL# 0 


(0713) 
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MISSION 

of 

Rome  Air  Development  Center 

RA VC  plan*  and  ex.ec utes  AeAe.aA.ch,  development,  test  and 
selected  acquisition  programs  in  support  o f  Command,  Cont/iol 
Communications  and  Intelligence  (C^I)  activities.  Technical 
and  engineering  support  within  areas  o f  technical  competence 
is  provided  to  ESP  Program  Offices  (PCs)  and  other  ESV 
elements.  The  principal  technical  mission  areas  are 
communications ,  electromagnetic  guidance  and  control,  sur¬ 
veillance  of  ground  and  aerospace  objects,  intelligence  data 
collection  and  handling,  information  system  technology , 
ionospheric  propagation,  solid  state  sciences,  microwave 
physics  and  electronic  reliability,  maintainability  and 
compatibility . 


